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ABSTRACT 

This  thesis  proposes  a  framework  by  which  optimal  system  renewal  decisions 
may  be  made  in  a  consistent  and  timely  manner  within  the  military,  under  realistic 
conditions  of  changing  environments.  A  unique  network  representation  of  the  system 
renewal  process  was  used  to  develop  a  prototype  version  of  the  analytical  core  model. 
Its  plausibility  and  usefulness  was  demonstrated  by  a  series  of  case  studies.  The  case 
studies  also  show  how  pertinent  staff  forecasts  can  be  organised  and  integrated  to 
provide  decision  makers  with  a  broad,  consistent,  long-term  perspective  of  the  issues 
relevant  to  system  renewal  planning.  They  are  presented  with  a  graphic  picture  of 
the  entire  solution  space  as  structured  by  the  scenarios  considered.  Various  solutions 
are  suggested  for  each  scenario  and  their  robustness  may  be  tested  by  thoroughly 
exercising  the  model  for  a  wide  range  of  scenarios.  Prediction  of  what  and  when 
renewals  are  likely  to  be  and  estimation  of  the  associated  cost  and  effectiveness 
allow  anticipatory  long-term  plans  to  be  formulated.  The  thesis  also  suggests  how 
bulk  versus  phased  procurement  decisions  and  force  level  and  mix  issues  could  be 
analysed  using  the  model. 
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I.  INTRODUCTION 

A.  INTRODUCTION 

Military  forces  continually  undergo  change  as  obsolete  or  uneconomical  assets 
are  retired,  new  weapons  are  introduced  and  existing  ones  modified  or  upgraded. 
This  process  of  system  renewal  is  affected  by  a  myriad  of  economic,  technological 
and  politico-military  considerations.  Changes  in  enemy  capabilities  or  doctrine,  for 
instance,  may  render  certain  systems  less  effective  or  even  obsolete.  New  technology 
may  offer  opportunities  to  improve  present  capabilities.  Most  of  the  available  open 
literature,  analyses  and  models  for  system  renewal  are  either  one-dimensional,  too 
general,  or  too  commercially  oriented  to  be  useful  to  military  planners.  There  is  a 
need  to  synthesize  a  coherent,  balanced  view  relevant  to  military  planners  to  provide 
comprehensive  policy  guidance  on  upgrading,  replacement  and  related  matters. 

B.  PROBLEM 

The  system  renewal  planning  problem  is  encountered  prior  to  the  last  milestone 
of  the  defence  acquisition  cycle.  It  is  basically  one  of  being  able  to  evaluate  and 
decide  in  a  timely  and  rational  manner,  under  conditions  of  changing  technological, 
economic  and  politico-military  environments,  both  the  following: 

1.  When  will  a  defence  asset,  in  its  present  form,  become  obsolete  or  no  longer 
cost-effective  to  retain  (relative  to  existing  or  emerging  alternatives)? 

2.  What  should  then  be  done?  What  are  the  alternatives  and  which  one  should 
be  chosen? 


C.  DEFINITION 

System  renewal  encompasses  all  changes  made  to  maintain  or  improve  the 
cost-effectiveness  of  the  present  defence  system.  It  includes  the  acquisition  of  new 
capabilities,  retirement  from  service  of  obsolete  ones,  relegation  to  reserves  of  sur- 
plus capabiHties  and  conversion  of  systems  to  new  roles.  This  view  streamlines  the 
treatment  of  changes  to  the  status  quo  and  involves  a  higher  level  of  system  repre- 
sentation. It  allows  us  to  consider  complex  changes,  such  as  when  a  new  support 
ship  displaces  a  couple  of  smaller  ones,  one  of  which  is  preserved  for  the  boneyard 
and  the  other  converted  to  a  hospital  ship. 

D.  PURPOSE 

The  purpose  of  this  thesis  is  to  develop  a  rational  and  consistent  decision 
model  for  system  renewal  planning  that  is  relevant  and  useful  to  military  planners. 
Various  issues  pertinent  to  system  renewal  are  integrated  into  a  common  analytical 
framework  to  make  for  a  more  comprehensive  and  better  formalised  concept  of 
obsolescence  in  military  systems. 

This  work  is  intended  as  a  background  study  for  later  development  of  a  com- 
prehensive policy  to  complement  life-cycle  management  initiatives  in  the  Singapore 
Ministry  of  Defence.  It  serves  as  a  rudimentary  framework  that  can  be  expanded 
later  for  actual  use. 

E.  IMPETUS 

The  choice  of  this  research  topic  was  motivated  by  the  following: 

•  Singapore's  defence  buildup  has  been  underway  for  two  decades  and  already 
system  renewal  decisions  are  due  for  some  of  the  aging  systems  that  were  ac- 
quired initially.  A  framework  that  facilitates  taking  the  long-term  perspective 


in  system  renewal  planning  will  help  ensure  that  such  optimal  decisions  are 
made  in  a  systematic  and  timely  manner.  As  the  assets  involved  are  consid- 
erable, improvements  in  system  renewal  planning  is  of  significant  importance. 

•  Milestone  5  (upgrading/replacement)  was  introduced  in  U.S.  DoD  Instruc- 
tion 5000.2  in  September  1987.  No  clear  guidance  has  been  provided  on  the 
approaches  to  be  adopted  for  this  new  milestone. 

•  Available  literature  generally  either  deals  with  replacements  aimed  at  pre- 
empting failure  (preventive  maintenance)  or  reduces  the  problem  to  a  purely 
economic  one  (commercial  approach).  There  is  a  dearth  of  literature  of  rele- 
vance to  military  system  renewal  planning. 

F.     SCOPE 

Replacement  policies  concerning  routine  preventive  maintenance  are  not  ex- 
amined in  this  thesis.  The  scope  is  limited  to  major  systems  or  forces  for  which 
the  considerable  effort  involved  in  identifying,  tracking  and  analysing  capability  re- 
quirements, costs  and  measures  of  effectiveness  (MOEs)  are  justifiable.  At  the  heart 
of  the  framework  to  be  presented  is  an  analytical  model,  the  workings  of  which 
will  be  illustrated  under  various  assumptions  based  on  a  hypothetical  weapon  sys- 
tem. Attention  is  focused  on  single  mission  systems  in  changing  environments.  The 
treatment  of  threat  evolution,  force  level  planning  and  higher  level  optimisation 
are  briefly  discussed.  Simple  sensitivity  analysis  is  illustrated  but  full  and  proper 
stochastic  modelling  and  risk  analysis  lie  outside  the  scope  of  this  thesis. 


II.  ISSUES  AND  LITERATURE  SURVEY 

Systems  are  renewed  as  they  become  obsolete.  Why?  When  is  a  system 
considered  obsolete?  What  are  the  issues  involved  in  planning  for  system  renewal? 
These  are  the  larger  questions  addressed  in  this  chapter.  We  begin  by  explaining 
obsolescence  and  its  consequences.  Next,  for  a  broad  overview  of  the  system  renewal 
problem,  various  issues  pertaining  to  the  selection  and  scheduling  of  alternatives  are 
discussed.  The  chapter  ends  with  a  survey  of  some  of  the  research  that  have  been 
carried  out  on  system  renewal. 

A.      OBSOLESCENCE 

Obsolescence  may  be  categorised  as  follows: 

•  Economic:  A  system  is  considered  economically  obsolete  when  it  becomes  less 
cost-effective  than  available  alternatives  because  of  cost  escalation  attributed 
to  physical  deterioration,  scarcity  of  spares,  inflation,  etc.  An  old  truck  that 
breaks  down  often,  guzzles  fuel  that  has  become  expensive,  or  has  little  spares 
support  is  economically  obsolete. 

•  Technological:  This  form  of  obsolescence  occurs  when  innovation  makes  avail- 
able a  new,  more  cost-effective  alternative  to  an  incumbent  system.  Fighter 
planes  with  propellers,  for  example,  became  technologically  obsolete  when  jet- 
propelled  ones  became  available. 

•  Operational:  Operational  obsolescence  arises  when  a  system  fails  to  meet  its  re- 
quirements because  of  deterioration  in  performance,  enemy  counter-measures 
or  shifts  in  own  strategic/tactical  doctrine.     Immobile  defence  systems,  for 


instance,  may  become  obsolete  if  it  is  known  that  they  can  be  effectively  en- 
gaged or  bypassed  by  the  enemy,  or  if  strategic  emphasis  shifts  from  attrition 
to  maneuver  warfare. 

Obsolescence  may  be  a  compound  effect  of  the  above.  A  gun  that  is  no  longer 
accurate  may  also  be  uneconomical  to  maintain. 

The  continued  operation  of  an  obsolete  system  is  always  a  waste  of  money, 
but  technological  and  operational  obsolescence  may  have  more  serious  consequences. 
Failure  to  exploit  technological  opportunities  or  remedy  operational  deficiencies  may 
lead  to  military  disasters.  Technology  is  not  exploited  merely  by  keeping  pace  with 
research  as  was  shown  by  the  French  experience  in  1939.  The  French  Army  then 
had  for  some  time  developed  prototypes  of  aircraft,  armor  and  antitank  weapons 
superior  to  those  of  the  German  Army.  Their  wait  for  still  better  ones  to  come  along 
before  introducing  them  into  the  inventory  was  a  reason  for  their  unpreparedness 
[Possony  and  Pournelle,  1970].  The  Polish  Army's  fatal  cling  to  cavalry  is  a  clear 
lesson  in  the  consequence  of  operational  obsolescence. 

Not  all  obsolescence  is  bad.  The  term  "built-in  obsolescence"  is  sometimes 
used  for  the  practice  of  deliberately  designing  machines  that  deteriorate  rapidly 
after  a  given  life;  the  user  being  compensated  for  this  by  a  lower  purchase  price. 
This  can  be  seen,  for  example,  in  commercial  batteries  with  different  lives. 

B.     STRATEGIC  ISSUES 

Defence  forces  exist  for  a  purpose,  namely  to  deter  aggression  and,  if  deterence 
fails,  to  successfully  wage  war.  The  strategy  to  achieve  this  purpose  should  guide 
the  formulation  of  system  renewal  plans.  Potential  adversaries  should  be  identified 
and  their  intentions,  war-proneness,  capabilities,  weaknesses  and  strategies  assessed 


to  arrive  at  a  strategy  that  is  congruent  with  foreign  pohcy.  Competitive  strat- 
egy [Englund,  1987],  for  example,  seeks  a  favorable  cost  exchange  ratio  by  making 
obsolete  major  components  of  the  enemy's  defence  system.  Stealth  bombers  and 
SDI  are  manifestations  of  this  potentially  destabilising  strategy.  Some  references  on 
strategic  analysis  and  forecasting  are  provided  in  the  bibliography. 

There  are  many  facets  of  strategy  that  affects  system  renewal  plans.  Resources 
for  modernisation  may  be  channelled  to  boost  readiness  if  conflict  seems  imminent. 
The  danger  of  provoking  or  worsening  an  arms  race  must  be  weighed  carefully  if  plans 
call  for  a  significantly  increased  capability.  The  impact  of  a  plan  on  the  economy  may 
also  be  significant.  The  B-2  program  recommended  by  the  USAF,  for  example,  asks 
for  a  whopping  seventy-five  billion  dollars.  The  impact  of  modest  levels  of  defence 
spending  on  the  economy  is  usually  seen  to  be  beneficial.  Fredericksen's  [1989] 
analysis,  for  example,  showed  a  favorable  feedback  relationship  between  defence 
spending  and  economic  growth  in  Indonesia  for  the  period  1964  to  1985.  Economic 
issues  and  methodology  in  arms  race  analysis  are  surveyed  in  Leidy  and  Staiger 
[1985]. 

The  link  between  strategy  and  procurement  is  historically  lacking  in  the  US 
(and  Singapore).  The  chairman  of  a  Defense  Resource  Management  Study  in  1979, 
for  example,  declared  that  there  is  broad  agreement  that  the  first  P  in  PPBS  is 
silent  [Rice,  1979].  The  Packard  Commission  in  1986  stated: 


. . .  there  is  a  great  need  for  improvement  in  the  way  we  think  through  and  tie 
together  our  security  objectives,  what  we  spend  to  achieve  them,  and  what  we 
buy.  The  entire  undertaking  for  our  nation's  defence  requires  more  and  better 
long-range  planning.  This  will  involve  concerted  action  by  our  professional 
military,  the  civilian  leadership  of  the  Department  of  Defense,  the  President, 
and  the  Congress. 


Today,  there  is  no  rational  systenn  whereby  the  executive  branch  and  the 
Congress  reach  coherent  and  enduring  agreement  on  national  military  strategy, 
the  forces  to  carry  it  out,  and  the  funding  that  should  be  provided  -  in  light  of 
the  overall  economy  and  competing  claims  on  national  resources.  The  absence 
of  such  a  system  contributes  substantially  to  the  instability  and  uncertainty 
that  plagues  our  defence  program.  These  cause  imbalances  in  our  military 
forces  and  capabilities  and  increases  the  cost  of  procuring  military  equipment. 


This  thesis  seeks  to  improvise  a  link  between  strategy  and  force/ resource  plan- 
ning by  incorporating  strategic  outlook  into  system  renewal  plans.  The  cost  and 
effectiveness  of  alternatives  of  interest  are  projected  into  the  future  to  make  for 
more  forward-looking  plans.  Figure  2.1  illustrates  the  links  between  strategic  issues 
and  systenn  renewal  planning.  It  is  hoped  that  at  some  future  time  system  renewal 
plans  can  be  aggregated  and  linked  to  arms  race  models,  strategic  analysis  models, 
econometric  models,  etc.  to  provide  a  higher  level  of  optimisation. 

Technology  pushes  system  renewal  because  it  holds  the  key  to  one's  edge  over 
the  enemy.  However,  excessive  emphasis  on  performance  exacts  an  inordinate  cost. 
Figures  2.2  and  2.3  [Gansler,  1989]  show  the  cost  histories  and  performance  growth 
respectively  of  a  series  of  U.S.  fighter  and  attack  aircraft.  The  performance  assess- 
ment methodology  used  in  Figure  2.3  is  explained  in  Chapter  4.  As  can  be  seen, 
the  difference  in  unit  cost  between  high-end  aircraft  (e.g.  F-15A,  F-18)  and  low-end 
ones  (e.g.,  F-16A,  A-7A)  had  grown  significantly  through  the  years.  Performance 
growth  was  however  much  more  modest,  indicating  that  cost-effectiveness  of  high- 
end  aircraft  had  deteriorated  relative  to  that  of  the  low-end  ones.  Figure  2.3  suggests 
that  the  performance  of  the  first  operational  units  of  a  low-end  aircraft  lags  that  of 
high-end  ones  by  only  five  years.  It  also  implies  that  technological  superiority  over 
Soviet  counterparts  could  have  been  maintained  using  lower-risk,  lower-cost  aircraft. 
The  lower  cost  translates  into  more  aircraft,  which  may  matter  more  in  terms  of 


military  results.  Expensive  high-performance  or  multi-role  assets  also  require  in- 
tricate defences  a^  they  are  too  valuable  to  lose.  Quality  must  also  be  traded  off 
against  quantity  when  selecting  alternatives  in  a  system  renewal  plan.  This  tradeoff 
can  be  accomplished  within  the  model  to  be  introduced  in  the  following  chapters  by 
denoting  quantity  a^  a  performance  parameter.  The  impact  of  future  innovation  on 
the  cost-effectiveness  of  alternatives  of  interest  is  projected  so  that  optimal  numbers 
of  the  right  alternatives  will  be  procured  at  the  right  time. 

The  question  of  how  far  to  push  technology  in  a  procurement  program  has 
major  cost  and  effectiveness  implications.  Figure  2.4  illustrates  how  the  technolog- 
ical frontier  for  a  certain  type  of  system  may  advance  with  R&D  effort.  Point  A  is 
typical  of  a  target  set  to  ensure  the  ability  to  counter  probable  threats  over  a  long 
horizon.  Due  to  the  risks  involved,  it  is  likely  in  this  case  to  end  up  somewhere  inside 
the  area  bounded  by  the  points  A,  B  and  C.  The  less  ambitious  and  better  plan 
would  be  to  seek  to  quickly  and  cheaply  arrive  at  say  point  D  first,  then  improve 
incrementally  to  point  A.  Having  such  pre-planned  product  improvements  (P^I) 
requires  good  forecasts  of  likely  future  upgrades.  The  model  to  be  introduced  facil- 
itates the  selection  and  scheduling  of  such  upgrades.  The  reader  may  be  interested 
in  the  1984  report  by  the  U.S.  Defense  Science  Board  entitled  "Improved  Defense 
through  Equipment  Upgrades:  The  U.S.  and  Its  Security  Partners". 

Requirement  specifications  largely  determine  the  technological  content  of  a 
program.  These  requirements  may  not  be  the  result  of  careful  study.  In  fact,  they 
understandably  tend  to  reflect  the  psychological  comfort  levels  of  those  who  have  to 
deal  with  the  exigencies  of  war.  After  all,  it  never  hurts  to  be  able  to  go  the  extra 
mile.  Accordingly,  the  best  of  various  candidate  systems  often  becomes  the  "re- 
quirement". The  importance  of  specifying  objective,  mission-oriented  requirements 
cannot  be  overstressed. 
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Figure  2.1  :  System  Renewal  and  Related  Issues 
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C.     ECONOMIC  ISSUES 

Cost-efFectiveness  of  alternatives  are  dependent  on  future  interest  rate,  foreign 
exchange  rates,  inflation  rate,  resale/salvage/scrap  market  condition  and  budget 
allocations.  The  U.S.  Department  of  Defense  publishes  forecasts  of  inflation  indices 
for  various  categories  of  military  hardware  which  can  be  used  in  system  renewal 
planning.  DoD  would  of  course  like  to  estimate  conservatively  (high  inflation)  so 
as  to  get  more  money  to  play  safe  in  case  it  actually  is  high.  However,  in  times  of 
high  inflation,  there  is  in  fact  considerable  political  pressure  for  lower  estimates  in 
order  to  suppress  defence  budget  growth  projections.  The  cumulative  inflation  from 
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1978  to  1981  was  45  percent,  while  the  inflation  estimated  in  the  DoD  budget  was 
only  28  percent  [Gansler,  1989].  Therefore  one  must  exercise  discretion  in  using  the 
inflation  indices. 

Assuming  a  common  interest  rate  for  all  programs  does  not  make  for  equitable 
comparison  as  alternatives  with  difi^erent  life-spans  and  cost  profiles  are  affected 
differently.  One  should  try  to  forecast  the  rate  as  accurately  as  possible  to  avoid 
biased  results.  As  economic  factors  such  as  interest  rates,  foreign  exchange  rates 
and  market  conditions  are  very  difficult  to  estimate,  sensitivity  of  solutions  to  their 
likely  fluctuations  should  be  tested  in  actual  system  renewal  planning. 

Adoption  of  a  longer  term  perspective  in  system  renewal  planning  should  be 
matched  by  the  adoption  of  a  multi-year  budget.  Procurement  of  needed  quantities 
of  a  system  is  usually  most  economical  at  a  certain  rate.  An  optimal  system  renewal 
plan  may  involve  humps  in  budgetary  outlays.  The  minimum  efficient  rate  of  pro- 
duction for  F-15  aircraft,  for  example,  is  120  per  year.  However,  DoD  was  allowed 
only  41  per  year  [Gansler,  1989].  This  U.S.  practice  of  stretching  out  programs  to 
smooth  out  the  humps  tends  to  lead  to  gross  inefficiencies  in  procurement.  The 
floating  of  treasury  bonds,  for  example,  may  be  a  better  way  to  achieve  a  more 
efficient  procurement  program.  For  multi-year  budgets  to  work,  accurate  cost  esti- 
mates are  necessary.  Buying  ofl'-the-shelf  items,  adopting  low-risk  approaches  (e.g., 
by  emphasizing  P^/)  and  designing  to  cost  are  collaborative  means  towards  this 
end. 

D.     MODELLING  ISSUES 

Military  system  renewal  planners  have  to  grapple  with  problems  of  inulti- 
facetedness,  vagueness,  uncertainty  and  subjectivity.  Even  if  the  efl^ectiveness  of  a 
given  alternative  can  be  summarised  by  a  single  MOE,  one  has  still  to  consider  its 
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cost  and  probable  service  life.  The  system  renewal  problem  has  therefore  at  least 
three  dimensions  -  effectiveness,  cost  and  useful  life.  Often,  service  life  is  simply 
assumed  and  net  cost  is  amortised  into  a  rate  over  this  period  for  evaluations.  How- 
ever, advocates  of  an  alternative  will  tend  to  assume  a  longer  life  than  opponents. 
Instead  of  assuming  when  an  alternative  will  become  obsolete,  the  model  to  be  pre- 
sented derives  the  service  life  as  an  output  based  on  economic,  technological  and 
military  (requirement)  forecasts.  This  should  make  for  more  realistic  and  equitable 
evaluations. 

Risks  due  to  uncertainty  can  be  assessed  by  testing  the  sensitivity  of  solutions 
to  possibilities  learned  from  experience.  Dominance  or  robustness  of  a  solution 
would  provide  some  assurance  of  its  validity  under  a  range  of  expected  scenarios. 
The  quality  of  experience  is  of  critical  importance  for  such  risk  assessment.  As  such, 
qualified  experts  should  be  consulted  in  the  formulation  of  scenarios. 

Besides  multi-facetedness  and  uncertainty,  a  planner  has  to  contend  with  prob- 
lems of  vagueness.  Strategic  analysis  may  highlight  the  need  to  heighten  readiness  in 
the  coming  years  but  such  qualitative  goals  do  not  readily  lend  themselves  to  math- 
ematical modelling.  Even  if  readiness  can  be  precisely  quantified  by  a  parameter, 
there  may  not  be  consensus  concerning  what  constitutes  the  proper  level  of  readi- 
ness. Fuzzy  set  and  logic  theory  offer  means  to  bridge  the  gap  between  quantitative 
and  qualitative  reasoning,  thus  linking  strategy  formally  to  procurement  planning. 
The  efforts  of  J.  Dockery  [1985]  is  an  example  of  work  done  in  this  area.  Karprzyk 
and  Orlovski  [1987]  and  Zimmerman,  et  a!.,  [1984]  also  discuss  new  paradigms  more 
suited  to  situations  where  vagueness  and  uncertainty  are  clearly  present.  The  prob- 
lem of  vagueness  is  not  considered  in  this  thesis  but  may  be  addressed  in  future 
embellishments. 
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Although  the  exercise  of  human  judgment  is  essential  in  system  renewal  plan- 
ning, the  subjectivity  involved  must  be  managed  well  to  reduce  the  possibility  of  bias 
(conscious  or  otherwise).  The  problems  encountered  in  assessing  utility  functions, 
for  example,  were  surveyed  by  Hershey  et  al.,  [1982]  and  Farquhar  [1984].  It  should, 
however,  be  noted  that  intuitive  (subjective)  judgment  of  experts  may  in  certain 
cases  prove  superior  to  analytical  studies  [Hammond,  1987].  Human  behavioral  as- 
pects are  important  in  a  real  application,  but  they  are  not  considered  in  this  thesis 
because  the  cases  studied  are  hypothetical  ones  meant  for  illustrating  the  model 
that  was  developed.  Utility  functions,  for  example,  were  assumed  to  be  explicitly 
defined  even  though  such  functions  are  difficult  to  assess  accurately  in  practice. 

The  horizon  of  optimisation  is  an  important  issue  in  modelling  system  renewal. 
Very  often,  this  horizon  is  either  fixed  at  some  arbitrary  value  or  taken  to  be  in- 
finity. Both  practices  suffer  from  serious  flaws.  Choosing  horizons  that  are  a  small 
non-integral  number  of  life-spans  would  produce  incorrect  results.  This  occurs,  for 
instance,  if  the  horizon  for  a  system  with  a  typical  lifespan  of  ten  years  is  set  to, 
say,  twenty-five  years.  Assuming  an  infinite  horizon  is  clearly  unrealistic.  There  are 
procedures  that  allow  one  to  find  the  horizon  beyond  which  the  decision  for  the  first 
change  remains  optimal.  However,  these  so-called  "planning  horizon  procedures" 
are  essentially  also  based  on  the  infinite  horizon  assumption.  This  thesis  takes  a 
different  approach,  namely  the  finite-stage  approach  (from  dynamic  programming), 
in  which  the  optimisation  horizon  is  a  variable  output  given  by  the  years  spanned 
by  a  sequence  of  systems. 

Modelling  the  effectiveness  of  a  military  system  is  sometimes  confused  with 
predicting  the  outcome  of  its  use  in  battle.  The  two  are  related  but  not  synonmous. 
As  an  example,  let  us  suppose  that  a  certain  armour  unit  has  the  MOE  structure 
shown  in  Figure  2.5.     The  outcome  in  a  given  battle  scenario  may  be  described 
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by  the  number  of  own  and  enemy  losses.  Neither  of  these  numbers  expresses  fully 
the  effectiveness  of  the  force.  An  abstract  figure  of  merit  would  have  to  be  com- 
puted from  these  numbers  based  on  judgment  concerning  what  constitutes  success. 
Therefore,  one  might  say  that  battle  outcomes  are  predicted  by  combat  modelling 
whereas  effectiveness  is  evaluated  by  "combat  value"  modelling.  The  abstract  figure 
of  merit  may  be  said  to  be  the  MOE,  and  performance  parameters  that  determine 
this  value  may  be  said  to  be  the  measures  of  performance  (MOPs).  This  distinction 
is,  however,  an  artificial  one  as  the  MOPs  may  collectively  be  considered  to  be  a 
vector-valued  MOE.  Accordingly,  MOE  and  MOP  in  this  thesis  are  both  taken  to 
mean  the  same  thing  -  key  performance  parameters  that  determine  the  effective- 
ness of  a  system.  In  the  example  shown  in  Figure  2.5,  the  MOEs  would  be  range, 
accuracy  and  lethality  if  these  parameters  are  the  ones  that  have  major  effects  on 
effectiveness. 

E.     ORGANISATIONAL  ISSUES 

System  renewal  plans  may  be  developed  for  many  different  levels,  ranging  from 
a  single  equipment  to  an  entire  task  force.  Unless  they  are  coordinated  such  that  the 
objectives  of  upper  echelons  are  considered  by  lower  echelon  planners,  sub-optimal 
or  even  counterproductive  plans  are  likely  to  be  made.  Most  resource  allocation 
models  do  not  consider  the  organisation  as  a  factor  in  the  decision  process;  that 
is,  they  assume  a  monolithic  decision-making  entity.  The  same  decision-makers  are 
assumed  whether  the  subject  of  planning  is  a  single  artillery  gun  or  an  infantry  divi- 
sion. Ruefli  [1974]  surveyed  several  economic  and  management  science  models  that 
provide  remedy  for  this  unrealistic  assumption  by  explicitly  reflecting  organisational 
structure  in  the  model.  The  D an tzig- Wolfe  model  [Dantzig  and  Wolfe,  1960]  is  the 
best  known  of  the  classical  models.  In  more  recent  years,  models  such  as  the  one  by 
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Figure  2.5  :  Hierarchy  of  MOEs 

Source  :  Modified  version  of  the  diagram  from 
"Notes  on  MOEs"  by  E.B.Rockower,  Aug.  1985. 


Cassidy  et  al,  [1971]  that  are  based  on  behavioral  theories  have  emerged.  Organ- 
isational issues  are  not  considered  in  this  thesis  as  the  model  is  still  a  conceptual 
one  that  is  not  ready  for  organisation-wide  implementation. 

F.      LITERATURE  SURVEY 

No  hterature  suitable  for  military  system  renewal  planners  could  be  found  as 
all  those  surveyed  were  commercially  oriented  ones  that  subjugated  performance 
considerations  under  cost.    Most  of  them  consider  only  replacement  alternatives. 
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The  earliest  of  these  were  based  on  simple  calculations  and  differential  and  integral 
calculus.  Examples  are  Preinrich's  Constant  Chain  Model  [1940]  and  Terborgh's 
MAPI  system  [1949].  In  more  recent  years,  methods  such  as  dynamic  program- 
ming, Markov  chains  and  Pontryagin's  maximum  principle  have  become  more  pop- 
ular. The  first  such  model  was  introduced  by  R.  Bellman  [1955]  using  dynamic 
programming.  Howard  [1960]  used  a  policy-iteration  method  to  solve  multi-stage 
decision  problems  in  connection  with  Markov  processes.  Numerous  articles  conern- 
ing  stochastic  models  were  published  in  the  1960's.  A  survey  of  the  treatment  of 
stochastically  failing  equipment  was  compiled  by  Jorgenson  et  al.,  [1967].  Naslund 
[1966]  was  the  first  to  apply  the  maximum  principle  to  the  replacement  problem. 
Thomson  [1968],  Rapp  [1969],  Sethi  [1970]  and  many  others  used  these  so-called 
"optimal  control"  theory  models  that  uses  the  maximum  principle.  There  are  nu- 
merous variants  in  replacement  models,  the  assumptions  of  which  were  grouped  by 
Luxhoj  and  Jones  [1986]  in  a  general  framework.  Some  of  the  works  not  referenced 
in  this  survey  may  be  found  in  the  bibliography. 

Some  of  the  more  noteworthy  recent  publications  that  explicitly  considered 
technological  advances  are  as  follows: 

1.  R.C.  Stapleton  et  al.,  [1972]  examined  the  effects  of  different  forms  of  tech- 
nological change  (change  in  purchase  price,  resale  value,  operating  costs  and 
rate  of  deterioration)  on  costs  using  dynamic  programming.  This  work  pro- 
vided the  author  of  this  thesis  with  basic  ideas  on  how  to  model  the  impact 
of  technological  progress  on  cost. 

2.  S.  Chand  and  S.  Sethi  [1979  and  1982]  applied  the  planning  horizon  procedure 
developed  by  H.M.  Wagner  and  T.M.  Whitis  [1958]  for  the  dynamic  lot  size 
model  to  the  replacement  problem.  It  allows  one  to  obtain  the  first  replacement 
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time  for  an  infinite  planning  horizon  based  on  the  optimal  policy  for  a  finite 
horizon.  Technological  advances  is  forecasted  for  this  equivalent  finite  horizon. 
The  1982  paper,  which  extends  the  first  to  the  case  of  multiple  alternatives, 
treated  upgrades  within  the  same  framework  as  replacement  alternatives.  This 
idea  led  to  the  adoption  of  a  common  framework  for  all  alternatives  in  this 
thesis. 

3.  T.  Goldstein  et  al.,  [1988]  examined  the  influence  of  a  forthcoming  new  and 
better  machine  (based  on  an  expected  technological  breakthrough)  on  the 
system  renewal  decision.  Each  year,  if  the  new  machine  does  not  appear,  a 
decision  is  made  to  either  continue  with  the  incumbent  and  wait  for  the  bet- 
ter machine,  or  buy  a  new  version  of  the  incumbent.  The  model  determines 
the  optimal  age  for  replacing  the  existing  machine.  This  model  ignores  evolu- 
tionary changes  that  characterises  technological  innovation  and  thus  was  not 
adopted  in  this  thesis. 

4.  Y.  Kusaka  [1985]  showed  a  method  for  deciding  whether  to  keep  or  replace 
the  existing  machine  at  the  present  time  without  determining  the  sequence  of 
subsequent  replacement  times.  A  joint  paper  with  H.  Suzuki  [1988]  expanded 
on  the  earlier  evaluation  system  by  introducing  the  concept  of  "control  limit" 
policy.  Kusaka's  method  was  not  adopted  in  this  thesis  as  the  provision  of 
estimates  concerning  future  renewals  is  deemed  to  be  a  useful  feature. 
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III.  THE  MODEL 

A.  INTRODUCTION 

As  evident  in  the  previous  chapter,  many  complex  issues  impinge  upon  de- 
cisions concerning  systems  renewal  planning.  Literature  searches  have  failed  to 
uncover  any  model  suitable  for  system  renewal  planning  in  the  military.  The  dy- 
namic system  renewal  planning  (DSRPM)  model  to  be  presented  in  this  chapter 
is  a  first-cut  attempt  to  weave  a  unifying  framework  spanning  some  of  the  more 
quantitative  issues  involved  in  systems  renewal  planning.  The  overall  approach  is 
outlined  first,  followed  by  theoretical  underpinnings  of  the  model. 

B.  APPROACH 

The  adopted  approach  is  to  use  a  versatile  methodology,  namely  dynamic  pro- 
gramming augmented  by  the  Lagrangian  relaxation  technique,  to  explicitly  model 
the  effects  of  aging,  technological  advancement,  and  changes  in  capability  require- 
ments. The  methodology  operates  on  an  acyclic  directed  network  of  the  kind  shown 
in  Figure  3.1.  This  network  representation  of  the  system  renewal  process  is  unique 
to  this  thesis.  As  this  network  lies  at  the  heart  of  the  model,  a  brief  description  is 
in  order  before  we  move  on  to  the  other  parts  of  the  model. 
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Note  1  :  Nodes  denote  possible  transitions  while  arcs 
represent  possible  alternatives  and  system  life  spans. 
The  dotted  arc  from  node  3  to  7,  for  instance,  stands 
forthe  possibility  of  acquiring  alternative  2  in  1990  and 
using  it  as  System  B  until  1992. 
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Figure  3.1   :  Network  Representation  of  System  Renewal  Process 
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The  network  represents  the  system  renewal  process  involving  the  selection  of  a 
succession  of  systems  to  fulfill  some  capability  requirements.  The  root  node  1  repre- 
sents the  present  system  (System  A)  in  the  current  year.  Each  of  the  arcs  emanating 
from  it  represent  a  possible  residual  life-span  of  this  system.  These  arcs  terminate 
on  a  column  of  nodes  representing  possible  points  in  time  for  transition  to  the  next 
system  (System  B).  Transitions  may  be  replacements  or  upgrades.  Similarly,  arcs 
extending  from  these  nodes  to  an  adjacent  column  of  nodes  represent  the  possible 
years  for  transition  from  System  B  to  System  C.  The  number  of  arcs  connecting  each 
pair  of  nodes  corresponds  to  the  number  of  alternatives  that  may  be  selected  for 
System  B;  likewise  for  System  C.  Thus  if  cumulative  cost  and  effectiveness  of  using 
a  given  alternative  over  a  given  life-span  are  assigned  as  attributes  to  the  associated 
arc,  optimisation  for  the  three  systems  reduces  to  one  of  finding  the  "best"  path 
through  the  network. 

Finding  the  "best"  path  constitutes  solving  a  three  stage  dynamic  program- 
ming problem.  If  cost  is  set  to  a  very  large  number  whenever  an  alternative  ages 
such  that  one  or  more  performance  parameters  at  the  end  of  a  given  life  span  fails 
to  meet  stipulated  requirements,  then  dynamic  programming  yields  the  solution 
that  meets  the  requirements  at  minimum  cost.  However,  if  one  seeks  the  optimal 
trade-off  between  cost  and  effectiveness,  then  the  problem  becomes  a  multicriteria 
one.  This  problem  can  be  solved  using  the  Lagrangian  relaxation  technique  by  for- 
mulating this  problem  as  one  of  maximising  effectiveness  given  a  certain  budgetary 
limit  on  cost.  Details  of  this  technique  are  given  in  section  D  of  this  chapter.  The 
"best"  path  can  also  be  solved  using  linear  programming  formulation  of  the  network 
flow  problem.  However,  the  dynamic  programming  formulation  has  more  intuitive 
appeal  and  lends  itself  easily  to  incorporation  of  lead-time,  discounting  and  other 
considerations. 
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The  aging  process  is  modeled  by  discounting  performance  and  cost  character- 
istics over  time,  based  on  empirical  parametric  relationships.  What  these  charac- 
teristics will  be  for  a  given  alternative  acquired  new  in  a  certain  year  is  projected 
based  on  technological  trends.  Forecasts  of  inflation  and  foreign  exchange  rates  are 
also  factored  into  the  projected  acquisition,  maintenance  and  other  costs.  All  out- 
lays are  converted  to  current  year  dollars  based  on  forecasted  interest  rates.  The 
projected  rates  may  vary  from  year  to  year  but  a  steady  trend  is  assumed  beyond  an 
arbitrary  forecast  horizon.  Performance  requirement  levels  may  also  be  stipulated 
for  each  year  based  on  analysis  of  politico-military  trends.  A  unique  set  of  life-cycle 
cost  and  effectiveness  measures  is  thus  associated  with  the  network  for  each  scenario 
expressing  a  given  set  of  expectations  concerning  relevant  future  trends.  The  mini- 
mum cost  solution  and  efficiency  boundary  for  each  scenario  can  then  be  obtained 
by  applying  dynamic  programming  and  Lagrangian  relaxation  to  the  corresponding 
network.  The  efficiency  boundary  expresses  the  upper  bound  on  effectiveness  over 
the  entire  spectrum  of  life-cycle  cost  spanned  in  a  given  scenario  (Figure  3.2).  The 
point  on  this  boundary  that  yields  the  highest  effectiveness- to-cost  ratio  is  the  so- 
lution that  gives  the  most  bang  for  the  buck.  As  the  boundary  is  piecewise  linear, 
all  points  along  a  line  segment  may  qualify  if  extrapolation  brings  the  line  through 
the  origin.  Otherwise,  a  Lagrangian  solution  point  at  the  intersection  of  two  line 
segments  will  qualify. 

The  reader  may  wonder  why  a  three  stage  process  is  adopted  in  this  model. 
The  model  itself  does  not  restrict  the  user  to  any  number.  Rather,  the  matter  is 
one  of  judgment  concerning  the  credibility  of  long-term  forecasts.  Although  it  is 
desirable  to  optimise  over  as  long  a  stream  of  systems  as  possible,  it  is  also  futile 
to  try  to  visualise  the  distant  future.  A  two-stage  model  comprising  the  present 
system  and  its  immediate  successor  is  satisfactory  as  far  as  a  reasonable  optimisation 
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horizon  goes.  A  three-stage  model  would  be  preferable  as  it  allows  realistic  bounds 
to  be  developed  for  the  immediate  successor's  economic-technological  life.  The  cost 
of  competing  alternatives  can  therefore  be  more  equitably  compared  since  their 
amortisation  periods  may  be  more  accurately  estimated.  Infinite  horizon  models 
meant  for  static  environment  problems  are  clearly  not  suitable  for  our  purposes. 
Accordingly,  assumptions  of  monotonicity,  contraction,  boundedness,  etc.,  which 
are  required  for  asymptotic  stability  of  such  models  [Henig  1978,  Beckmann  1968, 
Boudarel  1971]  need  not  be  made. 

It  may  be  noted  that  a  long  lead  time  often  precedes  transitions,  be  they 
upgradings  or  replacements.  Lead  time  of  an  alternative  is  defined  to  mean  the 
minimum  number  of  years  by  which  the  initiation  of  acquisition  must  precede  the 
achievement  of  initial  operating  capability  (IOC).  By  setting  to  infinity  the  cost  of 
alternatives  that  cannot  be  introduced  because  of  insufficient  preceding  lead  time, 
we  can  avoid  generating  solutions  contradicting  lead  time  requirements.  This  check 
is  made  not  only  for  transitions  from  the  present  system  to  its  immediate  successor, 
but  also  for  subsequent  ones.  The  latter  requires  the  assumption  that  initiation 
of  a  transition  never  begins  before  the  preceding  one  even  takes  place.  A  path  in 
the  network  is  thus  deemed  infeasible  if  there  is  any  arc  representing  a  system  life 
shorter  than  the  lead  time  required  by  its  successor. 

While  no  one  can  reliably  predict  the  future,  one  ignores  future  possibilities 
and  trends  only  to  one's  disadvantage.  Taking  calculated  risks  after  examining 
a  broad  range  of  scenarios  is  the  only  realistic  approach  to  handling  uncertainties 
about  the  future.  For  this  purpose,  the  flexibility  of  the  model  presented  is  a  valuable 
asset.  Long-range  projections  can  be  varied  more  than  near-term  ones  to  account 
for  greater  uncertainty  in  these  projections. 
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It  should  be  emphasized  that  the  focus  of  analysis  is  on  the  timing  and  choice 
of  alternatives  for  the  first  transition,  the  optimality  of  which  is  affected  less  by 
inaccuracies  in  long-range  forecasts  than  short-term  ones.  By  updating  the  scenarios 
annually  for  changes  in  the  perception  of  qualified  experts,  rolling  plans  can  be 
developed  for  system  renewal  decisions  to  be  made  optimally  and  in  a  timely  fashion. 
For  this  reason,  the  model  is  named  the  Dynamic  System  Renewal  Planning  Model 
(DSRPM). 
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Figure  3.2  :  Domination  of  Solution  Space  by  Efficiency  Boundary 
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C.     ASSUMPTIONS 

For  the  purposes  of  this  thesis,  the  following  assumptions  have  been  used: 

•  A  steady  state  exists  beyond  an  arbitrary  forecast  horizon  where  all  projections 
are  either  constant,  linear  or  harmonic.  This  assumption  expresses  the  lack  of 
detailed  knowledge  about  the  remote  future. 

•  Cost  and  performance  parameters  can  be  measured  deterministically  and  util- 
ity functions  can  be  obtained  explicitly.  (Note:  This  assumption  would  be 
relaxed  in  future  embellishments  incorporating  stochastic  models  and  implicit 
approaches  to  utility  modeling.) 

•  An  upgraded  system  is  a  new  system  with  acquisition  cost  equal  to  the  direct 
upgrading  costs  plus  the  salvage  or  resale  value  of  the  original  system  in  the 
previous  year.  A  new  acquisition  with  no  precedent  is  simply  one  in  which 
the  incumbent  system  is  a  null  or  dummy  system.  With  this  assumption, 
upgrading,  replacement  and  new  acquisition  alternatives  are  integrated  within 
a  common  framework. 

•  Acquisitions  (upgrading  or  replacement)  and  their  associated  outlays  occur 
at  the  beginning  of  a  transition  year.  The  old  system  is  disposed  of  at  the 
end  of  the  year.  This  definitive  assumption  provides  for  proper  costing  of 
the  metamorphic  process  by  which  a  new  system  achieves  initial  operating 
capability  (IOC)  while  the  old  is  prepared  for  disposal.  Disposal  does  not 
really  occur  in  an  upgrading,  but  the  new  capability  is  still  assumed  to  be 
achieved  at  the  end  of  the  transition  year. 

•  No  transition  is  initiated  before  the  preceding  one  takes  place,  as  was  explained 
in  the  previous  section. 
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One  should  note  the  small  number  and  unrestrictive  nature  of  the  assumptions 
made  in  the  model.  This  is  a  direct  consequence  of  the  versatility  of  the  methodology 
employed.  There  are  other  assumptions  needed  to  model  particular  alternatives  and 
scenarios  such  as  those  pertaining  to  utility  functions,  nature  and  impact  of  trends, 
etc.,  but  these  are  not  generic  to  the  model  discussed. 

D.     METHODOLOGY 

Dynamic  programming  and  Lagrangian  relaxation  are  the  cornerstones  of  the 
model's  integrity  and  the  source  of  much  of  its  power.  This  section  explains  these 
key  techniques  as  applied  to  finite  discrete  decision  processes.  As  methodologies 
concerning  forecasting,  cost  estimation  and  other  components  of  the  model  may 
vary  according  to  the  peculiarities  of  the  application,  these  will  not  be  discussed. 
The  reader  is  referred  to  the  bibliography  for  materials  on  these  aspects  of  the  model. 

1.      Dynamic  Programming 

Borrowing  the  notation  used  by  the  pioneering  developer  of  dynamic  pro- 
gramming, Richard  E.  Bellman,  [Bellman  and  Kalaba,  1965],  we  begin  by  examining 
the  important  notion  of  policy  in  the  context  of  system  renewal.  If  a  decision  o?,-  in 
the  r'"*  stage  results  in  a  transition  from  a  state  y,-  to  a  state  t/,+i,  then  an  n-stage 
process  may  be  represented  as 

y\     =    T{yo,do) 
2/2     =     T{yudi) 

Vn     =    T{yn-udn-i)  (3.1) 

In  the  3-stage  system  renewal  process  depicted  in  Figure  3.1,  for  instance, 
the  process  of  operating  the  present  system  for  a  year,  then  acquiring  alternative 
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1  to  use  for  two  years  before  acquiring  anew  another  alternative  1  for  a  further 
two  years  use  would  be  represented  by  a  path  Pi,3j,i2  of  unbroken  line  arcs  passing 
through  nodes  1,  3,  7  and  12. 

Suppose  we  are  concerned  with  maximizing  a  prescribed  scalar  function 
of  the  states  and  the  decisions: 

R{yo,y\,y2,--- ,yn\  c?o,</i,---,<^n-i)  (3.2) 

This  so-called  criterion  function,  in  the  case  of  a  system  renewal  process,  may  be 
the  cumulative  expected  return  accrued  from  acquiring  and  operating  a  particular 
succession  of  alternatives,  i.e.,  traversing  a  given  path  in  the  network.  Without  an 
optimal  rule  as  to  how  decisions  should  be  made,  the  only  way  to  obtain  the  max- 
imum of  R  and  its  associated  path  is  by  exhaustive  enumeration  along  all  possible 
paths.  Effort  is  much  reduced  with  the  use  of  rules  in  the  form 

dk  =  dk{yo,yu- ■•  ^yk',  do,di,-- •  ,dk-i)  (3.3) 

These  functions  are  called  policy  functions,  or  simply  policies.  A  pohcy  that  maxi- 
mizes the  criterion  function  R  is  called  an  optimal  pohcy. 

Let  us  consider,  in  particular,  policies  which  depend  only  on  the  current 
state,  k^  that  is  policies  of  the  form 

dk  =  dk{yk)  (3-4) 

This  additional  simplification  is  possible  only  when  the  structure  of  the  criterion 
function  R  has  the  important  property  of  divorcing  the  past  from  the  future.  Ex- 
amples of  criterion  functions  with  this  "separable  criterion"  property  include 
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k 
and        Ylr{yi,di)  (3.6) 

The  first  example  is  of  the  most  common  type  of  criterion.  The  usual  criterion  of 
cumulative  cost  or  cumulative  returns  used  in  system  renewal  is  of  this  variety.  The 
reliability  of  series  systems  would  perhaps  come  to  mind  in  the  case  of  the  second 
one.  An  example  of  a  criterion  function  that  does  not  have  the  separability  property 
is  that  of  an  overall  ratio,  that  is, 

'^ (3.7) 


EC(^n^.) 

i=0 

The  use  of  overall  cost- effectiveness  ratio  as  a  criterion  in  the  renewal  process  would 
thus  require  decisions  to  be  based  not  only  on  the  current  state  but  also  on  the  path 
leading  to  that  state. 

When  the  criterion  function  used  has  this  vital  separability  property,  the 
optimal  pohcy  is  characterized  very  simply  by  the  following  principle  of  optimality: 

(yo,?/i,- ••  ,?/m,ym+i,-- •  ,^n)   is  Optimal  =>  (^q,  i/i  T' "^  ^m )   is  optimal         (3.8) 

This  intuitive  result  can  be  proved  by  contradiction.  Suppose  that  for  the  (m  -|-  1) 
stage  process,  (j/o,  J/i,  ■  •  ■ ,  J/m)  is  not  optimal.  Specifically  say  [xq^Xi,-  •  •  ,Xm)  is 
optimal.  Then  (xo,  .ti,  •••,  Xm^Vm+i,  •••,  yn)  is  better  than  (yo,yi,  ••• ,  t/m,  Vm+i, 
■  ■  ■  ,yn)-  But  this  contradicts  the  original  premise  that  (yo^yi?  ■••  ?  Vmi  Vm+i,-  ■  •  iVn) 
is  optimal.  Therefore,  expression  3.8  must  be  true.  Using  the  principle,  one  obtains 
for  the  examples  given  by  expressions  3.5  and  3.6  the  following  optimal  policies: 

fk+\{T{yk,dk))     =     max[c{yk,dk)  -\-  fkiVk)]  (3.9) 
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fk+i{T{yk,dk))    =    max{c{yk,dk)fk{yk)]  (3.10) 

where  Ji{yi,di)  is  the  criterion  function  in  state  i  obtained  using  an  optimal  policy. 

Dynamic  programming  hinges  around  the  formulation  of  optimal  policies 
to  enable  solutions  to  be  found  by  systematic  recursion  [Beckmann,  1968].  A  solution 
may  be  given  in  terms  of  the  sequence  of  criterion  functions  {/t(t/,)}  or  the  sequence 
of  policy  functions  {dk{yk)]  •  Each  determines  the  other,  although  one  should  note 
that  there  may  be  several  optimal  policies  yielding  the  same  sequence  of  criterion 
functions.  A  policy  may  also  be  satisfied  by  one  or  more  solutions.  Only  when  the 
solution  is  unique  is  the  policy  a  necessary  and  sufficient  condition  for  optimality. 
Otherwise,  it  is  necessary  to  compare  the  solutions  using  additional  considerations 
such  a^  risk,  expected  life-span,  etc.,  to  define  the  optimal  solution. 

Being  a  flexible  and  potent  tool,  dynamic  programming  has  many  vari- 
ations. In  stochastic  dynamic  programming,  the  next  state  is  chosen  ba^ed  on 
conditional  transition  probabilities  and  the  criterion  function  is  given  as  a  random 
variable,  usually  the  mean  of  possible  rewards  [Ross,  1983].  Variations  with  vector- 
valued  criterion  functions  have  been  developed  [Henig,  1978],  but  these  multicriteria 
dynamic  programming  models  are  still  not  as  widely  applied  as  the  stochastic  vari- 
ants. Although  these  multicriteria  models  could  be  used  in  system  renewal  problems 
to  obtain  dominating  (Pareto-efficient)  solution  sets,  it  is  simpler,  more  useful,  nat- 
ural and  elegant  to  use  the  Lagrangian  relaxation  technique  in  such  problems,  as 
shall  be  shown  in  the  following  section. 

2.      Lagrangian  Relaxation  Technique 

Two  conceptual  approaches  used  in  systems  analysis  are  either  to  fix 
effectiveness  and  minimize  cost  or  fix  the  budget  and  maximize  effectiveness  [Fisher, 
1970].  Suppose  we  choose  the  latter  and  formulate  the  system  renewal  problem  as 
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max  YlE{yi,di)  (3.11) 

subject  to 

i:C(y.,rf.)<  Budget  (3.12) 

1=0 

where  E{yi,di)  and  C(r/,, c/,)  are  the  life-cycle  effectiveness  and  cost  respectively 
of  making  a  decision  <i,  in  state  t/,.  Units  for  the  former  are  utiles,  while  that  of 
costs  are  dollars.  In  the  system  renewal  process  shown  in  Figure  3.1,  for  example, 
if  the  current  state  is  represented  by  node  3,  and  the  life- cycle  effectiveness  and 
cost  of  acquiring  alternative  1  and  operating  it  for  two  years  are  2,100  utiles  and 
$7m  respectively,  then  £'(3,2)  =  2,100  utiles  and  C(3,2)  =  $7m.  The  constraint 
in  expression  3.12  may  be  relaxed  by  moving  it  to  the  objective  function  with  the 
inclusion  of  a  Lagrangian  multiplier  A,  thus 


max   \j2^{y^,di)  -A 


^C(y.,<)- Budget 

t  =  0 


(3.13) 


Noting  that  expressions  3.5  and  3.13  are  of  the  same  form,  we  know  that  the  optimal 
policy  for  maximizing  this  criterion  function  is  given  by  the  principle  of  optimaUty 


as: 


fk^i{r{yk,d,))  =  max  {[E{yk,dk)  -  AC(y,,4)]  +  fk{yk)]  (3.14) 

Therefore,  for  any  fixed  A,  conventional  dynamic  programming  may  be  used  to 
obtain  the  optimal  solution  to  expression  3.13.  Lagrangian  relaxation  is  thus  a 
technique  for  converting  constrained  optimisation  problems  into  unconstrained  max- 
imization (or  minimization)  problems  [Everett,  1963].  It  may  also  be  interpreted  as 
a  technique  to  solve  the  dual  of  the  original  problem  given  by  Equations  3.11  and 
3.12,  as  will  be  apparent  in  the  following  discussion. 

To  consider  the  effect  of  varying  A,  let  us  define  the  Lagrangian  function 
Z/(A)  to  be  equal  lo  the  objective  function  given  in  expression  3.13.  That  is, 
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L(A)  =   max   Ij^Eivi^d,)  -  \ 

do,-,dn     [,^Q 


X]C'(?/,,c/,)  -  Budget 


Lt=0 


(3.15) 


This  piecewise  linear  function  is  illustrated  in  Figure  3.3. 
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Figure  3.3:   Lagrangian  Function  L{X) 


When  A  =  0,  the  budgetary  constraint  is  ignored  and  we  obtain  the 
maximum  effectiveness  solution  to  the  (71  +  1)  system  renewal  problem  given  in 
expression  3.13.  Denoting  the  cumulative  life-cycle  effectiveness  and  cost  of  this 
solution  as  LCEq  and  LCCo  respectively,  we  note  that  L{0)  =  LCEq  and  that  the 
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slope  at  this  vertical  axis  intercept  is  —{LCCo  —  Budget).  For  marginal  increases  in 
A  beyond  zero, 

L{\)  =  LCEo  -  \{LCCo  -  Budget)  (3.16) 

which  is  the  equation  of  the  leftmost  line  segment.  When  A  reaches  a  threshold 
value,  the  solution  changes  to  one  with  smaller  cumulative  life-cycle  effectiveness 
and  cost.  The  next  Hne  segment  will  thus  have  a  lower  intercept  value  and  an 
increased  slope.  This  therefore  implies  that  the  Lagrangian  function  will  always  be 
convex  for  system  renewal  problems  formulated  as  in  expression  3.13. 

The  Lagrangian  multiplier  A  has  a  physical  interpretation,  namely  that 
of  the  relative  value  of  cost  and  effectiveness,  i.e.,  shadow  or  dual  price.  Its  unit  is 
utiles  per  dollar.  As  A  is  increased  from  zero  to  infinity,  the  Lagrangian  function 
steps  through  the  entire  solution  space,  from  the  maximum  effectiveness  solution  to 
the  minimum  cost  solution.  The  latter  corresponds  to  the  rightmost  hne  segment 
which  extends  to  infinity.  What  then  is  the  optimal  value  of  A?  To  see  this,  we  first 
note  that  for  A  >  0,  L{\)  is  an  upper  bound  on  the  optimal  value  F*  of  the  original 
objective  function  given  in  expression  (3.11).  This  is  because  minL(A)  is  the  dual 
of  the  original  problem.  The  minima  of  L{\)  therefore  yields  F* .  That  is, 

minL(A)  =  L(A*)  =  F*  (3.17) 

This  is  akin  to  seeking  the  constrained  stationary  points  in  analog  applications  of 
Lagrange's  method  to  differentiable  functions. 

The  minima  of  L{\)  do  not  necessarily  exist.  In  a  commonly  observed 
case,  LCCo  is  higher  than  the  budget  and  the  slope  at  the  vertical  intercept  is 
negative.  As  L{\)  is  convex,  A*  is  always  finite  and  non-zero  in  such  cases.  If, 
however,  LCCo  is  lower  than  the  budget,  that  is  to  say  the  maximum  effectiveness 
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solution  is  affordable,  then  the  slope  at  the  intercept  will  be  positive  and  we  obtain 
A*  =  0.  Should  LCCq  be  higher  than  the  budget  and  there  are  no  affordable 
solutions,  then  there  will  be  no  line  segments  of  non-negative  slopes,  no  minima  and 
hence  no  feasible  solution. 

The  value  of  A*  can  be  found  most  expeditiously  and  accurately  by  the 
supporting  plane  method  illustrated  in  Figure  3.4.  Given  the  upper  and  lower 
bounds  on  A*,  the  method  involves  the  extrapolation  of  the  line  segments  corre- 
sponding to  these  bounds  to  obtain  the  intersection  value  A,.  The  intercept  and 
slope  of  these  line  segments  are  given  by  the  cumulative  life-cycle  effectiveness  and 
cost  respectively  of  the  solution  obtained  by  solving  expression  (3.13)  for  the  bound- 
ing values  of  A.  The  lower  or  upper  bound  is  replaced  by  A,  if  the  slope  of  the  line 
segment  corresponding  to  A,  if  the  slope  of  the  line  segment  corresponding  to  A,  is 
negative  or  non- negative  respectively.  This  process  is  iterated  until  A^  equals  one  of 
the  bounds.  This  equality  occurs  because  L{\)  is  piecewise  linear  and  convex,  which 
enables  A*  to  be  obtained  precisely.  The  starting  lower  bound  is  set  to  be  zero.  Any 
suitably  large  value  of  A  that  yields  a  line  segment  with  non-negative  slope  may  be 
used  as  the  initial  upper  bound  Xmax,  where 

Amax  =  No.    of    nodes    X    LCEmax  -^  LCCmin  (3.18) 

Experience  of  the  author  with  similar  problems  shows  this  supporting  plane  method 
to  be  faster  and  more  accurate  than  the  commonly  used  bisection  method,  and  A* 
may  be  found  in  only  a  few  iterations  even  for  very  large  networks  with  thousands 
of  nodes  and  arcs. 

For  each  scenario,  using  the  methodology  described,  we  are  able  to  obtain 
the  solution  points  marking  the  efhciency  boundary  shown  in  Figure  3.2.  This 
is  begun  by  artificially  setting  a  sufficiently  high  budget  to  obtain  the  maximum 
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Figure  3.4:  Iterations  for  Derivation  of  A* 

efFectiveness  solution,  which  is  marked  1  in  Figure  3.5.  The  corresponding  value 
of  A*  is  obviously  zero.  The  budget  is  then  set  to  just  under  the  value  of  the 
cumulative  life-cycle  cost  C\  of  the  maximum  effectiveness  solution.  This  produces 
the  next  solution  point,  marked  2  in  the  diagram,  and  the  process  is  repeated  until 
we  reach  the  minimum  cost  solution,  which  is  the  sixth  solution.  No  solutions  will  be 
obtained  for  budgets  less  than  Cg.  The  optimal  Lagrangian  multiplier  A*  obtained  in 
each  case  determines  the  slope  of  the  line  segment  to  the  right  of  the  corresponding 
solution  point.  It  may  be  noted  that  the  efficiency  boundary  consists  of  bounding 
hyperplanes  which  envelope  the  solution  space  but  do  not  pass  through  all  the 
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Pareto-efficient  points.  There  are,  for  instance,  three  such  points  denoted  by  circles 
for  budgets  between  C3  and  C4  which  will  not  be  discovered  by  the  Lagrangian 
relaxation  method.  The  basic  cause  of  these  inaccessible  regions  (called  gaps)  is 
non-concavity  in  the  relationship  between  the  optimal  effectiveness  and  cost.  These 
gaps  will  not  exist  if  the  law  of  diminishing  marginal  utility  holds  strictly.  H. 
Nakayama,  et  al.,  [Nakayama,  1975],  has  a  good  survey  of  methods  to  deal  with  this 
characteristic  of  the  Lagrangian  relaxation  method.  The  existence  of  these  gaps  does 
not  create  a  problem  as  far  as  the  application  of  the  method  in  systems  analysis  is 
concerned.  This  is  because  the  solution  offering  a  maximum  overall  effectivenss- to- 
cost  ratio  is  always  a  Lagrangian  solution  point.  Also,  a  straightforward  application 
of  dynamic  programming  will  produce  the  minimum  cost  solution  even  if  it  does  not 
lie  on  the  efficiency  boundary.  Note  that  as  system  lives  are  modeled  explicitly  for 
economic,  technological  and  politico-military  changes,  the  net  cost  of  alternatives 
can  be  amortized  over  the  proper  period  to  reveal  the  true  minimum  cost  solution 
as  defined  by  amortized  annual  cost,  as  will  be  amplified  in  the  following  paragraph. 
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We  have  seen  how  characterization  of  the  systems  renewal  problem  as  a 
constrained  optimization  problem  (in  which  one  seeks  to  maximize  effectiveness  for  a 
given  budget)  enables  us  to  obtain  the  efficiency  boundary  corresponding  to  a  given 
scenario.  A  band  of  these  boundaries  is  found  for  each  set  of  scenarios  investigated. 
If  probabilities  are  attached  to  scenarios,  our  "expected"  efficiency  boundary  and 
its  associated  optimal  solutions  can  be  obtained.  Furthermore,  bearing  in  mind  that 
the  cost  of  an  arc  in  the  network  is  set  to  infinity  if  there  is  insufficient  lead  time 
preceding  it  or  if  any  performance  parameter  deteriorates  below  its  requirement  level, 
we  know  that  any  solution  with  finite  cost  is  feasible  from  the  effectiveness  point  of 
view.  The  minimum  cost,  in  current  year  dollars,  for  the  various  feasible  solution 
paths,  as  well  as  the  associated  period  of  use,  can  then  be  translated  into  amortized 
costs  flows  for  given  interest  rate  assumptions.  On  this  basis,  we  can  find  the  solution 
that  allows  us  to  minimize  equivalent  annual  cost  while  meeting  a  given  set  of  lead 
time  and  performance  requirements.  Therefore,  both  the  fixed  budget  and  the  fixed 
effectiveness  approaches  are  embodied  within  the  model,  delivering  the  maximum 
overall  effectiveness-to-cost  ratio  solution  and  the  minimum  cost  satisficing  solution 
respectively.  Furthermore,  it  should  be  noted  that  A*  has  the  interpretation  of 
shadow  or  dual  prices  in  economic  analysis.  Therefore,  efficiency  boundaries  of 
systems  with  similar  MOEs  can  be  meaningfully  compared.  Also,  choice  of  system 
mix  for  various  capabilities  may  be  guided  by  the  fact  that  shadow  prices  should 
be  roughly  the  same  for  all  the  systems  in  an  optimal  mix.  The  combined  use 
of  dynamic  programming  and  the  Lagrangian  relaxation  technique  thus  provide  us 
with  a  conceptually  simple,  yet  potent  and  useful  methodology  that  is  completely 
congruent  with  and  derivable  from  the  tenets  of  systems/economic  analysis. 
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IV.  IMPLEMENTATION 

A.     INTRODUCTION 

The  dynamic  systems  renewal  planning  model  (DSRPM)  presented  in  the  pre- 
vious chapter  provides  a  formal  basis  by  which  diverse  economic,  technological  and 
politico-military  concerns  may  be  integrated  within  a  common  framework  (Figure 
4.1)  to  make  for  more  rational  management  of  defense  assets.  Credible  and  valid 
results,  however,  can  be  obtained  only  with  the  active  participation  of  relevant  logis- 
tics, engineering,  finance,  intelligence  and  operational  agencies  in  the  development 
of  proper  effectiveness  and  cost  sub-models  as  well  as  realistic  scenarios.  High-level 
management  support  is  thus  essential  for  the  successful  implementation  of  the  model. 
To  foster  a  better  understanding  of  the  model  and  to  demonstrate  its  plausibility 
and  usefulness,  a  small  prototype  was  developed  based  on  a  simple  hypothetical 
weapon  system.  This  prototype  is  documented  in  the  next  chapter  along  with  case- 
studies.  Although  the  addition  of  selected  successful  applications  would  give  greater 
credence  to  the  model,  the  limited  time  available  for  work  on  this  thesis  precluded 
full-scale  implementation  for  any  real  systems.  Nevertheless,  effort  shall  be  made 
to  outline  conceptually  how  implementation  is  envisaged.  Only  aspects  of  cost  es- 
timation, effectiveness  modelling  and  requirement  formulation  peculiar  to  DSRPM 
is  discussed.  For  a  full  discussion,  the  reader  is  referred  to  relevant  materials  in  the 
bibliography. 
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B.     CONCEPTUAL  FRAMEWORK 

Implementation  of  the  model  is  best  explained  in  terms  of  its  constituent  el- 
ements and  their  relationships  to  the  processes  involved  in  defense  planning.  As 
depicted  in  Figure  4.1,  DSRPM  ha^  five  elements,  all  of  which  are  part  of  the 
force/resource  planning  process.  A  model  with  these  five  elements  is  developed 
for  each  entity  of  force  planning,  which  may  be  weapon  systems,  task  forces,  bat- 
tle groups,  etc.  In  renewal  planning  parlance,  all  these  entities  are  simply  called 
systems,  and  candidate  systems  are  called  alternatives. 

The  cost  estimation  element  projects  the  life-cycle  costs  for  both  the  present 
system  as  well  as  future  alternatives.  The  effectiveness  modelling  element  concerns 
the  appraisal  of  the  value  of  a  system  as  a  function  of  its  measures  of  effectiveness. 
The  deterioration  modelling  element  accounts  for  the  change  in  effectiveness  and 
cost  over  the  life  of  a  system  caused  by  wear-out,  operational  attrition  and  enemy 
response.  The  other  two  elements  are  self-explanatory.  All  of  these  elements  and 
their  implementation  are  elaborated  in  sections  C,  D,  and  E.  Only  deterioration 
modelling  and  the  selection  and  scheduling  of  programs  are  organic  to  the  DSRPM. 
The  other  elements  are  on-going  force/ resource  planning  processes  adapted  with 
future  projections  for  use  in  the  DSRPM. 

Formulation  of  future  alternatives  is  a  part  of  the  force/resource  planning  pro- 
cess that  lies  wholly  outside  the  province  of  the  DSRPM.  It  involves  definition  of 
the  physical  and  functional  forms  of  alternatives  identified  in  the  strategy  formu- 
lation process  and  is  shaped  by  technological  outlook.  Alternatives  encompass  the 
following: 

1.  Conversion  of  incumbent  systems  to  serve  different  roles 

2.  Upgrading  of  incumbent  systems  to  serve  expanded  or  higher  roles 
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3.  Modernisation  of  incumbent  systems  to  better  serve  their  assigned  roles 

4.  Replacement  of  incumbent  systems  with  newer  and  probably  better  alterna- 
tives 

5.  Activation  of  reserve  systems  for  active  deployment  (incumbent  systems  are 
the  ones  in  dormant  or  preserved  status) 

6.  Reduction  of  number  of  incumbent  systems  through  disposal/sale 

7.  Retirement  of  incumbent  systems  into  the  reserves 

8.  Acquisition  of  new  (no  precedent)  systems  to  meet  qualitative  and/or  quanti- 
tative deficiencies  or  new  (threat)  requirements. 

No  present  system  exists,  i.e.,  system  A  in  the  network  shown  in  Figure  3.1 
does  not  exist,  in  the  case  of  new  acquisitions. 

C.      EFFECTIVENESS  MODELING 

The  effectiveness  modeling  adopted  in  the  DSRPM  is  completely  congruent 
with  the  way  military  systems  are  selected  and  managed.  Each  system  has  one  or 
more  missions  that  represent  the  raison  d'etre  of  that  system.  These  missions  are 
usually  phrased  in  terms  that  are  too  broad  for  concrete  guidance.  What  it  takes  to 
carry  out  the  mission(s)  are  stipulated  clearly  as  required  operational  (capabilities) 
requirements.  These  are  usually  lower  bounds  set  on  key  performance  parameters 
such  as  speed,  payload,  survivability,  operational  range,  etc.  The  effectiveness  of 
candidate  systems  is  evaluated  based  on  these  key  parameters  and  some  criterion. 
The  decision  to  choose  a  particular  system  need  not  be  based  exclusively  on  objec- 
tive parameters,  and  the  criterion  used  is  not  necessarily  a  formal  one.   The  views 
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of  qualified  experts  and  decision  makers  are  also  important.  In  fact,  the  politi- 
cal process  leading  to  the  decision  may  even  be  colored  by  elements  that  detract 
from  the  mihtary  contribution  to  national  security,  such  as  amicable  relations  with 
another  state.  However,  assuming  rational  arguments  hold  sway,  an  effectiveness 
model  representing  the  embodiment  of  the  values,  preferences  and  possibly  biases 
of  the  decision  makers,  tempered  with  experience  gained  from  performance  audits, 
can  be  constructed.  Justification  papers,  test  and  evaluation  plans,  trial  reports, 
etc.  are  good  sources  of  information  needed  to  construct  an  effectiveness  model. 
Combat  models  and  wargames  also  provide  valuable  input.  The  aim  of  effectiveness 
modelling  for  this  model  is  to  describe  value  judgements  of  the  decision-makers  and 
not  to  predict  battle  outcomes  (although  predictions  do  tend  to  affect  judgement). 
Therefore,  the  same  effectiveness  model  used  for  evaluation  of  alternatives  in  the 
acquisition  phase,  and  possibly  altered  by  performance  audits  during  deployment, 
is  incorporated  in  the  DSRPM.  The  only  difference  in  the  DSRPM  case  is  that 
prospective  alternatives  are  evaluated  alongside  the  incumbent,  and  performance 
parameters  are  projected  both  for  future  generations  of  each  alternative  as  well  as 
through  the  years  of  use  after  an  alternative  is  acquired.  Figure  4.2  shows  how  the 
j"*  performance  parameter  of  an  alternative  as  acquired  in  year  y  may  change  over 
the  years  of  usage  t. 

As  an  example,  let  us  suppose  that  a  weapon  system  was  evaluated  based  on 
an  effectiveness  model  with  the  following  additive  structure: 


TE,    =    EHV7,{X;,a)    provided  X,-„  >  X,  for  V; 

=     0  otherwise  (4.1) 
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where 

•  TEa  is  the  total  effectiveness  of  alternative  a. 

•  Wj  is  the  weight  of  the  j      performance  parameter  expressing  its  importance 
relative  to  the  other  parameters. 

•  Xj^a  is  the  value  of  the  j      parameter  for  alternative  a. 

•  Uj  {Xj^a}  is  the  utihty  function  that  transforms  Xj^a  into  a  utility  value. 

•  2Lj  is  the  minimum  level  required  of  the  j"'  parameter  (operational  require- 
ment). 

In  this  case,  the  effectiveness  sub-model  used  in  the  application  of  DSRPM  to  the 

system  will  be: 

m 

TEa'M  =  E  "'.^i  {^,,<yda'At)]  (4.2) 

i=i 

where 

•  TEa',y{t)  is  the  total  effectiveness  of  alternative  a'  (not  necessarily  the  same  as 
a)  if  acquired  new  in  year  y  during  the  t^^  year  of  deployment. 

•  Xjy^y  is  the  value  of  the  j"*  performance  parameter  for  alternative  a'  if  ac- 
quired in  year  y. 

•  daj{t)  is  the  deterioration  function  for  the  j"*  parameter  of  alternative  a'  in 
the  t^^  year  of  deployment. 

The  alternatives  a  and  a'  need  not  be  the  same.  The  latter  typically  comprises 
not  only  the  replacement  alternative,  but  also  upgrading  options.  Note  that  the 
structure  of  the  effectiveness  model,  weights  and  utility  functions  are  retained  in 
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the  DSRPM.  TEa>,y{t)  is  the  effectiveness  associated  with  all  arcs  representing  al- 
ternative a',  leading  out  of  nodes  of  year  y,  with  life-spans  of  t  years  in  the  network 
shown  in  Figure  3.1.  If  Xa',j{t)  <  2Lj,y+t  ^^r  any  j,  i.e.,  if  any  parameter  fails  to 
meet  the  requirement  level  stipulated  for  the  year  of  operation  {y  -f  t),  then  the  cost 
of  the  corresponding  arc  is  set  to  infinity. 

The  effectiveness  model  used  here  is  similar  to  that  of  the  TASCOM  method- 
ology developed  by  The  Analytic  Sciences  Corporation  for  OSD  over  several  years 
in  the  late  1970's  and  the  early  1980's.  The  TASCOM  model  provides  a  first-order 
indication  of  an  aircraft  performance  potential  by  computing  a  score  based  on  at- 
tributes such  as  range,  speed,  maneuverability,  payload,  basing  mode,  target  acqui- 
sition and  fire  control  capabilities.  These  factors  are  normalised  against  a  consistent 
standard  or  assigned  a  graduated  evaluation  score.  The  relative  weights  are  based 
on  the  consensus  judgment  of  experienced  pilots  as  well  as  military  analysts.  The 
final  performance  are  scores  based  on  weighted  calculations.  The  reader  is  referred 
to  Hildebrandt,  et  al.,  1986,  and  Congressional  Budget  Office,  ^^  Tactical  Combat 
Forces  of  the  USAF,  Issues  and  Alternatives,^^  April  1985,  for  further  details  on  the 
TASCOM  methodology. 

The  values  of  the  parameters  for  future  versions  of  each  alternative  are  pro- 
jected as  Xj^a',y  These  are  discounted  (for  deterioration)  through  the  years  of  usage 
by  the  function  dj{t).  Technical  expertise  is  required  for  such  forecasting.  Require- 
ment levels  2Lj,y+t  need  not  be  static  but  may  be  varied  through  the  future  years 
in  accordance  with  projected  threats.  Figure  4.2  illustrates  how  this  requirement 
level  may  change  with  time.  Projections  are  required  only  for  the  span  of  the  un- 
derlying network  (Figure  3.1)  and  a  steady  state  should  be  assumed  beyond  the 
forecast  horizon.  These  projections  form  the  effectiveness  aspect  of  DSRPM  scenar- 
ios. Examples  of  Uj  {    ]  and  (la',j{    )  arc  shown  in  Figures  4.3  and  4.4  respectively. 
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The  utility  function  Uj  {  }  is  the  instrument  that  allows  us  to  express  non-linear 
relationships  such  as  that  of  diminishing  marginal  utility.  Adjusting  it  regularly  for 
the  effect  of  changes  in  enemy  force  structure,  size,  tactics,  strategy,  etc.,  also  serves 
to  account  for  the  interactive  nature  of  defense  "posturing".  The  example  shown 
here  is  of  explicit  utihty  functions  and  single  mission  systems.  Implicit  approaches 
may  be  used  if  utility  functions  can  be  specified  explicitly  only  with  considerable 
difficulty.  Multiple  mission  systems  may  be  modelled  by  combining  utility  values 
using  weights  or  other  multi-criteria  modelling  schemes. 
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Figure  4.2  a  :  Change  in  MOE  Value 


Figure  4.2  b  :  Change  in  Requirement  Level 
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Figure  4.2  c  :  Change  in  MOE  Value  and  Requirement  Level 
Figure  4.2  :  Illustration  of  Change  in  MOE  Value  and  Requirement  Level 
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Figure  4.3  :  Examples  of  Utility  Functions 
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Figure  4.4:  Examples  of  Deterioration  Functions 

D.      COST  ESTIMATION 

Resource  requirements  are  evaluated  for  both  incumbent  systems  and  prospec- 
tive alternatives.  Cost  models  for  the  former  deal  with  the  projection  of  ownership 
costs  based  on  actual  cost  experience  as  well  as  forecasts  of  disposal  cost  or  salvage 
value.  Those  models  developed  for  future  alternatives  are  either  statistical  ones 
based  on  the  physical  and  functional  characteristics  of  these  alternatives  (Figure 
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4.1)  or  analogy  types  based  on  comparison  of  similar  systems  and  judgement  con- 
cerning complexity  factors.  Procurement  costs  are  included  in  cost  estimation  of 
alternatives  but  are  considered  sunk  costs  in  the  case  of  incumbent  systems. 

Cost  models  used  in  DSRPM  differ  from  the  ones  typically  used  in  evaluations 
in  that  various  types  of  costs  such  as  procurement  cost,  operating  cost,  etc.,  are 
estimated  not  just  for  one  possible  life-cycle  but  for  multiple  possible  ones  in  the 
future.  The  projection  of  costs  needed  to  adapt  cost  models  for  use  in  DSRPM  is 
usually  accomplished  using  time-series  or  regression  analysis.  Adjustments  may  be 
made  if  significant  and  relevant  technological  breakthroughs  are  anticipated.  For  ex- 
ample, if  major  advances  in  certain  automation  technologies  are  expected  five  years 
from  now,  then  manpower  costs  may  be  adjusted  accordingly.  Such  adjustments 
should  be  made  conservatively  and  in  a  not  too  abrupt  manner.  Figure  4.5  shows 
an  example  of  the  cost  parameter  forecasts  that  are  used  in  the  model. 

To  examine  how  costs  may  be  projected  into  the  future,  let  us  consider  the 
following  simplified  life-cycle  cost  model: 

LCC.,y[i)  =  ACa,y  +  MC„,,r-''  -f  OCa,ye*-^('-^)  +  DC^Ai)  (4-3) 

DC^M  =  U,,y  -  Va,ye-"''^'  (4.4) 

where: 

•  LCCa,y{t)  is  the  hfe-cycle  cost  of  alternative  a  if  procured  in  year  y  and  used 
for  t  years 

•  ACa,y  is  the  acquisition  cost  of  alternative  a  if  procured  in  year  y 

•  MCa,y  is  the  maintenance  cost  of  alternative  a  in  the  first  year  if  procured  in 
year  y 
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Figure  4.5:  Examples  of  Cost  Parameter  Forecast 

•  OCa,y  is  the  operating  cost  of  alternative  a  in  the  first  year  if  procured  in  year 

y 

•  DCa,y{t)  is  the  disposal  cost  of  alternative  a  if  procured  in  year  y  and  used  for 
t  years  (Note:  This  may  be  negative  if  asset  has  resale  value.) 

•  Ua,y  is  the  upper  bound  of  DCa,y{t) 

•  Va,y  is  a  proportionality  constant  dictated  by  the  resale  market 

•  7a,y5^a,y  ^jid  7]a^y  are  exponential  constants 
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If  inflation  and  technological  change  have  been  (and  presumably  will  be)  gradual, 
the  individual  costs  may  be  forecasted  using  non-linear  regression  models  such  as 
the  following: 

ACa,y  =  PAC,{1  +  tacY  +  Cy  (4.5) 

where: 

•  PACa  is  the  present  acquisition  cost  of  alternative  a 

•  Vac  is  the  compound  inflation-technological  improvement  rate  for  acquisition 
cost 


•  Cy  is  the  error  term 


Solutions  are  obtained  using  either  the  maximum  likelihood  approach  or  Zellner's 
method  [Judge,  et  al.,  1988]  based  on  historical  data.  If  major  breakthroughs  are 
expected  in  the  production  of  alternative  a,  the  forecasted  figures  for  ACa,y  may 
be  adjusted  to  some  extent.  Anticipated  foreign  exchange  rates  may  also  be  fac- 
tored into  cost  projections  if  they  are  relevant  to  a  system  renewal  problem.  The 
cost  aspect  of  DSRPM  is  manifested  in  these  projections.  There  are  many  differ- 
ent methods  of  forecasting  both  cost  and  efi^ectiveness  [Ayres,  1969,  Stewart  and 
Wyskida,  1987,  and  Ascher  and  Overbold,  1988],  the  choice  of  which  depends  very 
much  on  the  degree  of  understanding  one  has  of  the  dynamics  of  the  change  process. 

E.     OTHER  ELEMENTS 

The  deterioration  function  dj{t)  introduced  in  the  section  on  effectiveness  mod- 
elling is  an  example  of  how  tiic  change  in  performance  parameters  over  the  life  of  a 
system  may  be  modelled.  Audit  programs  such  as  0 T&E,  inspections  and  exercises 
as  well  as  technical  trials  and  laboratory  tests  provide  the  data  needed  to  derive 
these  functions.  To  see  how  deterioration  may  be  incorporated  in  a  cost  model,  let 
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us  consider  the  following  maintenance  cost  model: 

dM 
"df 

where: 


=  cArt^  (4.6) 


•  M  is  the  cumulative  maintenance  cost  incurred 

•  t  is  the  age  of  the  system 

•  c  is  a  proportionality  constant 

•  a  and  /3  are  exponentiality  constants 

This  model  links  the  rate  of  change  in  maintenance  cost  to  the  age  t  as  well  as  past 
maintenance  effort  measured  in  terms  of  cumulative  cost.  Given  an  initial  condition 
of  M  =  0  at  <  =  0,  we  have: 

M  =  AP  (4.7) 

where: 

7  ^  P^  (4.8) 


(1-a) 


A  = 


(4.9) 


Q  /  1  (4.10) 

Hence  the  rate  of  increase  in  maintenance  cost  is  given  by: 

^  =  yl7r-'  (4.11) 

The  parameters  A  and  7  are  obtained  by  using  regression  analysis  on  historical 
maintenance  cost  date.  A  is  the  maintenance  cost  incurred  in  the  first  year  of  use. 
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As  7  is  always  greater  than  1,  the  rate  of  increase  in  maintenance  cost  ^  will 
always  be  increasing  over  time.  The  impact  of  inflation  and  technological  advances 
in  maintenance  cost  is  modelled  by  varying  both  A  and  7  for  future  alternatives.  A 
reduction  in  operating  cost  for  future  generations  may  also  be  modelled  in  a  similar 
fashion. 

Operational  requirements  for  present  systems  are  projected  by  defining  the 
minimum  requirement  level  in  each  year  for  various  performance  parameters  of  a 
system.  These  are  denoted  as  2Lj,y-i-t '"  section  B.  For  instance,  shifts  in  enemy  beach 
defence  tactics  may  result  in  significant  increases  in  the  speed  required  of  amphibious 
assault  craft.  Requirement  projections  form  an  integral  part  of  a  DSRPM  scenario 
and  are  driven  mainly  by  the  net  assessment  process  (Figure  4.1). 

Force  levels  decline  even  in  peacetime  with  attrition  due  to  accidents,  wear 
and,  in  the  case  of  human  systems,  retirement  and  resignation.  If  the  number  of  an 
alternative  a'  acquired  in  year  y  is  taken  to  be  the  performance  parameter  Xj^a',y, 
the  decline  due  to  attrition  can  be  modelled  by  the  deterioration  function  daj{t). 
The  effectiveness  of  various  force  levels  can  then  be  modelled  using  the  associated 
utility  function  Uj  {},  and  the  force  level  requirement  in  the  year  {y  -\- 1)  is  simply 
2Lj,y+f  Treating  the  number  acquired  as  a  performance  parameter  in  this  way 
allows  us  to  check  the  adequacy  of  the  size  of  an  acquisition.  The  attractiveness  of 
a  phased  build-up  versus  the  economies  of  a  bulk  buy  may  be  studied  by  defining 
alternatives  representing  different  systems  and  force  levels.  DSRPM  produces  an 
acquisition  program  that  meets  projected  force  level  requirements  by  selecting  and 
scheduling  these  alternatives.  Figure  4.6  illustrates  the  kind  of  solution  that  may  be 
generated  by  DSRPM  in  such  studies.  Force  mix  can  be  optimised  by  ensuring  that 
the  lagrangian  multiplier  associated  with  the  levels  of  the  various  force  components 
are  similar.  In  this  way,  force  level  and  mix  issues  can  be  dealt  with  by  the  model. 
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then  one  of  12  for  the  build-up  followed 
by  one  of  6  to  moke  up  for  attrition  loss. 


10 
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15 


20 


Figure  4.6:   Phased  Acquisition  Solution  Example 

The  last  element  we  look  at  concerns  the  selection  and  scheduling  of  alterna- 
tives for  implementation.  The  inputs  for  these  are  the  cost  and  effectiveness  models 
discussed  in  sections  B  and  C  of  this  chapter  as  well  as  budgetary,  temporal,  legal 
and  other  constraints.  The  combined  dynamic  programming  -  Lagrangian  relax- 
ation discussed  in  Chapter  3  is  applied  to  these  inputs  to  produce  solutions  for  the 
fixed  budget  and  fixed  effectiveness  approaches.  Constraints  in  the  form  of  quanti- 
tative rules  can  be  incorporated  within  DSRPM.  An  example  is  that  of  lead  time 
requirements  for  various  alternatives  which  is  implemented  in  the  prototype  that 
was  developed.  The  model  is  exercised  by  iterating  it  for  various  scenarios  to  make 
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explicit  uncertainties  faced  by  the  decision  maker.  Strategic  considerations  with 
the  formal  quantitative  results  from  DSRPM  are  then  the  basis  of  system  renewal 
decisions  which  are  the  output  of  the  force/resource  planning  process. 

DSRPM  is  easily  implemented  in  phases.  The  model  could  first  be  applied  to 
a  couple  of  incumbent  systems  for  which  systems  renewal  decisions  are  expected  in 
the  next  one  or  two  years.  This  would  allow  enough  time  for  a  full  implementation 
of  the  model  for  these  systems.  This  phase  serves  to  validate  the  model  with  a 
real  application  and  the  experience  gained  would  guide  implementation  in  the  next 
phase.  The  model  should  however  be  applied  as  early  as  possible  in  a  system's 
life-cycle  to  provide  enough  time  for  the  cost  and  effectiveness  models  to  be  refined 
with  inputs  from  operational  trials,  exercises,  wargames,  etc.  Therefore,  the  second 
phase  should  see  the  model  applied  to  some  new  acquisitions.  When  management  is 
comfortable  with  the  model  and  is  ready  to  introduce  it  on  a  broader  scale,  a  master 
implementation  plan  can  then  be  formulated  to  formalise  the  use  of  the  model. 
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V.  PROTOTYPE  DESIGN 

A.      GENERAL 

To  foster  a  better  understanding  of  the  DSRPM,  a  prototype  was  developed 
based  on  a  hypothetical  weapon  system.  This  chapter  documents  the  prototype 
which  was  developed  on  an  IBM  AT-compatible  personal  computer  using  PC  SIM- 
SCRIPT  II. 5.  The  language  was  chosen  for  its  superlative  power  and  ease  of  use, 
as  well  as  its  provisions  for  incorporation  of  simulation  in  later  development.  The 
program  is  designed  to  produce  for  each  run  a  single  solution  path  corresponding  to 
a  point  on  Figure  3.2.  The  program  reads  data  from  an  input  file  "IN. DAT"  and 
channels  results  to  an  output  file  "OUT. DAT".  The  type  of  solution  depends  on  the 
approach  selected  for  that  run,  namely  Lagrangian  relaxation,  maximum  effective- 
ness, minimum  satisficing  cost  or  maximum  system  cost  effectiveness  ratio.  Solution 
approach  is  specified  by  the  user  in  the  input  file.  The  budget  is  relevant  only  when 
the  Lagrangian  method  is  used.  The  overall  structure  of  the  prototype  is  described 
in  Figure  5.1  and  the  program  is  listed  in  Appendix  A. 
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Figure  5.1  :  Prototype  DSRPM  Block  Diagram 
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B.     ASSUMPTIONS 

The  following  assumptions  are  peculiar  to  the  hypothetical  weapon  system 
upon  which  the  prototype  was  developed: 

•  The  efTectiveness  sub-model  is  given  by  Equation  4.3  which  was  discussed  as  an 
example  in  the  previous  chapter.  Deterioration  functions  daj{t)  are  assumed 
to  be  exponential  for  all  a  and  j.  Furthermore,  utility  functions  Uj  {X}  are 
assumed  to  be  linear  and  bounded  (maximum  value  of  1)  as  follows: 

max  |A  I 

Therefore,  the  total  effectiveness  of  an  alternative  a  acquired  in  year  y  in  the 
<"*  year  of  use  is  given  by: 


TE^M    =    j:^V,Uj{X,,^J^Jt)]  (5.2) 

"      w  y 


The  cost  sub-model  used  is  also  similar  to  the  example  discussed  in  the  previ- 
ous chapter.  Annual  maintenance  cost  AMCa,j{t)  is  assumed  to  be: 

AMCa,y{t)  =  Ka,yla,yt'"'''-'  (5.3) 

At  the  end  of  a  system's  life,  all  parts  that  can  be  salvaged  are  sold  and  those 
parts  that  could  not  be  sold  are  disposed  of  at  a  certain  cost.  Therefore,  the 
model  has  both  disposal  cost  and  salvage  value.  Operating  cost,  disposal  cost 
and  salvage  value  are  all  assumed  to  vary  exponentially  through  the  years, 
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Therefore,  the  total  life-cycle  cost  incurred  by  an  alternative  a  acquired  in 
year  y  and  used  for  t  years  is: 

+  i:{/<'a..7a.v<^--^+0Q,e*-'}  (5.4) 

t=\ 

C.     MODULES 

This  section  briefly  explains  the  function  and  workings  of  each  component 
module  of  the  prototype. 

1.  PREAMBLE 

This  module  is  peculiar  to  the  SIMSCRIPT  II. 5  language.  It  provides 
for  declaration  of  entities  and  global  variables.  Entities  are  structured  data  items 
that  represent  elements  of  a  model.  Definition  of  an  element  as  a  permanent  entity 
entails  the  automatic  creation  of  separate  arrays  for  each  declared  attribute  when  the 
entity  is  created,  the  length  of  these  arrays  being  the  number  of  that  element.  For 
example,  alternatives  were  defined  to  be  the  permanent  entity  ALT  with  attributes 
ALT. NAME  and  LEAD.TIME.  When  the  entity  ALT  is  created,  two  arrays  with 
the  names  of  the  attributes  will  be  created  simultaneously  and  their  length  would  be 
equal  to  the  number  of  alternatives  N.ALT  (a  system  variable).  Temporary  entities 
would  be  relevant  only  in  the  modeling  of  dynamic  systems  where  elements  enter 
and  leave. 

2.  MAIN 

As  the  name  suggests,  this  module  controls  the  overall  computational 
flow,  which  is  shown  graphically  in  Figure  5.2.  Note  that  the  modules  BEST. PATH 
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and  SUM.LCC.LCE  are  executed  several  times  if  the  Lagrangian  relaxation  ap- 
proach is  selected. 

3.  INITIALIZE 

This  module  reads  the  input  data,  controls  the  modules  generating  the 
system  renewal  network  and  creates  the  permanent  entities  and  associated  arrays.  It 
also  reserves  the  memory  space  needed  by  the  arrays.  The  flowchart  for  this  module 
is  shown  in  Figure  5.3.  The  input  data  is  read  in  free  format  in  the  order  shown  in 
the  listing. 

4.  CALCl.NODE 

The  structure  of  the  network  is  determined  in  this  module  based  on  the 
number  of  alternatives  and  the  minimum  and  maximum  life-span  of  each  system. 
The  number  of  permanent  entities  NODE  and  ARC  to  be  created  is  also  computed. 

5.  CALC2.NODE 

This  module  computes  for  each  node  information  such  as  the  year  and 
system  (transition)  it  belongs  to,  as  well  as  the  arcs  entering  it.  The  hierarchical  list 
data  structure  illustrated  in  Figure  5.4  was  used  to  represent  the  network.  Therefore, 
the  indices  of  the  arcs  entering  node  n  runs  from  EP{n)  to  [EP{n  +  1)  —  1]. 

6.  CALC.ALT.CODE 

The  code  of  the  alternative  represented  by  each  arc  is  assigned  in  this 
module. 

7.  CALC.ARC 

This  module  determines  the  tail  node,  life-cycle  cost  and  effectiveness 
associated  with  each  arc.  The  one-time  cost  items  in  Equation  5.4  are  discounted  and 
summed  separately  from  the  cumulative  items  such  as  maintenance  and  operating 
costs.   Acquisition  cost  of  the  present  system  is  ignored  as  sunk  cost.   The  benefit 
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accrued  in  each  year  is  added  to  yield  the  cumulative  or  life-cycle  elTectiveness  of 
the  arc  representing  a  given  alternative  and  life-span.  If  the  value  of  any  MOE 
falls  below  that  required  in  any  year  spanned  by  an  arc,  the  life-cycle  cost  of  that 
arc  is  set  to  an  unfavorably  high  value  to  mark  it  as  infeasible  and  to  prevent  its 
inclusion  in  the  selection  path.  The  highest  index  of  the  MOE(s)  that  failed  to  meet 
its  requirement(s)  is  recorded  as  the  limiting  MOE  in  such  cases.  The  flowchart  for 
this  module  is  shown  in  Figure  5.5. 
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(       START       ) 

CalllNITIAI.IZE 

•  Input  control  information 

•  Create  network  representing 
system  renewal  process 

•  Input  projected  cost  and 
effectiveness  parameters 

•  Compute  life-cycle  cost  and 
effectiveness  associated  with 
each  arc  in  the  network 


Call  BEST.PATH 
•    Uses  dynamic  program 
ming  to  solve  for  path 
that  is  optimal  under 
given  approach 


CallSUM  LCC  l-CK 
Compute  cumulative 
life-cycle  cost  and 
efTectivencss  along 
optimal  path 


CallLAGKANGIAN 
Determine  initial  upper 
and  lower  bounds  on  X 
Determine  A*  by  iterations 
of  the  supporting  plane 
method 

Note:  Rach  iteration 
entails  calling 
BKST.PATIIand 
SUM. LCCLCE  to  solve 
for  L(A)  and  slope  of  L(A). 


Call  RESULTS  PRINT 
Output  of  optimal  solution  for 
given  approach  based  on 
given  scenario  (and  budget  in 
the  case  of  the  Lagrangian 
relaxation  approach) 


END 


Figure  5-2:  Flowchart  or  MAIN  Program 
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Figure  5.3:  Flowchart  for  INITIALIZE  Module 
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Figure  5.4:  Hierarchical  List  Data  Structure  for  System  Renewal 

Network 
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START 


Initialize  life-cycle  eff  and 
limiting  MOE  to  zeto 


Compute  cost  &  effectiveness 
of  present  system  in  the 
current  year  as  well  as 
disposal  cost  and /or 
salvage  value 


Accumulatecost  and  eff  to 
obtain  cumulative  cost  and 
eff  for  various  posible 
residual  life  spans  of  the 
present  system 


Compute  cost  8  effectiveness 
of  alternatives  for  later 
systems  m  their  initial  year  of 
acquisition  as  well  as 
disposal  cost  and/or  salvage 
value 


Accumulatecost  &  eff  of 
alternatives  in  later  years  to 
obtain  the  life-cycle  cost  & 
eff  for  various  possible  life- 
spans of  the  alternatives 


includes  checks  of 
IVIOE  value  (which 
deteriorates)  in  each 
year  to  the  level 
required  in  that  year 
Failure  to  meet 
requirement  for  any 
MOE  in  any  year 
spanned  by  an  arc  results 
in  the  life  cycle  cost 
being  set  to  a  very  high 
value 


Return 


Figure  5.5:  Flowchart  for  CALC.ARC  Module 
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8.  ECHO.INPUT 

Key  inputs  may  be  echoed  for  debugging  purposes  by  setting  the  flag 
variable  ECHO. ON  in  the  input  file  to  "on".  This  module  will  then  be  activated 
and  intermediate  results  such  as  the  life-cycle  cost  and  effectiveness  associated  with 
each  arc  are  printed  along  with  the  key  input  data. 

9.  LAGRANGIAN 

This  module  produces  the  optimal  path  for  a  given  budget  using  the 
Lagrangian  relaxation  method.  The  value  of  A*  corresponding  to  the  optimal  path 
is  found  by  iterations  of  the  supporting  plane  method.  Each  of  these  iterations 
entails  calhng  the  modules  BEST.PATH  and  SUM.LCC.LCE  to  solve  for  L(A.)  and 
the  slope  of  L{X)  at  A,.  The  limitations  of  binary  representation  will  be  exceeded, 
even  with  double  precision  variables,  if  A,  and/or  the  cost  values  are  too  high.  The 
initial  upper  bound  Xmax  given  by  Equation  3.18  works  well  for  cost  values  scaled 
to  no  more  than  10^.  Figure  5.6  shows  the  flow  of  the  computations  in  this  module. 
Note  that  even  though  the  slope  at  the  minimum  of  L{X)  itself  may  be  positive  or 
negative  depending  on  whether  A3  =  A2  or  A3  =  Aj  respectively,  the  optimal  solution 
path  is  given  by  the  positively  sloped  line  segment  that  is  adjacent  to  the  minimum 
point.  When  A3  =  Ai,  this  line  segment  corresponds  to  A2  (otherwise,  the  intercept 
of  the  line  segments  at  Aj  and  A2  would  not  have  yielded  A3  =  Aj).  Therefore,  the 
slope  is  readjusted  when  A3  =  Ai. 

10.  BEST. PATH 

This  module  produces  the  optimal  solution  path  for  the  specified  approach 
using  forward  dynamic  programming.  The  criterion  function  L{J)  associated  with 
each  arc  J  when  the  approach  is  to  maximize  the  effectiveness  to  cost  ratio  of  the 
individual  systems  is  given  by: 
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_  LCEjJ) 

^   '  ~  {LCC{J)  -  MIN.LCC  +  1)  ^^■^' 


where: 


•  LCC{J)  and  LCE{J)  are  the  hfe-cycle  cost  and  efTectiveness  respectively  of 
arc  J 

•  MIN.LCC  =  min  LCC(J) 

This  is  to  allow  for  the  possibility  of  LCC{J)  being  zero  or  negative. 
Provisions  are  also  made  to  ensure  that  the  solution  path  has  sufficient  lead  time 
for  each  alternative  chosen.  The  computational  flow  is  illustrated  in  Figure  5.7. 
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Figure  5-6:  Flowchart  for  LAGRANGIAN  Module 
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Figure  5.7:  Flowchart  for  BEST.PATH  Module 
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11.  SUM.LCC.LCE 

The  solution  path  is  traced  backward  from  the  terminal  node  with  the 
highest  value  of  the  criterior  function  CRIT.FUNC(I)  and  the  life-cycle  cost  and 
effectiveness  of  the  systems  are  summed  along  the  way  in  this  module.  Note  that 
all  the  terminal  nodes  will  have  very  low  values  (negative)  of  CRIT.FUNC(I)  if  all 
paths  violate  either  lead  time  or  MOE  requirements  or  both,  i.e.,  no  feasible  solution 
exists.  Conversely,  if  CRIT.FUNC(I)  of  a  terminal  node  is  not  a  very  low  value, 
below  10~^,  then  a  feasible  solution  path  exists  that  leads  to  that  node. 

12.  RESULTS.PRINT 

Computational  results  are  printed  by  this  module. 

D.      COMPLEXITY  AND  EXECUTION  TIME 

In  terms  of  time  complexity,  the  dominating  piece  of  code  is  the  CALC.ARC 
module.  Let  the  difference  between  the  maximum  and  minimum  life-span  of  the 
present  system  be  rii  and  that  of  successor  systems  be  n2,n3,  •  •  •  ,71,,  where  s  is  the 
number  of  stages  or  systems  in  the  systems  renewal  network.  If  we  assume  n  —  ni 
for  Vi  >  1,  then  the  structure  of  the  network  is  given  by  Table  5.L 
The  total  number  of  arcs  n  is  thus  given  by: 

s-2 

n     =     1 +n,  -|-a(5  -  1)(1 -f  ni)(l +71) -|-an(l -f  n)^i 

=     1  -f  77,  +  a{s  -  1)(1  +  7ri)(l  4-  7i)  +  a77(l  +  n){s  -  2){s  -  l)/2      (5.6) 

If  we  let  the  average  life  of  successor  systems  be  denoted  as  /,  and  the  number  of 
MOEs  evaluated  be  e,  then  the  growth  rate  of  execution  time,  ET,  is  of  the  order 
7n(/  +  e).  In  other  words, 
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TABLE  5.1:   System  Renewal  Network  Structure  (ni  =  n  for  V,  >  1) 


Transition  # 

System  #i 

No.  of  Nodes 

No.  of  Arcs 

0 

1 

1 

(1  +  ni) 

1 

(1+nO 

2 

a(l  +  ni)(l4-n) 

2 

(1  +  7li  +  7l) 

3 

a(l  +  rii  +  7i)(l  +  n) 

3 

(l4-ni  +  2n) 

4 

a(l  -f  7li  +  27l)(l  +  70 

4 

(l  +  7ii  +  3n) 

And  so  on  until  i  =  s 

Note:  a  denotes  the  number  of  alternatives  evaluated 


ET    =     ^[m(/  +  e)] 


=    0 


14-n,  -\-a{s-  l)(l+n,)(l  +  n) + 


an^{s-2){s-  1) 


(/  +  e) 


=     0    max(an,ns,  a7i^5^)max(/,e) 


(5.7) 


Therefore,  execution  time  for  the  prototype  DSRPM  is  of  polynomial  order,  the  key 
determinants  of  which  are  n  and  s.  We  can  also  see  that  time  complexity  is  linear 
in  a,n,l  and  e.  For  a  run  with  a  =  2,ii\  =  5,n  =  5,s  =  3,1  =  7.5  and  e  =  3,  the 
execution  time  of  the  prototype  DSRPM  on  the  Zenith  AT-compatible  is  typically 
about  43  seconds. 
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VI.  CASE  STUDIES 

Case  studies  were  conducted  using  the  prototype  to  demonstrate  the  plausi- 
bility and  usefulness  of  the  DSRPM.  To  better  serve  the  purposes  of  illustration  and 
plausibility  testing,  the  case  studies  were  designed  to  comprise  a  series  of  scenarios 
of  incremental  complexity,  all  of  which  are  based  on  the  same  hypothetical  system 
discussed  in  Chapter  V.  Although  effort  was  made  to  incorporate  as  much  realism  as 
possible,  the  case  studies  and  the  weapon  system  were  kept  simple  to  better  portray 
the  workings  of  the  model.  A  good  understanding  of  and  confidence  in  the  model  is 
necessary  before  one  confronts  the  complexities  involved  in  a  real  application. 

In  the  case  studies,  the  key  effectiveness  parameters  or  MOEs  are  assumed  to 
be  range  (nautical  miles),  accuracy  (percentage  hit)  and  availability  (percentage) 
with  relative  weight  \Vj  of  0.4,  0.3  and  0.3  respectively.  The  cost  and  effectiveness 
characteristics  of  the  incumbent  system  assumed  are  shown  in  Table  C.l.  Disposal 
cost  is  assumed  to  be  zero  for  the  incumbent  system  and  alternatives,  i.e.,  DCa,y  =  0 
for  Va. 

TABLE  6.1:   Characteristics  of  the  Incumbent  System 


MOE 

Current  Value 

Po 

Cost  Parameters 

Range  (NM) 

104.6 

-0.03 

Maintenance 

Operating 

Salvage 

Accuracy  (%) 

87.2 

-0.02 

ko 

70 

OCo 

^0 

Slo 

oo 

Availability  (%) 

93.5 

-0.03 

30 

1.3 

7 

0.05 

-20 

-0.2 
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As  a  starting  point,  we  shall  assume  that  two  alternatives,  1  and  2,  are  available 
and  that  their  minimum  and  maximum  life  are  five  years  and  ten  years  respectively. 
They  are  assumed  to  require  no  development  or  acquisition  lead  time  (presumably 
available  "off-the-shelf").  The  incumbent  system  (alternative  0)  is  assumed  to  have  a 
minimum  and  maximum  life  of  zero  years  and  five  years  respectively,  i.e.,  alternatives 
may  be  introduced  immediately,  or  at  the  latest,  five  years  hence. 

Four  solution  approaches  are  used  in  each  of  the  case  studies,  namely  the 
minimum  average  annual  cost  (MAAC)  approach,  the  minimum  annualised  cost 
(MAC)  approach,  the  system  effectiveness  to  cost  ratio  (SECR)  approach  and  the 
Lagrangian  relaxation  (LR)  approach.  The  average  annual  cost  is  obtained  by  di- 
viding the  total  cost  of  a  succession  of  three  systems  (in  current  dollars)  by  the  years 
spanned  by  the  systems.  The  MAAC  approach  thus  ignores  monetary  interest.  The 
MAC  solution  is  obtained  by  amortizing  the  total  cost  of  the  system  over  the  years 
spanned.  It  is  the  solution  obtained  when  one  wishes  to  minimise  cost  needed  to 
meet  effectiveness  requirements.  The  SECR  approach  maximizes  the  cumulative 
total  of  the  effectiveness- to-cost  ratio  of  each  of  the  three  systems  in  a  path.  The 
LR  approach  was  explained  in  Chapter  3.  Solution  points  obtained  using  the  LR 
approach  are  finked  as  in  Figure  3.5  to  form  the  efficiency  boundary.  The  point  on 
the  boundary  that  yields  the  maximum  overall  effectiveness  to  cost  ratio  would  be 
the  preferred  solution  if  one  adopts  the  approach  of  maximising  returns  on  cost. 

A  notational  point  to  note  is  that  the  transition  years  shown  in  the  results  of 
the  case  studies  are  all  abbreviated  to  the  last  two  digits.  That  is,  97  and  05  mean 
1997  and  2005  respectively.  The  current  year  is  1989. 


A.      CASE  STUDY  1:  STATIONARY  ENVIRONMENT 

The  first  scenario  we  examine  is  one  in  which  the  only  change  occurring  is 
that  due  to  deterioration.  Technology  and  requirement  levels  are  at  a  standstill  and 
no  inflation  occurs.  The  value  of  money  is  constant  over  time,  i.e.,  interest  rate  is 
zero.  This  highly  improbable  scenario  is  implied  by  the  commonly  used  "equal  life" 
assumption  whereby  an  alternative,  once  selected,  will  be  selected  over  and  over 
again  and  will  always  have  the  same  service  life. 

Supposing  alternatives  1  and  2  have  the  following  cost  and  effectiveness  char- 
acteristics shown  in  Table  6.2  and  6.3.  Alternative  2  is  thus  three  times  as  expensive 
to  acquire  but  has  a  longer  range  and  offers  lower  maintenance  and  operating  costs 
than  Alternative  1.  Requirement  levels  are  assumed  constant  at  95  nm  for  range, 
80%  for  accuracy  and  70%  for  availability.  Results  of  analysis  performed  using  the 
prototype  DSRPM  for  this  case  is  presented  in  Table  6.4  below. 

TABLE  6.2:  Characteristics  of  Alternative  1 


MOE 

Value  \\"heu  New 

P\ 

Cost  Parameters 

Ilauge(nm) 

120 

-0.03 

Acquisition 

Mainteuauce 

Operating 

Salvage 

Accuracy  (%) 

95 

-0.02 

AC, 

h 

71 

OCx        (^1 

5V, 

o\ 

!  Availability  (7c)                  95 

-0.03 

100 

30 

1.3 

7    0.05 

-50    -0.2  1 

The  value  of  the  total  life-cycle  cost  and  effectiveness  are  plotted  in  Figure  6.1. 
As  interest  rate  is  zero,  MAAC  and  MAC  approaches  yielded  the  same  solution.  The 
results  show  that  an  alternative  that  is  selected  will  be  selected  repeatedly  in  the 
future,  the  service  life  each  time  being  the  same.  This  is  only  to  be  expected  as, 
urKJcr  stat ic  conditions,  tiicre  arc  no  reasons  for  change.  Note  that  the  solution  that 
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TABLE  6.3:  Characteristics  of  Alternative  2 


MOE 

Value  When  New 

P7 

Cost  Parameters 

Range  (run) 

140 

-0.03 

Acquisition 

Maintenance 

Operating 

Salvage 

Accuracy  (%) 

95 

-0.02 

ACi 

^1 

72 

OC2 

<^2 

s\\ 

a-i    \ 

Availability  (%) 

95 

-0.03 

300 

15 

13 

4 

0.05 

-90     -0.2 

TABLE  6.4:  Results  for  Stationary  Environment 


Solution 

Total  LCC 

Total  LCE 

LCE/LCC 

A  t.  Chosen 

Transition  Year 

K 

Approach 

(3  Systems) 

(3  Systems) 

Ratio 

1 

2 

3 

1 

2 

3 

MAAC 

1206.05 

16.550 

0.01372 

0 

2 

2 

91 

99 

07 

MAC 

1206.05 

16.550 

0.01372 

0 

2 

2 

91 

99 

07 

SECR 

1097.35 

15.062 

0.01373 

0 

2 

2 

89 

97 

05 

LR 

1270.84 

17.265 

0.01359 

0 

2 

2 

89 

00 

08 

0 

LR 

1200.05 

16.550 

0.01372 

0 

2 

2 

91 

99 

07 

0.01104 

LR 

1147.33 

15.816 

*0.01379 

0 

2 

2 

90 

98 

06 

0.01251 

LR 

801.84 

10.934 

0.01364 

0 

1 

1 

90 

95 

00 

0.01413 

LR 

751.86 

10.179 

0.01354 

0 

1 

1 

89 

94 

99 

0.01509 

gives  the  maximum  return  per  unit  cost  for  the  three  systems  overall  (marked  with 
an  asterisk  in  Table  6.4  and  a  square  in  Figure  6.1)  is  different  from  that  obtained 
using  the  SECR  approach.  The  overall  LCE  to  OCC  ration  in  the  case  of  the  SECR 
solution  is  not  as  high  even  though  the  solution  also  lies  on  the  efTicicncy  boundary. 
The  implication  of  this  is  obvious  -  maximizing  return  per  unit  cost  for  individual 
systems  considered  separately  leads  to  a  good  but  not  best  solution  in  terms  of 
overall  return  per  unit  cost. 
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Figure  6.1:  Stationary  Environment  Case  Study  Results 

B.     CASE  STUDY  2:  INFLATION  AND  INTEREST  RATES 

The  next  scenario  we  look  at  is  probably  the  one  most  commonly  encountered 
in  a  typical  system  renewal  study.  Future  inflation  and  interest  rates  are  assumed 
to  be  non-zero  and  constant,  and  future  technological  advancements  and  changes 
in  requirement  levels  are  ignored.  For  illustrative  purposes,  let  us  say  the  inflation 
and  interest  rates  were  projected  to  be  constant  at  four  and  five  percent  per  annum 
respectively.  Cost  projections  based  on  these  rates  are  shown  in  Figure  6.2  in  the 
sample  input  file  for  this  case  study  (Appendix  B).  Analysis  based  on  this  scenario 
produced  the  results  shown  in  Table  6.5. 
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Figure  6  2a:  Cost  Evolution  for  New  Alternative  1 
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Figure  6  2b:  Cost  Evolution  for  New  Alternative  2 

Figure  6.2  :  Evolution  of  Cost  Parameters  for  Alternatives  1  and  2 
With  Inflation  (  No  Technological  Progress  ) 


79 


TABLE  6.5:  Results  for  Inflation  and  Interest  Rates 


Solution 

Total  LCC 

Total  LCE 

LCE/LCC 

Alt. 

Chosen 

Transit 

on  \'ear 

a: 

Approach 

3  Systems 

3  Systems 

Ratio 

1 

2 

3 

1 

2 

3 

MAAC 

1031.26 

17.266 

0.01674 

0 

2 

2 

92 

00 

08 

MAC 

811.98 

13.375 

0.01647 

0 

2 

1 

90 

98 

03 

SECR 

921.03 

14.396 

0.01563 

0 

2 

2 

91 

97 

04 

LR 

1031.26 

17.265 

0.01674 

0 

2 

2 

92 

00 

OS 

0 

LR 

986.44 

16.550 

•0.01678 

0 

2 

2 

91 

99 

07 

0.01596 

LR 

944.29 

15.816 

0.01675 

0 

2 

2 

90 

98 

06 

0.01742 

LR 

811.98 

13.375 

0.01647 

0 

2 

1 

90 

98 

03 

0.01845 

LR 

683.34 

10.934 

0.01600 

0 

1 

1 

90 

95 

00 

0.01898 

LR 

643.73 

10.179 

0.01581 

0 

1 

1 

89 

94 

99 

0.01904 

Total  LCC  and  LCE  from  this  table  are  plotted  in  Figure  6.3.  The  efficiency 
boundary  is  shifted  towards  lower  cost  since  the  interest  rate  is  higher  than  the 
inflation  rate.  An  opposite  movement  would  occur  if  the  converse  was  true.  Note 
that  the  shift  is  greater  for  the  higher  cost  end  because  inflation  and  interest  rates 
have  a  bigger  impact  there.  As  the  interest  rate  is  significantly  non-zero,  we  see 
that  the  solutions  obtained  using  the  MAAC  and  MAC  approaches  are  different. 
Note  also  that  the  SECR  solution  is  clearly  inferior  as  it  lies  somewhat  below  the 
efficiency  boundary.  From  the  above  table,  we  can  see  that  Alternative  1  is  the 
low-end  alternative.  The  last  line  in  the  table  shows  the  solution  that  meets  the 
specified  minimum  requirements  at  the  lowest  possible  total  cost,  which  is  the  one 
where  the  service  lives  of  all  three  systems  are  at  their  minimums. 
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Figure  6.3:  Inflation  and  Interest  Rate  Case  Study  Results 

C.      CASE  STUDY  3:  TECHNOLOGICAL  ADVANCEMENT 

For  this  scenario,  let  us  suppose  that  technological  progress  has  the  elTect  of 
suppressing  cost  growth  due  to  inflation  from  four  to  two  percent  per  annum.  The 
revised  cost  projections  are  shown  in  Figure  6.4.  The  sample  input  file  for  tliis  case 
study  is  shown  in  Appendix  B.  In  addition,  technological  advances  are  also  reflected 
in  improvement  in  performance  characteristics.  For  the  purpose  of  this  case  study, 
let  us  suppose  that  the  various  MOEs  evolve  as  shown  in  Figure  6.5  and  the  sample 
input  file  in  Appendix  B.  Improvements  are  primarily  in  range,  with  Alternative  1 
becoming  marginally  better  than  Alternative  2  about  ten  years  from  now.  Other 
conditions  are  as  assumed  in  Case  Study  2.  The  results  for  this  case  study  are  shown 
in  Figure  6.6  and  Table  6.6. 
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Contrasting  Table  6.5  and  Table  6.6,  we  note  that  Alternative  1  is  now  the 
dominantly  favored  alternative  because  it  became  competitive  performance-wise  as 
well  as  cost-wise.  Technological  advancement  resulted  not  only  in  a  further  shift 
in  the  efficiency  boundary  towards  lower  cost  but  also  an  extension  into  higher 
levels  of  cumulative  effectiveness.  The  lower  end  of  the  boundary  is  not  affected 
as  much  as  solutions  in  the  higher  end  that  extends  further  into  the  future.  Note 
that  technological  advancement  also  results  in  lengthening  of  the  maximum  service 
life  (three  systems),  a^  evident  in  the  maximum  effectiveness  solution  identified  by 
X*  =  0.  In  Case  Studies  1  and  2,  no  solutions  extend  beyond  the  year  2008  whereas 
in  this  case  study,  some  solutions  extend  to  the  year  2009. 

TABLE  6.6:  Results  for  Technological  Advancement 


Solutiou 

Total  LCC 

Total  LCE 

LCE/LCC 

Alt. 

Choseu 

Transit 

ion  Year 

Ar 

1  Approach 

3  Systems 

3  System.'? 

Ratio 

1 

2 

3 

1 

2 

3 

MAAC 

888.11 

18.491 

0.02082 

0 

1 

1 

92 

00 

09 

MAC 

736.01 

14.509 

0.01971 

0 

1 

1 

91 

99 

04 

SECR 

804.94 

16.220 

0.02015 

0 

1 

1 

92 

99 

06 

LR 

907.83 

18.775 

0.02068 

0 

2 

2 

92 

00 

09 

0 

LR 

888.11 

18.491 

+0.02082 

0 

1 

1 

92 

00 

09 

0.01444 

LR 

634.53 

12.219 

0.01926 

0 

1 

1 

91 

96 

01 

0.02473 

LR 

601.75 

11.306 

0.01889 

0 

1 

1 

90 

95 

00 

0.02601 

LR 

574.70 

10.506 

0.01828 

0 

1 

1 

89 

94 

99 

0.03180 
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Figure  6.4  :  Evolution  of  Cost  Parameters  for  Alternatives  1  and  2 
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Figure  6.6:  Teclinological  Advancement  Case  Study  Results 
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D.      CASE  STUDY  4:  REQUIREMENT  ESCALATION 

Let  us  suppose  that  everything  is  as  before,  except  that  requirement  levels  set 
for  range  and  availability  increase  over  time  as  shown  in  Figure  6.7.  The  results  for 
this  scenario  are  shown  in  Table  6.7  and  illustrated  in  Figure  6.8. 

TABLE  6.7:  Results  for  Requirement  Escalation 


Solution 

Total  LCC 

Total  LCE 

LCE/LCC 

At. 

Chosen 

Transit 

on  Year 

A* 

Approach 

3  Systems 

3  Systems 

Ratio 

1 

2 

3 

1 

2 

3 

MAAC 

741.62 

14.498 

0.01955 

0 

1 

91 

96 

04 

MAC 

741.C2 

14.498 

0.01955 

0 

1 

91 

96 

04 

SECR 

771.83 

14.865 

0.01926 

0 

2 

91 

98 

04 

LR 

872.43 

16.431 

0.01883 

0 

2 

91 

98 

06 

0 

LR 

835.93 

16.395 

•0.01961 

0 

2 

91 

98 

06 

0.00097 

LR 

741.G2 

14.498 

0.01955 

0 

1 

91 

96 

04 

0.02011 

LR 

034.53 

12.219 

0.01926 

0 

1 

91 

96 

01 

0.02129 

LR 

574.70 

10.506 

0.01828 

0 

1 

89 

94 

99 

0.02863 

Comparing  Figures  6.6  and  6.8,  we  can  see  that  the  escalation  in  requirement 
levels  resulted  in  a  downward  depression  of  the  higher  cost  end  of  the  efficiency 
boundary.  This  happens  because  the  escalation  in  requirement  levels  results  in  the 
shortening  of  service  lives  only  if  the  solution  path  extends  far  enough  to  be  affected. 
Note  that  the  maximum  effectiveness  solution,  for  instance,  now  extends  only  until 
the  year  2006  whereas  in  the  previous  case  study  it  reaches  the  year  2009.  An 
interesting  result  in  the  equality  of  the  MAAC  and  MAC  solutions  even  though 
the  interest  rate  is  five  percent  per  annum.  Note  also  that  solutions  calling  for  the 
first  transition  to  be  made  in  1991  are  clearly  dominating. 
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Figure  6.7:  Escalation  in  Requirement  Levels 
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Figure  6.8:   Requirement  Escalation  Case  Study  Results 
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E.      CASE  STUDY  5:  SCHEDULING  OF  UPGRADES 

Supposing  Alternative  1  is  an  upgrading  option,  the  exercise  of  which  is  sig- 
nificantly cheaper  if  carried  out  along  with  major  scheduled  overhauls.  If  such 
maintenance  is  carried  out  every  three  years,  (e.g.,  refit  of  a  ship),  and  the  savings 
is  say  twenty  percent,  then  the  acquisition  cost  of  a  new  Alternative  1  varies  as 
illustrated  in  Figure  6.9.  In  this  scenario,  everything  else  is  a^  assumed  for  Case 
Study  4. 

The  results  for  this  scenario  are  given  in  Table  6.8  and  Figure  6.10.  As  evident 
from  the  table,  there  are  no  saving  realised  in  the  higher  cost  end  where  the  discounts 
are  not  taken  up.  The  lower  end  of  the  efficiency  boundary  shifts  to  the  left  as  a 
result  of  the  savings.  We  see  that  the  solution  offering  maximum  effectiveness  per 
unit  cost  suggests  that  the  incumbent  should  be  replaced  with  Alternative  2  in 
1990,  then  upgraded  with  Alternative  1  in  1997  (a  non-maintenance  year  with  no 
discount)  for  use  until  the  year  2005.  The  SECIl  solution  also  indicates  an  upgrade 
carried  out  in  a  non-maintenace  year  (1998).  The  overall  desirability  of  matching  a 
major  upgrading  program  to  some  maintenance  schedule  to  minimise  downtime  or 
to  secure  some  savings  can  thus  be  evaluated  using  the  model. 
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Figure  6.10:  Scheduling  of  Upgrading  Case  Study  Results 
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TABLE  6.8:  Results  for  Upgrade  Scheduling 


Solution 

Total  LCC 

Total  LCE 

LCE/LCC 

Alt. 

Chosen 

Transit 

on  Year 

A* 

Approach 

3  Systems 

3  Systems 

Ratio 

1 

2 

3 

1 

2 

3 

MA  AC 

741.62 

14.498 

0.01955 

0 

1 

1 

91 

96 

04 

MAC 

653.39 

12.734 

0.01949 

0 

1 

1 

89 

94 

02 

SECR 

771.83 

14.865 

0.01926 

0 

2 

1 

91 

98 

04 

LR 

872.43 

16.431 

0.01883 

0 

2 

2 

91 

98 

06 

0 

LR 

835.93 

16.395 

0.01961 

0 

2 

1 

91 

98 

06 

0.00097 

LR 

791.73 

15.532 

•0.01962 

0 

2 

1 

90 

97 

05 

0.01953 

LR 

538.06 

10.506 

0.01953 

0 

1 

1 

89 

94 

99 

0.01981 

Note:  In  this  case  study,  a  1  for  alternative  chosen  represents  an  upgrade  of  the 
preceding  system.  Thus  the  sequence  0,2,1  means  the  incumbent  is  replaced  with 
alternative  2  and  then  upgraded  with  alternative  1. 

F.      CASE  STUDY  6:  DEVELOPMENT  LEAD  TIME 

Let  us  suppose  that  everything  is  as  in  Case  Study  5  except  that  Alternative 
1  is  an  upgrade  that  requires  one  year  of  development  and  cannot  be  introduced 
immediately.  The  impact  of  such  lead-time  requirements  can  be  seen  in  Table  6.9 
and  Figure  6.11.  From  the  table,  we  can  see  that  the  lowest  cost  solution  (bottom  of 
table)  still  involves  an  immediate  transition  as  in  the  other  case  studies.  However, 
the  transition  is  to  Alternative  2  in  this  case  as  the  model  now  bars  the  immediate 
execution  of  Alternative  1.  The  diagram  clearly  shows  that  the  lead-time  require- 
ment pushes  up  cost  at  the  lower  cost  end  of  the  efficiency  boundary.  The  MAAC 
solution,  the  SECR  solution  and  the  maximum  effectiveness  per  unit  cost  solution 
are  all  unaffected.  The  MAC  solution,  however,  became  equal  to  the  MAAC  solu- 
tion (curiously).  Note  that  the  MAC  solution  in  Case  Study  5  involved  immediate 
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introduction  of  Alternative  1  (upgrade  of  alternative  0  preceding  it);  a  solution  that 
is  not  admissible  in  this  case  study. 

TABLE  6,9:  Results  for  Development  Lead-Time 


Solution 

Total  LCC 

Total  LCE 

LCE/LCC 

Alt.  Chosen 

Transition  "\'ear 

a; 

Approach 

3  Systems 

3  Systems 

Ratio 

1 

2 

3 

1 

2 

3 

MAAC 

741.62 

14.498 

0.01955 

0 

1 

91 

96 

04 

MAC 

741.62 

14.498 

0.01955 

0 

1 

91 

96 

04 

SECR 

771.83 

14.865 

0.1926 

0 

2 

91 

98 

02 

LR 

872.43 

16.431 

0.01883 

0 

2 

91 

98 

06 

0 

LR 

835.93 

16.395 

0.01961 

0 

2 

91 

98 

06 

0.00097 

LR 

791.73 

15.532 

•0.019G2 

0 

2 

90 

97 

05 

0.01953 

LR 

741.62 

14.498 

0.01907 

0 

1 

91 

96 

04 

0.02063 

LR 

634.53 

12.219 

0.01926 

0 

1 

91 

96 

01 

0.02129 

LR 

621.31 

10.754 

0.01731 

0 

2 

89 

94 

99 

0.11078 
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Figure  6.11:   Development  Lead  Time  Case  Study  Results 

G.      GENERAL  OBSERVATIONS 

From  the  series  of  six  case  studies  analysed  using  the  prototype  DSRPM,  we 
saw  the  impact  of  various  issues  as  they  were  integrated  into  the  model.  DSRPM  of- 
fers a  common  framework  for  decision  makers  to  explore  the  consequences  of  various 
possibilities.  It  does  not  advocate  any  single  solution  for  a  given  scenario,  but  rather 
maps  out  the  efficiency  boundary  spanning  the  entire  solution  space  and  indicates 
the  solutions  obtained  using  various  approaches.  A  plot  of  the  results  for  various 
scenarios,  as  is  done  in  Figure  6.12  for  the  case  studies,  gives  the  decision  maker  a 
good  overview  of  the  problem.  The  picture  can  be  enriched  by  varying  the  shade  of 
the  result  for  each  scenario  according  to  its  likelihood  of  occurrence.  Furthermore, 
a  third  axis  indicating  overall  service  time  for  the  three  systems  can  be  added  as 


91 


in  Figure  6.13  to  give  fuller  expression  to  the  problem.  As  explained  in  Chapter  2, 
the  system  renewal  problem  is  basically  a  three-dimensional  one  -  cost,  effectiveness 
and  time. 
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Figure  6.12;  Consolidated  Results  for  Case  Studies 


Although  three  systems  are  analysed  in  the  model,  the  focus  is  really  on  when 
the  first  transition  should  take  place  and  what  alternative  should  be  selected  to 
succeed  the  incumbent  system.  Anal^'sis  may  reveal  a  dominating  (good  under 
all  scenarios)  or  robust  (favored  under  most  scenarios)  solution  to  the  timing  and 
choice  involved  in  the  first  transition.  For  example,  the  escalation  of  requirement  in 
Case  Study  4  resulted  in  hastening  of  the  first  transition  by  a  year  for  the  maximum 
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Figure  6.13:  3-D  Plot  of  Efficiency  Boundaries 

overall  effectiveness  to  cost  ratio  solution.  In  Case  Study  6,  however,  we  saw  that  this 
solution  was  not  affected  by  the  development  lead  time  of  one  year.  By  varying  the 
scenarios  for  various  requirement  escalation  trends,  development  lead  time  estimates, 
etc.,  we  can  test  the  sensitivity  of  various  solutions  to  variations  in  estimates  and 
trends.  Even  if  analysis  proved  inconclusive,  a  rough  idea  of  how  long  the  incumbent 
is  likely  to  be  kept  will  emerge.  No  commitment  needs  to  be  made  if  the  time  to 
the  likely  date  of  the  first  transition  is  more  than  the  lead  time  required  to  effect 
the  likely  transition.  For  example,  if  analysis  shows  that  an  EW  suite  upgrade  of  a 
frigate  should  take  place  three  years  down  the  road  and  such  an  upgrade  has  a  lead 
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time  of  two  years,  we  needn't  make  a  commitment  yet,  but  we  know  the  decision 
is  due  a  year  from  now.  At  that  time,  tlie  model  can  be  run  again  with  updated 
information.  In  this  way,  we  can  ensure  decisions  are  made  in  a  timely,  fashion  using 
the  latest  information. 

Besides  facilitating  timely,  consistent  and  rational  decision-making,  the  model 
also  looks  forward  to  the  future  to  provide  estimates  of  cost,  effectiveness  and  service 
life  that  facilitates  good  planning.  If  aggregated  to  sufficiently  high  levels,  the  model 
allows  future  year  procurement  budgets  to  be  forecasted.  Designs  and  support  plans 
for  succeeding  systems  can  be  optimised  to  service  lives  that  are  more  realistic 
in  terms  of  broad  economic,  technological  and  operational  outlook.  Amortization 
periods  would  no  longer  need  to  be  arbitrarily  assumed  and  budgets  can  be  better 
planned  knowing  cost  estimates  of  succeeding  systems.  The  anticipated  system 
selections  also  serves  to  better  guide  planning  for  defense  infrastructure,  R&D  and 
other  strategic  concerns.  The  explicitness  of  the  model  fosters  better  consensus 
building  and  coordination  among  the  various  agencies  responsible  for  system  renewal 
planning.  As  discussed  in  Chapter  4,  Section  E,  force  level  and  mix  issues  as  well 
as  phased  build-up  versus  bulk  procurement  decisions  can  also  be  dealt  with  by 
the  DSRPM.  The  model  thus  provides  a  comprehensive  and  systematic  framework 
for  system  renewal  decisions  to  be  made  optimally  and  in  a  consistent  and  timely 
manner  based  on  current  forecasts. 

Lastly,  it  should  be  noted  that  inflation  and  interest  rates  were  assumed  con- 
stant and  linear  utility  functions  were  used  in  the  case  studies  for  reason  of  conve- 
nience and  not  necessity.  The  model  allows  these  rates  and  effectiveness  requirement 
levels  to  be  varied  for  future  years  in  a  scenario.  Long-range  projections  can  be 
varied  more  than  near-term  ones  to  account  for  greater  uncertainty  in  these  projec- 
tions.  The  impact  of  technological  advancement  is  also  freely  interpreted  in  terms 
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of  changes  in  cost  and  effectiveness  parameter  projections  and  no  restrictions  are 
placed  on  the  definition  of  utiHty  functions.  The  DSRPM  thus  affords  very  flexible 
means  to  model  the  factors  influencing  systems  renewal  problems. 
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VII.  CONCLUSIONS  AND 
RECOMMENDATIONS 

A.      CONCLUSIONS 

The  dynamic  system  renewal  planning  model  (DSRPM)  provides  a  system- 
atic framework  within  which  optimal  system  renewal  decisions  may  be  made  in  a 
consistent  and  timely  manner,  under  realistic  conditions  of  changing  environments. 
Conceptually,  it  embodies  a  concept  of  obsolescence  that  is  more  comprehensive, 
better  formalised  and  more  relevant  to  military  planners  than  is  available  in  current 
literature.  Prediction  of  what  future  renewals  are  likely  to  be  made  and  when,  and 
estimation  of  the  associated  cost  and  effectiveness  would  allow  plans  to  be  made 
anticipatively  with  a  long-term  perspective.  Through  simple  reformulation  and  use 
of  shadow  prices,  bulk  versus  phased  procurement  decisions  as  well  as  force  level 
and  mix  issues  can  also  be  examined  using  the  model. 

The  plausibility  and  usefulness  of  the  model  was  demonstrated  by  a  scries  of 
case  studies  using  the  prototype  that  was  developed.  The  equal  life  assumption 
often  used  in  system  renewal  models  was  illustrated  by  the  first  case  study.  Other 
case  studies  showed  the  impact  of  interest  and  inflation  rates,  technological  progress, 
requirement  escalation  and  development  lead  time.  Application  of  the  model  to  the 
scheduling  of  upgrades  that  are  cheaper  to  carry  out  during  overhaul  periods  was 
also  shown. 

The  prototype  DSRPM  illustrates  how  the  model  may  be  used  as  a  descriptive 
decision  aid  to  better  understand  the  complex  effects  of  a  diverse  range  of  issues  on 
the  system  renewal  decision.  It  was  shown  how  the  model  organises  financial,  oper- 
ational, intelligence  and  other  staff  j)rojcctions  within  an  integrated  framework  to 
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provide  decision  makers  with  a  broad,  consistent,  long-term  perspective  of  relevant 
issues.  The  decision  maker  is  presented  with  a  graphic  picture  of  the  entire  solu- 
tion space  as  structured  by  the  system  renewal  scenarios  considered.  Scenarios  are 
flexibly  modelled  and  are  used  to  point  to  solutions  that  are  robust  in  the  face  of  un- 
certainty. By  continually  refining  scenarios  based  on  lessons  learned  from  exercises, 
audits,  and  other  new  information  on  trends,  past  experience  is  brought  to  bear  on 
the  present  so  that  the  basis  of  decisions  is  kept  current  and  realistic.  The  commu- 
nication among  all  parties  involved  in  this  process  is  facilitated  by  the  explicitness 
of  the  model.  This  in  turn  fosters  consensus  building  necessary  for  cohesive  effort 
in  defense  planning. 

The  model  that  was  developed  is  a  conceptual  prototype  that  is  as  yet  untested 
by  a  real  application.  Much  remains  to  be  done  before  it  is  ready  for  organisation- 
wide  implementation.  However,  it  holds  much  promise  and  should  be  examined 
further  for  future  implementation. 

B.     RECOMMENDATIONS 

The  prototype  that  was  developed  is  a  rudimentary  one  that  could  be  extended 
in  many  ways  as  future  research  topics.  It  is  recommended  that  research  proceed  in 
the  following  order  of  tractability: 

1.  Factoring  foreign  exchange  rates  into  the  model  is  a  desirable  goal.  The  treat- 
ment of  interest  and  inflation  rates  would,  however,  require  some  adjustments. 

2.  Optimal  force  mix  can  be  determined  if  shadow  prices  of  different  components 
of  a  task  force  are  made  compatible.  Means  of  integrating  the  utility  functions 
of  these  components  under  a  common  mission  structure  needs  to  be  examined. 
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3.  Simulation  and  forecasting  capabilities  can  be  incorporated  to  better  explore 
the  range  of  system  renewal  scenarios  that  may  be  encountered  in  the  future. 

4.  Stochastic  models  can  be  constructed  to  better  portray  uncertainty  in  cost, 
effectiveness  and  requirement  levels. 

5.  Utility  modelling  may  be  facilitated  by  the  use  of  implicit  approaches  that  do 
not  require  utility  functions  to  be  explicitly  defined  over  the  entire  domain  of 
the  MOE  value. 

6.  System  renewal  models  developed  for  various  systems  could  be  linked  and  ag- 
gregated to  higher  entities  using  schemes  that  incorporate  organisation  struc- 
ture into  the  solution  algorithm  to  avoid  sub-optimisation.  Plans  developed 
for  ships  could,  for  example,  be  tied  to  plans  for  their  operational  task  forces 
and  administrative  units.  The  ideal  would  be  to  achieve  a  level  of  aggregation 
high  enough  to  interface  with  arms  race  and  econometric  models  so  that  very 
high  levels  of  optimisation  may  be  achieved. 

7.  Fuzzy  set  and  logic  theory  may  be  introduced  to  better  model  elements  that 
are  inherently  vague.  Requirement  levels,  for  example,  may  be  made  elastic  by 
using  what  is  known  as  fuzzy  numbers.  This  may  allow  some  form  of  abstract 
qualitative  reasoning  to  be  represented  in  the  model. 
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APPENDIX  A:  PROGRAM  LISTINGS 


1  preamble 
2 

3  normally,  mode  is  real,  dimension  is  0 
4 

5  permanent  entities 
6 

7        every  YEAR  has  a  DISCOUNT. FACTOR  and  a  CUM.DF 
8 

9  normally,  mode  is  integer 
10 

11  every  NODE  has  a  TRANS. YEAR,  a  SYS. CODE, 

12  a  BEST. PRED. NODE  and  a  BEST. PRED. ARC 

13  every  ARC  has  a  TAIL  and  an  ALT. CODE 

14  every  SYS  has  a  MIN.LIFE,  a  MAX. LIFE,  a  NUM. NODES,  a  NUM. ARCS  and 

15  a  BASE. YEAR 

16  every  ALT  has  an  ALT. NAME  and  a  LEAD. TIME 
17 

18  define  MOE,  MOE.DF  as  3-dimensional  real  arrays 

19  define  AC,  MC.Pl,  MC.P2,  OC,  OC.EF,  DC,  DC.EF,  SV,  SV.EF  and 

20  MIN.MOE.REQ  as  2-dimensional  real  arrays 

21  define  MOE. NAME  as  a  1-dimensional  text  array 

22  define  L  and  CRIT.FUNC  as  1-dimensional  double  arrays 

23  define  LCE,  LCC,  MOE.WT,  O.MOE  and  O. MOE.DF  as  1-dimensional 

24  real  arrays 

25  define  MAX . PRED. NODES ,  LIM.MOE  and  EP  as  1-dimensional  integer  arrays 

26  define  ARGMAX.CF,  ARC. MAC,  ARG.AAC,  CURR.YEAR,  NUM. MOE  and  NUM. ALT  as 

27  integer  variables 

28  define  LAMDA,  LAMDAl,  LAM DA 2 ,  LAM DA 3 ,  SUM. LCC,  SUM . LCC . 1 ,  SUM. LCC. 2, 

29  SUM. LCE,  SUM. LCE. 1,  SUM. LCE. 2,  MAX. LCE,  MIN.LCC,  MOE. VALUE, 

30  SLOPEl,  SLOPE2,  SL0PE3,  LI,  L2 ,  L3 ,  BUDGET,  AMC,  AOC,  DCY , 

31  SVY,  O. MC.Pl,  0.MC.P2,  O.OC,  O. OC.EF,  O.DC,  O. DC.EF,  O.SV, 

32  O. SV.EF,  MAX. MAC,  MAX.AAC,  AAC.CF,  INT. RATE  and  INT.DF  as 

33  real  variables 

34  define  MAX. OF  as  a  double  variable   •'  Double  means  double  precision. 

35  define  APPROACH,  ECHO. ON,  CASE. NAME  and  ALT. NAME  as  text  variables 
36 

37  end 
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1  main 

2 

3     call  INITIALIZE 

4 

5  if  APPROACH  =  "Lagrangian" 

6  call  LAGRANGIAN 

7  else 

8  call  BEST. PATH 

9  call  SUM.LCC.LCE 
10    always 

11 

12     call  RESULTS. PRINT 

13 

14     stop 

15 

16  end 
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1  routine  INITIALIZE 
2 

3  open  1  for  input,  name  Is  "IN. DAT" 

4  use  1  for  input 

5  open  2  for  output,  name  is  "OUT. DAT" 

6  use  2  for  output 
7 

8  read  CASE. NAME 

9  read  N.SYS,  N.ALT,  NUM.MOE,  CURR.YEAR,  INT. RATE 
10  read  ECHO. ON 

11 

12  create  every  SYS  and  ALT 

13 

14  N.YEAR  =  1 

15  for  I  =  1  to  N.SYS 

16  do 

17  read  MIN.LIFE(I)  and  MAX.LIFE(I)   '•  Computing  number  of  years  to 

18  N.YEAR  =  N.YEAR  +  MAX.LIFE(I)      '•  be  covered  in  the  analysis. 

19  loop 
20 

21  create  every  YEAR 
22 

23  for  I  =  1  to  N.ALT 

24  read  ALT. NAME (I)  and  LEAD. TIME (I) 

25  read  DISCOUNT. FACTOR( 1) 

2  6  CUM.DF(l)  =  DISCOUNT. FACTOR (1) 

27  for  I  =  2  to  N.YEAR 

28  do 

29  read  DISCOUNT. FACTOR (I) 

30  CUM.DF(I)  =  CUM.DF(I-l)  *  DISCOUNT. FACTOR ( I ) 

31  loop 

32  read  APPROACH      • '  Determines  the  type  of  solution  to  be  generated 

33  read  BUDGET        ' •  Used  only  in  the  Lagrangian  approach 
34 

35  ''  Allocating  memory  to  arrays 

36  reserve  MAX . PRED. NODES ( *)  as  N.SYS 

37  reserve  AC (*,*) ,  MC.P1(*,*),  MC.P2(*,*),  OC(*,*),  OC.EF(*,*), 

38  DC(*,*),  DC.EF(*,*),  SV ( * , * ) ,  SV.EF(*,*)  as  N.ALT  by  N.YEAR 

39  reserve  MOE(*,*,*)  and  MOE. DF(* , * , * )  as  N.ALT  by  NUM.MOE  by  N.YEAR 

40  reserve  MOE.NAME(*),  MOE.WT(*),  O.MOE(*),  O.MOE.DF(*)  as  NUM.MOE 

41  reserve  MIN.MOE. REQ( * , *)  as  NUM.MOE  by  N.YEAR 
42 

43  call  CALCl.NODE      '•  Structuring  the  network. 
44 

45  create  every  NODE  and  ARC 

46  reserve  CRIT.FUNC(*)  as  N.NODE 

47  reserve  EP(*)  as  (N.NODE  +  1) 

48  reserve  LCE(*),  LCC(*),  LIM.MOE(*)  and  L(*)  as  N.ARC 
49 

50  call  CALC2.N0DE      •'  Computing  the  year  &  system  of  each  node 

51  '•  &  the  start  index  of  arcs  entering  it. 
52 

53  call  CALC. ALT. CODE   •'  Assigning  code  of  alternatives. 

54 

55  read  CMC. PI,  0.MC.P2,  O.OC,  O.OC.EF,  O.DC,  O.DC.EF,  O.SV,  O.SV.EF 
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56  for  I  =  1  to  N.ALT 

57  for  J  =  1  to  N.YEAR 

58  readAC(I,J),  MC.P1(I,J),  MC.P2(I,J),  OC(I,J),  OC.EF(I,J), 

59  DC(I,J),  DC.EF(I,J),  SV(I,J)  and  SV. EF ( I , J ) 

60  for  J  =  1  to  NUM.MOE 

61  read  MOE.NAME(J),  MOE.WT(J),  O.MOE(J)  and  O.MOE.DF(J) 

62  for  J  =  1  to  NUM.MOE 

63  for  K  =  1  to  N.YEAR 

64  read  MIN.MOE.REQ(J, K) 
6  5     for  I  =  1  to  N.ALT 

66  for  J  =  1  to  NUM.MOE 

67  for  K  =  1  to  N.YEAR 

68  read  MOE(I,J,K)  and  MOE. DF(I , J , K) 
69 

70     call  CALC.ARC        • •  Computes  tail  node,  LCC  and  LCE  for  each  arc. 
71 

72  print  3  line  with  CASE. NAME  thus 

73  Dynamic  End-Of-Lif e-Cycle  Planning  Model 
74 

75  Case  of  ************************************** 
76 

77  skip  1  line 

78  print  4  line  with  APPROACH  and  BUDGET  thus 

79  Approach  -  ******************      Budget  =  $*********.** 
80 

81  Note  :  Budget  is  used  only  for  the  Lagrangian  relaxation  approach  and  is 

82  ignored  in  the  other  approaches. 
83 

8  4     if  (  ECHO. ON  =  "On"  ) 
8  5        call  ECHO. INPUT 
86 

87  skip  1  line 

88  print  9  line  thus 

89  Candidate  Arcs  For  Solution  Path 
90 

91  Note  -  LCC  is  made  to  be  unfavorably  high  if  the  minimum  requirement 

92  level  for  any  component  MOE  of  LCE  has  not  been  met.  The  limiting  MOE 

93  is  given  as  the  highest  index  of  the  MOE(s)  that  failed  to  meet  the 

94  minimum  requirement  level (s)  during  the  life-span  covered  by  the  arc. 
95 

96  Arc   Alt.   Tail   Head  Life-Cycle      ■   Limiting 

97  No.   Code   Node   Node    Effectiveness    Cost       MOE 
98 

99     skip  1  line 

100  for  I  =  2  to  N.NODE 

101  for  K  =  EP(I)  to  (  EP(I+1)  -  1  ) 

102  print  1  line  with  K,  ALT.CODE(K),  TAIL(K) ,  I,  LCE(K),  LCC(K) 

103  and  LIM.MOE(K)  thus 

104  ***   ***    ***    ***     ******.****   *****,**      *** 

105  always 
106 

107     return 

108 

109  end 


102 


1  routine  CALCl.NODE 
2 

3  ' '  Structuring  the  network  by  determining  the  number  of 

4  ' '  nodes  and  arcs  in  each  system 
5 

6  BASE.YEAR(l)  =  CURR.YEAR  +  MIN.LIFE(l) 

7  MAX. PRED. NODES (1)  =  1 

8  NUM. NODES (1)  =  MAX.LIFE(l)  -  MIN.LIFE(l)  +  1 

9  NUM. ARCS (1)  =  NUM. NODES (1) 

10  N.NODE  =  NUM.NODES(l)  +  1 

11  N.ARC  =  NUM. ARCS (1) 

12  for  I  =  2  to  N.SYS 

13  do 

14  J  =  I  -  1 

15  BASE.YEAR(I)  =  BASE.YEAR(J)  +  MIN.LIFE(I) 

16  MAX.PRED.NODES(I)  =  NUM . NODES ( I-l ) 

17  NUM.NODES(I)  =  NUM.NODES(J)  +  MAX.LIFE(I)  -  MIN.LIFE(I) 

18  NUM. ARCS (I)  =  N.ALT  *  NUM. NODES (J)  * 

19  (  MAX.LIFE(I)  -  MIN.LIFE(I)  +  1  ) 
2  0         N.NODE  =  N.NODE  +  NUM. NODES (I) 

21  N.ARC  =  N.ARC  +  NUM. ARCS (I) 

22  loop 
23 

24  return 

25 

2  6  end 
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1  routine  CALC2.N0DE 
2 

3  ' •  Computes  for  each  node  information  such  as  the  year  and  system 

4  ' '  It  belongs  to,  as  well  as  the  lead  index  (  contained  in  entry 

5  ' '  point  array  EP  )  of  the  arcs  entering  the  node 
6 

7  TRANS. YEAR (1)  =  CURR.YEAR 

8  SYS.CODE(l)  =  0 

9  EP(1)  =  0 

10  K  =  2 

11  JJ  =  1 

12  for  J  =  1  to  NUM. NODES (1) 

13  do 

14  TRANS. YEAR(K)  =  BASE.YEAR(l)  +  J  -  1 

15  SYS.CODE(K)  =  1 

16  EP(K)  =  EP(K-l)  +  1 

17  K  =  K  +  1 

18  loop 

19  for  I  =  2  to  N.SYS 

20  do 

21  for  J  =  1  to  NUM.NODES(I) 

22  do 

2  3  TRANS. YEAR (K)  =  BASE.YEAR(I)  +  J  -  1 

24  SYS.CODE(K)  =  I 

25  EP(K)  =  EP(K-l)  +  JJ 

26  M  =  MAX.LIFE(I)  -  MIN.LIFE(I)  +  1 

27  JJ  =  N.ALT*min.f (J,MAX.PRED.NODES(I) ,M, (NUM.NODES(I)-J+l) ) 

28  K  =  K  +  1 

29  loop 

30  loop 

31  EP(N.NODE+l)  =  EP(N.NODE)  +  N.ALT 
32 

33     return 

34 

35  end 
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1  routine  CALC. ALT. CODE 

2 

3     ''  Assigns  alternative  codes  to  each  arc. 

4 

5  for  K  =  1  to  NUM. NODES (1) 

6  ALT. CODE (K)  =  0 

7  for  I  =  (  NUM. NODES (1)  +  2  )  to  N.NODE 

8  for  K  =  EP(I)  to  (  EP(I+1)  -  1  ) 
9 

10  ALT.CODE(K)  =  (  mod . f  (  K  -  EP(I)  ,  N.ALT  )  +  1  ) 

11 

12     return 

13 

14  end 


105 


1  routine  CALC.ARC 
2 

3  ' '  Determines  the  tail  node,  life-cycle  cost  &  effectiveness 

4  ' •  associated  with  each  arc. 
5 

6  ••  Initialising  LCE  &  Limiting  MOE  arrays 

7  for  K  =  1  to  N.ARC 

8  do 

9  LCE(K)  =  0 

10  LIM.MOE(K)  =  0 

11  loop 
12 

13  ' •  Computations  for  the  present  system  (  system  1  ) . 

14  '•  All  arcs  here  have  the  alternative  code  0. 

15  for  K  =  1  to  NUM.NODES(l) 

16  do 

17  TAIL(K)  =  1 

18  M  =  TRANS. YEAR (K+1)  -  CURR.YEAR 

19  LCC(K)  =  O.DC  *  exp.f(  O.DC.EF  *  M  )  -  O.SV  *  exp.f(  O.SV.EF  *  M  ) 

20  for  I  =  1  to  NUM. MOE 

21  do 

22  MOE. VALUE  =  O.MOE(I)  ••  LCC  is  made  very  high 

23  if  (  MOE. VALUE  <  MIN . MOE . REQ ( I , M+1 )  )   ••  if  any  MOE  require- 

24  LCC(K)  =  9999999  ••  ment  level  is  not  met 

25  LIM.MOE(K)  =  I 

26  always 

27  ••  Note  :  Utility  fn.  is  Value  /  Max. Value 

28  if  (  I  =  1  ) 

29  LCE(K)  =  LCE(K)  +  MOE.WT(I)  *  MOE. VALUE  /  180 

30  else 

31  LCE(K)  =  LCE(K)  +  MOE.WT(I)  *  MOE. VALUE  /  100 

32  always 

33  loop 
34 

35  for  J  =  1  to  M 

36  do 

37  AMC  =  O.MC.Pl  *  O.MC.P2  *  (  J  **  (  O.MC.P2  -  1  )  ) 

38  AOC  =  O.OC  *  exp.f(  O.OC.EF  *  J  ) 

39  LCC(K)  =  LCC(K)  +  (  AMC  +  AOC  )  *  CUM.DF(J) 

40  for  I  =  1  to  NUM. MOE 

41  do 

4  2  MOE. VALUE  =  O.MOE(I)  *  exp.f(  O.MOE.DF(I)  *  J  ) 

43  if  (  MOE. VALUE  <  MIN . MOE . REQ ( I , J+ 1 )  )      •'  Checking  If  MOE 

44  LCC(K)  =  9999999  '•  req.  were  met. 

45  LIM.MOE{K)  =  I 

46  always 

47  •'  Note  :  Utility  fn.  is  Value  /  Max. Value 

48  if  (  I  =  1  ) 

49  LCE(K)  =  LCE(K)  +  MOE.WT(I)  *  MOE. VALUE  /  180 

50  else 

51  LCE(K)  ■=  LCE(K)  +  MOE.WT(I)  *  MOE. VALUE  /  100 

52  always 

53  loop 

54  loop 

55  loop 
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56 

57  ••  Computations  for  later  systems,  represented  by 

58  '•  arcs  with  alternative  code  higher  than  0. 

59  for  S  =  2  to  N.SYS 

60  for  J  =  3  to  N.NODE 

61  with  SYS.CODE(J)  =  S 

62  for  I  =  2  to  N.NODE 

63  with  SYS.CODE(I)  =  S  -  1 

64  do 

65  M  =  TRANS. YEAR(J)  -  TRANS . YEAR ( I )    •'  Number  of  years  used. 

66  if  (  (  M  >=  MIN.LIFE(S)  )  and  (  M  <=  MAX.LIFE(S)  )  ) 

67  for  N  =  1  to  N.ALT 

68  do 

69  TAIL(K)  =1  •'  Start  node  of  the  arc  K. 

70  DCY  =  DC(N,I)  * 

71  exp.f(  DC.EF(N,I)  *  M  )   •'  Disposal  cost  & 

72  SVY  =  SV(N,I)  *  ••  salvage  value  in  the 

73  exp.f(  SV.EF(N,I)  *  M  )   '•  last  year. 
74 

75  YR  =  TRANS. YEAR(I)  -  CURR.YEAR     •'  No.  of  years 

76  '•  from  the  present. 

77  LCC(K)  =  AC(N,I)  *  CUM.DF(YR+1)  + 

78  (  DCY-SVY  )  *  CUM.DF(YR+1) 

79  for  E  =  1  to  NUM.MOE 

80  do 

81  MOE. VALUE  =  MOE (N, E, YR+1) 

82  '•  Note  :  Utility  fn.  is  Value  /  Max. Value 

83  if  (  E  =  1  ) 

84  LCE(K)  =  LCE(K)  +  MOE.WT(E)  *  MOE. VALUE  /  180 

85  else 

86  LCE(K)  =  LCE(K)  +  MOE.WT(E)  *  MOE. VALUE  /  100 

87  always 

8  8  if  (  MOE. VALUE  <  MIN .MOE. REQ(E, YR+1 )  ) 

89  LCC(K)  =  9999999  •  '  Checking  if  MOE 

90  LIM.MOE(K)  =  E  •'  req.  were  met. 

91  always 

92  loop 

93  for  Y  =  1  to  M 

94  do 

95  AMC  =  MC.P1(N,I)  *  MC.P2(N,I)  *  '•  Annual  maint. 

96  (  Y  **  (  MC.P2(N,IJ  -  1  )  )    •'  t  operating 

97  AOC  =  OC(N,I)  *  exp.f(  OC.EF(N,I)  *  Y  )   ''  costs. 

98  LCC(K)  =  LCC(K)  +  (  AMC  +  AOC  )  * 

99  CUM.DF(YR+M) 

100  for  E  =  1  to  NUM.MOE 

101  do 

102  MOE. VALUE  =  MOE (N , E, YR+ 1 )  * 

103  exp.f  (  MOE. DF(N,E, YR+1)  *  Y  ) 

104  ••  Note  :  Utility  fn.  is  Value  /  Max. Value 

105  if  (  E  =  1  ) 

106  LCE(K)  =  LCE(K)  +  MOE.WT(E)  *  MOE. VALUE  /  180 

107  else 

108  LCE(K)  ■=  LCE(K)  +  MOE.WT(E)  *  MOE. VALUE  /  100 

109  always 

110  if  (  MOE. VALUE  <  MIN . MOE . REQ ( E, YR+Y+1 )  ) 
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Ill 

LCC(K)  =  9999999 

112 

LIM.MOE(K)  =  E 

113 

always 

114 

loop 

115 

loop 

116 

compute  MAX.LCE  as  the 

117 

maximum  of  LCE(K) 

118 

if  (  LCC(K)  <  9999999  ) 

119 

compute  MIN.LCC  as  the 

120 

minimum  of  LCC(K) 

121 

always 

122 

K  =  K  +  1 

123 

loop 

124 

al 

ways 

125 

loop 

126 

127 

return 

128 

129 

end 

' •  Checking  if  MOE 
' '  req.  were  met. 


MAX.LCE  &  MIN.LCC  are 
needed  for  computing 
the  upper  bound  on  the 
Lagrangian  multiplier. 
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1  routine  ECHO. INPUT 

2 

3  ' '  Echo  of  important  inputs 

4 

5  skip  1  line 

6  print  1  line  thus 

7  Sys.Code   Min.Life   Max. Life   Base.Yr.   No. Nodes   No. Arcs 

8  For  I  =  1  to  N.SYS 

9  print  1  line  with  I,  MIN.LIFE(I),  MAX.LIFE(I),  BASE. YEAR (I) , 

10  NUM.NODES(I) ,  NUM.ARCS(I)  thus 

11  ***        ***        ***        ****       ****       **** 

12 

13  skip  1  line 

14  print  1  line  thus 

15  Alternative  Name         Alt. Code         Lead  Time 

16  For  I  =  1  to  N.ALT 

17  print  1  line  with  ALT.NAME(I),  I,  LEAD.TIME(I)  thus 

18  ********************         ***  *** 

19 

20  skip  1  line 

21  print  2  line  thus 

22  Measure  of  Effectiveness  Present  System 

23  Code       Name        Weightage  Level    Exp. RateOfChange 

24  for  I  =  1  to  NUM.MOE 

25  print  1  line  with  I,  MOE.NAME(I),  MOE.WT(I),  O.MOE(I),  O.MOE.DF(I) 

26  ***   ***************    *.***  ******.***       ****.*** 

27 

28  return 

29 

30  end 
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1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 


routine  LAGRANGIAN 

' '  Solution  by  Lagrangian  relaxation  method 

The  starting  lower  bound  for  the  Lagrangian  multiplier  lamda  is  0. 
The  upper  bound  is  determined  such  that  It  is  barely  sufficient  to 
ensure  that  the  corresponding  slope  will  be  positive.  Variables 
may  be  overwhelmed,   even  when  double  precision  is  used,  if  lamda 
or  the  cost  values  are  too  high. 

LAMDA  1  =  0 

LAMDA2  =  (  N.NODE  *  MAX . LCE  /  MIN.LCC  )  +  1 

skip  1  line 

print  3  line  thus 

Iterations  of  the  Lagrangian  Relaxation  Method 


Lamda 1 
Skip 


Lan>da3 
1  line 


Lamda2 


Slopel   Slope3   Slope2 


LI 


L3 


L2 


• '  Iterations  by  intercept  of  supporting  planes  method, 
•Iterate' 


LAMDA  =  LAMDAl 
call  BEST. PATH 
call  SUM.LCC.LCE 
SUM.LCC. 1  =  SUM.LCC 
SUM. LCE. 1  =  SUM. LCE 
SLOPEl  =  -1  *  (  SUM.LCC.l 
LI  =  SUM. LCE. 1  -  LAMDAl  * 


-  BUDGET  ) 

(  SUM.LCC.l  -  BUDGET  ) 


if  {  SLOPEl  >  0  ) 

skip  1  line 

print  2  line  thus 
The  budget  is  sufficient  for  cost  to  be  no  longer  a  constraint, 
ie.  the  solution  is  that  of  maximum  effectiveness  over  time. 

return 
always 

LAMDA  =  LAMDA2 
call  BEST. PATH 
call  SUM.LCC.LCE 
SUM.LCC. 2  =  SUM.LCC 
SUM. LCE. 2  =  SUM. LCE 
SLOPE2  =  -1  *  (  SUM.LCC. 2 


L2 


if 


SUM. LCE. 2  -  LAMDA2  * 


) 


-  BUDGET  ) 

(  SUM.LCC. 2  -  BUDGET  ) 


(  SLOPE2  <  0 
skip  1  line 
print  4  line  with  BUDGET  thus 
No  solution  is  possible  because  of  one  of  the  following  reasons! 

(1)  the  given  budget  of  ********.**  is  insufficient 

(2)  one  or  more  MOE  requirement  level  is  too  high 

(3)  cost  values  are  too  high  and  needs  to  be  scaled  down. 
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56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
90 
91 
92 
93 
94 
95 
96 
97 
98 
99 
100 
101 
102 
103 
104 
105 
106 


stop 
always 

If  (  (  SLOPEl  =  SLOPE2  )  and  (  LI  =  L2  )  ) 

skip  1  line 

print  1  line  thus 
The  cost  values  are  too  high  for  computation  &  should  be  scaled  down. 

stop 
always 

•'  Intercept  of  the  line  segments  at  lamdal  6  lamda2  yields  lamda3. 
LAMDA3  =  (  SUM.LCE.l  -  SUM.LCE.2  )  /  (  SUM.LCC.l  -  SUM.LCC.2  ) 

LAMDA  =  LAMDA3 

call  BEST. PATH 

call  SUM.LCC.LCE 

SLOPE3  =  -1  *  (  SUM.LCC  -  BUDGET  ) 

L3  =  SUM.LCE  -  LAMDA3  *  (  SUM.LCC  -  BUDGET  ) 

print  1  line  with  LAMDAl,  LAMDA3,  LAMDA 2 ,  SLOPEl,  SLOPE3 ,  SLOPE2, 
LI,  L3,  L2  thus 

A**.***  «***,***  A***,***    ****,**  «***,**  ****.**    *•**,*  •***,*   *«**.* 

•'  The  minima  of  the  piece-wise  linear  Lagrangian  function  L 

•'  is  reached,  ie.  the  solution  is  found,  when  the  intercept 

' '  lamda3  is  equal  to  any  of  the  two  bounds 

if  (  LAMDA3  =  LAMDA2  )    •'  Minima  is  found  &  slope  there  is  positive 

return  *'  ie.  minima  Itself  Is  the  optimal  solution 

always 
if  (  LAMDA3  =  LAMDAl  )    '•  Minima  is  found  but  slope  is  negative. 

LAMDA  =  LAMDA2  "'  Lanida2  lies  on  the  solution  line  segment, 

call  BEST. PATH  '•  which  is  the  one  with  positive  slope  next 

call  SUM.LCC.LCE       '•  to  the  minima. 

SLOPE3    =    -1    *     (    SUM.LCC    -    BUDGET    ) 

L3  =  SUM.LCE  -  LAMDA 3  *  (  SUM.LCC  -  BUDGET  ) 

print  1  line  with  LAMDAl,  LAMDA3,  LAM DA 2 ,  SLOPEl,  SL0PE3,  SL0FE2, 
LI,  L3,  L2  thus 
***.***  -****,**«  ****,***   ****,«*  *«**,**  «««*.**   ****.*  ****,*  «***,* 

return 
always 

"'  Lamda3  replaces  the  redundant  bound  for  the  next  iteration, 
if  (  SL0PE3  >=  0  ) 

let  LAMDA2  =  LAMDA3 
else 

let  LAMDAl  =  LAMDA3 
always 


go  to  'Iterate' 


end 
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1  routine  BEST. PATH 
2 

3  • •  Forward  dynamic  programming  with  criteria  function  dependent  on 

4  •'  the  given  approach.  Provisions  are  made  to  ensure  that  the  solu- 

5  • '  tion  path  has  sufficient  lead  time  for  each  alternative  chosen. 
6 

7  CRIT.FUNC(l)  =  0 

8  for  I  =  2  to  N.NODE  ''  Initialising  crlt.  func.  with  a  very 

9  CRIT.FUNC(I)  =  -99999999    ''  low  value  as  0  is  too  high  a  default 

10  •'  value  to  prevent  adverse  selections. 

11  for  I  =  2  to  N.NODE 

12  do 

13  MAX.CF  =  -999999999 

14  for  J  =  EP(I)  to  (EP(I+1)-1) 

15  do 

16  • '  Lagrangian  relaxation  approach  to  obtain  the 

17  '•  max.  eff.  for  the  given  budget. 

18  if  (  APPROACH  =  "Lagrangian"  ) 

19  if  (  LCC{J)  <  9999999  ) 

20  L(J)  =  (  LCE(J)  -  LAMDA  *  LCC(J)  ) 

21  else 

22  L(J)  =  -9999999 

23  always 

24  always 
25 

26  ' '  The  approach  of  maximising  the  cost-eff . 

27  ' •  ratio  while  meeting  given  requirements  for 

28  ' '  effectiveness.  The  budget  is  ignored. 

29  if  (  APPROACH  =  "Cost-Eff . Ratio"  ) 

30  if  (  LCC(J)  <  9999999  ) 

31  L(J)  =  LCE(J)  /  (  LCC(J)  -  MIN.LCC  +  1  ) 

32  else 

33  L(J)  =  -9999999 

34  always 

35  always 
36 

37  • •  The  approach  of  satisfying  effectiveness 

38  '•  requirements  at  minimum  cost. 

39  '•  Yields  the  lower  bound  on  budget. 

40  if  (  APPROACH  =  "Satisficing"  ) 

41  L(J)  =  -1  *  LCC(J) 

42  always 
43 

44  • *  The  approach  of  maximising  effectiveness  without 

45  •'  concern  for  budget.  Yields  the  upper  bound  on  budget. 

46  if  (  APPROACH  =  "Max. Eff."  ) 

47  if  (  LCC(J)  <  9999999  ) 

48  L(J)  =  LCE(J) 

49  else 

50  L(J)  =  -9999999 

51  always 

52  always 
53 

54  ' '  Preventing  the  selection  of  an  arc  that  cannot 

55  • •  be  used  because  of  Insufficient  lead  time. 
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56  if  (  ALT.CODE(J)  >  0  ) 

57  if  (  BEST.PRED.NODE(TAIL(J) )  =  0  ) 

58  L(J)  =  -9999999 

59  else 

60  K  =  TRANS. YEAR(TAIL(J) )  - 

61  TRANS. YEAR(BEST.PRED.NODE(TAIL(J) ) ) 

62  If  (  K  <  LEAD,TIME(ALT.CODE(J) )  ) 

63  L(J)  =  -9999999 

64  always 

65  always 

66  always 
67 

68  If  (  CRIT.FUNC(TAIL(J) )  >  -9999999  ) 

69  If  (  CRIT.FUNC(TAIL(J) )  +  L(J)  )  >  MAX.CF 

70  MAX.CF  =  CRIT.FUNC(TAIL(J) )  +  L(J) 

71  CRIT.FUNC(I)  =  MAX.CF 

72  BEST.PRED.NODE(I)  •=    TAIL(J) 

73  BEST.PRED.ARC(I)  =  J 

74  always 

75  always 
76 

77  loop 

78  loop 
79 

80  • '  The  terminal  node  with  the  highest  value  in  the 

81  ••  criteria  function  yields  the  solution  path  that 

82  • '  is  optimal  for  the  given  approach. 

83  MAX.CF  =  -999999999 

84  MAX. MAC  =  -999999999 

85  MAX.AAC  =  -999999999 

86  for  I  =  1  to  N.NODE 

87  with  SYS.CODE(I)  =  SYS. CODE (N.NODE) 

88  do 

89  if  (  APPROACH  =  "Satisf icing"  ) 

90  N  =  TRANS.YEAR{I)  -  CURR.YEAR  +  1 

91  INT.DF  =  (  (  INT. RATE  /  100  )  +  1  )  **  (  N  -  1  ) 

92  AAC.CF  =  (  INT.DF  *  ( INT. RATE/100)  )  /  (  INT.DF  -  1 

93  If  (  (  CRIT.FUNC(I)  /  N  )  >  MAX. MAC  ) 

94  MAX. MAC  =  CRIT.FUNC(I)  /  N 

95  ARC. MAC  =  I 

96  always 

97  If  (  (  CRIT.FUNC(I)  *  AAC.CF  )  >  MAX.AAC  ) 

98  MAX.AAC  =  CRIT.FUNC(I)  *  AAC.CF 

99  ARG.AAC  =  I 

100  always 

101  always 

102  if  (  CRIT.FUNC(I)  >=  MAX.CF  ) 

103  MAX.CF  =  CRIT.FUNC(I) 

104  ARGMAX.CF  =  I 

105  always 

106  loop 
107 

108     return 

109 

110  end 
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1  routine  SUM.LCC.LCE 
2 

3  ' '  Totals  up  the  life-cycle  cost  &  effectiveness 

4  •'  for  the  solution  path. 
5 

6  I  =  ARGMAX.CF 

7  SUM.LCC  =  0 

8  SUM.LCE  =  0 

9  until  1=1 

10  do 

11  if  (  BEST.PRED.ARC(I)  =  0  ) 

12  print  2  lines  thus 

13  All  paths  violates  either  lead  time  or  MOE  level  requirements 

14  or  both.  No  feasible  solution  can  be  found. 

15  stop 

16  always 

17  SUM.LCC  =  SUM.LCC  +  LCC (BEST. PRED. ARC (I) ) 

18  SUM.LCE  =  SUM.LCE  +  LCE ( BEST . PRED. ARC ( I )  ) 

19  I  =  BEST. PRED. NODE (I) 

20  loop 
21 

22     return 

23 

24  end 
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1  routine  RESULTS , PRINT 

2 

3  ' '  Print-Out  of  Results 

4 

5  skip  1  line 

6  print  6  line  thus 

7  Network  of  Possible  Solutions  in  Nodal  Form 
8 

9  Note  -  The  criterion  function  of  a  node  will  be  very  low  if  all  paths 

10  leading  to  It  fall  to  satisfy  lead  time  req.  arid/or  minimum  MOE  levels. 
11 

12  No.   Sys.   Year   EP   Crit.Fn.   Best . Pred. Node   Best . Pred . Arc   Best. Alt. 

13  skip  1  line 

14  print  1  line  with  SYS.CODE(l),  TRANS . YEAR ( 1) ,  EP(1)  and 

15  CRIT.FUNC(l)  thus 

16  1   ***    **     *    ****.*** 

17  for  I  =  2  to  N.NODE 

18  print  1  line  with  I,  SYS.CODE(I),  TRANS . YEAR (I) ,  EP(I), 

19  CRIT.FUNC(I) ,  BEST. PRED. NODE(I) ,  BEST . PRED . ARC ( I )  and 

20  ALT. CODE(max.f(l, BEST. PRED. ARC(I) ) )  thus 

21  ***   ***    **     *    ****.***  ***               ***           *** 

22 

23  if  (  CRIT.FUNC(ARGMAX.CF)  <=  -99999  ) 

24  print  1  line  thus 

25  Warning  :  Lead-time  requirements  cannot  be  met. 

26  always 
27 

28  if  (  SUM.LCC  >=  9999999  ) 

29  print  1  line  thus 

30  Warning  :  Minimum  MOE  level  requirements  cannot  be  met. 

31  always 
32 

33  if  (  APPROACH  =  "Satisf icing"  )   ••  There  are  two  possible  solutions 

34  skip  1  line  ' '  for  the  satisflclng  approach. 

35  print  2  line  thus 

36  Minimum  Average  Annual  Cost  Solution  Path 

37  {  Time  value  of  money  is  ignored,  ie.  interest  rate  is  assumed  to  be  0  ) 

38  skip  1  line 

39  print  2  line  thus 

40  Transition         Alternative   Life-Cycle   Life-Cycle     Criterion 

41  Node     Year     System    Employed       Cost     Effectiveness   Func. Value 

42  skip  1  line 
4  3  I  =  ARC. MAC 
44  until  1=1 
4  5  do 

46  print  1  line  with  I,  TRANS . YEAR (I ) ,  SYS.CODE(I), 

47  ALT. CODE(BEST. PRED. ARC(I) ) ,  LCC ( BEST . PRED. ARC ( I ) ) , 

48  LCE(BEST.PRED.ARC(I) ) ,  CRIT.FUNC(I)  thus 

49  ***  *•        *  ***        *****  ^  ***  *****.***     *****.♦** 

50  I  =  BEST. PRED. NODE (I) 

51  loop 

52  skip  1  line 

53  print  2  line  with  INT. RATE  thus 

54  Minimum  Annualised  Cost  Solution  Path 

55  (  Monetary  Interest  assumed  constant  at  ***.***  %  per  annum  ) 
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56  skip  1  line 

57  print  2  line  thus 

58  Transition  Alternative   Life-Cycle   Life-Cycle     Criterion 

59  Node     Year     System    Employed  Cost     Effectiveness   Func. Value 

60  skip  1  line 

61  I  =  ARG.AAC 

62  until  1=1 

63  do 

64  print  1  line  with  I,  TRANS . yEAR(I) ,  SYS.CODE(I), 

65  ALT.CODE(BEST.PRED.ARC(I) ) ,  LCC ( BEST . PRED. ARC ( I ) ) , 

66  LCE(BEST.PRED.ARC(I) ) ,  CRIT.FUNC(I)  thus 

67  ***  **        *          ***                  *****.***       *****.***     *****.*** 

68  I  =  BEST. PRED. NODE (I) 

69  loop 

70  return 

71  always 
72 

73  skip  1  line 

74  print  4  line  thus 

75  Solution  Path 
76 

77  Transition        Alternative   Life-Cycle   Life-Cycle     Criterion 

78  Node     Year     System    Employed       Cost     Effectiveness   Func. Value 
79 

80  skip  1  line 

81  I  =  ARGMAX.CF 

82  until  1=1 

83  do 

84  print  1  line  with  I,  TRANS . yEAR(I) ,  SYS.CODE(I), 

85  ALT.CODE(BEST.PRED.ARC(I) ) ,  LCC ( BEST. PRED. ARC ( I )) , 

86  LCE(BEST.PRED.ARC(I) ) ,  CRIT.FUNC(I)  thus 

87  ***         -*  *  *  ***  «****,***         ***AA,*A*       *****.*** 

8  8        I  =  BEST. PRED. NODE (I) 

8  9  loop 

90 

91  skip  1  line 

92  if  (  (  CRIT,FUHC(ARGMAX.CF)  >  -99999  )  and  (  SUM. LCC  <  9999999  )  ) 

93  print  2  line  with  SUM. LCC  and  SUM.LCE  thus 

94  The  total  life-cycle  cost  and  effectiveness  for  the  optimal  solution 

95  are  *******,**  and  ******.***  respectively. 

96  always 
97 

98  skip  1  line 

99  if  (  APPROACH  =  "Lagrangian"  ) 
300        print  1  line  with  LAMDA3  thus 

101  Optimal  value  of  the  Lagrangian  multiplier  is  *****.***** 

102  skip  1  line 

103  always 
104 

105  return 

106 

107  end 
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APPENDIX  B:  SAMPLE  INPUT  FILES 


stationary 

Env 

Ironment 

3  2  3 

1989 

0 

Off 

0  5 

5  10 

5  10 

Alt. A 

0 

Alt.B 

0 

1111111111111 
1111111111111 


Lag 

rangl 

an 
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30 
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0 

0 

20  - 

.2 
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0 
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.05 

0 
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0 
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0 
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2 
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15 

.05 

0 

0 
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2 
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15 

.05 

0 

0 
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2 
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15 

.05 

0 

0 
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2 
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15 

.05 

0 
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2 
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15 
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0 
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2 
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15 
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15 
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2 
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15 
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0 

0 
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2 
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300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

300  15  1.3  4  .05  0  0  90  -.2 

Range(niti)  0.4  104.6  -.03 
Accuracy(%)  0.3  87.2  -.02 
Availability(%)  0.3   93.5    -.03 

95  95  95  95  95  95  95  95  95  95  95  95  95 
95  95  95  95  95  95  95  95  95  95  95  95  95 

80  80  80  80  80  80  80  80  80  80  80  80  80 
80  80  80  80  80  80  80  80  80  80  80  80  80 

70  70  70  70  70  70  70  70  70  70  70  70  70 
70  70  70  70  70  70  70  70  70  70  70  70  70 

120  -.03  120  -.03  120  -.03  120  -.03  120  -.03  120  -.03 

120  -.03  120  -.03  120  -.03  120  -.03  120  -.03  120  -.03 

120  -.03  120  -.03  120  -.03  120  -.03  120  -.03  120  -.03 

120  -.03  120  -.03  120  -.03  120  -.03  120  -.03  120  -.03 

120  -.03  120  -.03 

95  -.02  95  -.02  95  -.02  95  -.02  95  -.02  95  -.02 

95  -.02  95  -.02  95  -.02  95  -.02  95  -.02  95  -.02 

95  -.02  95  -.02  95  -.02  95  -.02  95  -.02  95  -.02 

95  -.02  95  -.02  95  -.02  95  -.02  95  -.02  95  -.02 

95  -.02  95  -.02 

95  -.03  95  -.03  95  -.03  95  -.03  95  -.03  95  -.03 

95  -.03  95  -.03  95  -.03  95  -.03  95  -.03  95  -.03 

95  -.03  95  -.03  95  -.03  95  -.03  95  -.03  95  -.03 

95  -.03  95  -.03  95  -.03  95  -.03  95  -.03  95  -.03 

95  -.03  95  -.03 

140  -.03  140  -.03  140  -.03  140  -.03  140  -.03  140  -.03 

140  -.03  140  -.03  140  -.03  140  -.03  140  -.03  140  -.03 

140  -.03  140  -.03  140  -.03  140  -.03  140  -.03  140  -.03 

140  -.03  140  -.03  140  -.03  140  -.03  140  -.03  140  -.03 

140  -.03  140  -.03 

95  -.02  95  -.02  95  -.02  95  -.02  95  -.02  95  -.02 

95  -.02  95  -.02  95  -.02  95  -.02  95  -.02  95  -.02 

95  -.02  95  -.02  95  -.02  95  -.02  95  -.02  95  -.02 

95  -.02  95  -.02  95  -.02  95  -.02  95  -.02  95  -.02 

95  -.02  95  -.02 

95  -.03  95  -.03  95  -.03  95  -.03  95  -.03  95  -.03 

95  -.03  95  -.03  95  -.03  95  -.03  95  -.03  95  -.03 

95  -.03  95  -.03  95  -.03  95  -.03  95  -.03  95  -.03 

95  -.03  95  -.03  95  -.03  95  -.03  95  -.03  95  -.03 

95  -.03  95  -.03 


1J8 


Inflation. &. Interest .Rates 

3  2  3  1989  5 

Off 

0  5 

5  10 

5  10 

Alt. A   0 

Alt.B   0 

.95  .95  .95  .95  .95  .95 

.95  .95  .95  .95  .95  .95 

.95  .95  .95  .95  .95  .95 

.95  .95  .95  .95  .95  .95 

,95  .95 


Satisf icing 

900 

30.0 

,  3 

7. 

.0 

.05 

0 

0 

20 

-.2 

100. 

,0 

31.2 

,  3 

7, 

,3 

.05 

0 

0 

50. 

.0 

-,2 

104  . 

.  0 

32.4 

.  3 

7. 

.6 

.05 

0 

0 

52. 

.0 

-,2 

108, 

.  2 

33.7 

,  3 

7. 

.9 

.05 

0 

0 

54  . 

.  1 

-,  2 

112. 

.  5 

35.1 

.  3 

8. 

.2 

.05 

0 

0 

56. 

,  2 

-,2 

117. 

,0 

36.5 

.3 

8. 

.  5 

,05 

0 

0 

58. 

.  5 

-,2 

121. 

.  7 

38.0 

.3 

8. 

.9 

,05 

0 

0 

60. 

.8 

-.2 

126. 

.5 

39.5 

.3 

9. 

.  2 

,05 

0 

0 

63  , 

,2 

-,2 

131. 

.6 

41.1 

.3 

9, 

.6 

,05 

0 

0 

65. 

.8 

-.2 

136. 

,9 

42.7 

,3 

10. 

.0 

.05 

0 

0 

68. 

.4 

-.2 

142. 

,3 

44  .4 

.3 

10. 

.4 

.05 

0 

0 

71. 

.2 

-.2 

148. 

.0 

46.2 

.3 

10. 

.8 

.05 

0 

0 

74. 

.0 

-.2 

153, 

.9 

48.0 

,3 

11. 

.2 

.05 

0 

0 

77, 

.0 

-.2 

160. 

.1 

50.0 

.3 

11. 

.7 

.05 

0 

0 

80, 

.1 

-.2 

166. 

.5 

52.0 

,3 

12. 

.  1 

.05 

0 

0 

83, 

,3 

-.2 

173, 

.2 

54.0 

,3 

12. 

.6 

.05 

0 

0 

86, 

.6 

-.2 

180, 

.1 

56.2 

,3 

13, 

.1 

.05 

0 

0 

90, 

.0 

-.2 

187, 

.3 

58.4 

,3 

13, 

.6 

.05 

0 

0 

93, 

.6 

-.2 

194, 

.8 

60,8 

,  3 

14, 

.2 

.05 

0 

0 

97, 

.4 

-.2 

202. 

.6 

63.2 

,  3 

14, 

.7 

.05 

0 

0 

101, 

.  3 

-.2 

210, 

.7 

65.7 

,3 

15, 

.3 

.05 

0 

0 

105, 

.  3 

-.2 

219, 

.1 

68.4 

,3 

16, 

.0 

,05 

0 

0 

109, 

.6 

-.2 

227, 

.9 

71.1 

,3 

16, 

.6 

,05 

0 

0 

113, 

.9 

-.2 

237  , 

.0 

73.9 

,3 

17, 

.2 

.05 

0 

0 

118, 

.5 

-.2 

246, 

.5 

76.9 

,3 

17, 

.9 

,05 

0 

0 

123, 

.2 

-,2 

256, 

.3 

80.0 

,3 

18, 

.7 

,05 

0 

0 

128, 

.2 

-,2 

266, 

.5 

83.2 

.3 

19, 

.4 

,05 

0 

0 

133, 

.3 

-,2 

300, 

.0 

15,6 

,  3 

4  , 

.2 

,05 

0 

0 

90, 

.0 

-,2 

312, 

.0 

16,2 

.  3 

4  . 

.  3 

,05 

0 

0 

93, 

.6 

-.2 

324  . 

.5 

16,9 

,3 

4  . 

.5 

,05 

0 

0 

97. 

.3 

-,2 

337, 

.5 

17.5 

.3 

4. 

.7 

,05 

0 

0 

101. 

.2 

-.  2 

351. 

.0 

18.2 

.3 

4  . 

.9 

,05 

0 

0 

105. 

.3 

-,2 

365. 

.0 

19,0 

.3 

5. 

.1 

.05 

0 

0 

109. 

.5 

-.2 

379. 

,6 

19,7 

,3 

5, 

,3 

,05 

0 

0 

113. 

.9 

-,2 

394. 

.8 

20.5 

.3 

5, 

.5 

,05 

0 

0 

118. 

,4 

-,2 

410. 

,6 

21.3 

3 

5, 

.7 

,05 

0 

0 

123. 

,2 

-,2 

427, 

,0 

22.2 

,3 

5. 

,9 

,05 

0 

0 

128. 

.1 

-,2 

444, 

,1 

23.1 

,3 

6. 

,2 

.05 

0 

0 

133. 

,2 

-.2 

461, 

.8 

24.0 

3 

6. 

,4 

.05 

0 

0 

138. 

,6 

-.2 

480, 

.  3 

25.0 

,  3 

6. 

.7 

.05 

0 

0 

144  . 

1 

-.2 

499, 

.5 

26.0 

,3 

6. 

,9 

.05 

0 

0 

149, 

9 

-.2 

519. 

.5 

27.0 

3 

7. 

,2 

.05 

0 

0 

155. 

,9 

-.2 

119 


540, 

.3 

28, 

.1 

1.3 

7. 

.5 

.05 

0 

0 

162.1 

-.2 

561, 

.9 

29, 

.2 

1.3 

7. 

,8 

.05 

0 

0 

168.6 

-.2 

584, 

.4 

30, 

.4 

1.3 

8. 

.1 

.05 

0 

0 

175.3 

-.2 

607, 

.7 

31. 

.6 

1.3 

8. 

,4 

.05 

0 

0 

182.3 

-.2 

632, 

.1 

32, 

.9 

1.3 

8. 

.8 

.05 

0 

0 

189.6 

-.2 

657. 

,3 

34. 

,2 

1.3 

9. 

,1 

.05 

0 

0 

197.2 

-.2 

683. 

.6 

35, 

,5 

1.3 

9. 

.5 

.05 

0 

0 

205.1 

-.2 

711. 

.0 

37, 

,0 

1.3 

9. 

,9 

.05 

0 

0 

213.3 

-.2 

739. 

,4 

38. 

,4 

1.3 

10. 

,3 

.05 

0 

0 

221.8 

-.2 

769. 

.0 

40. 

,0 

1.3 

10. 

.7 

.05 

0 

0 

230.7 

-.2 

799. 

,8 

41. 

,6 

1.3 

11. 

,1 

.05 

0 

0 

239.9 

-.2 

Range(nin)  0.4  104.6  -.03 

Accuracy(%)  0.3   87.2  -.02 

Availability(%)  0.3   93.5  -.03 

95  95  95  95  95  95  95  95  95  95  95  95  95 

95  95  95  95  95  95  95  95  95  95  95  95  95 

80  80  80  80  80  80  80  80  80  80  80  80  80 

80  80  80  80  80  80  80  80  80  80  80  80  80 

70  70  70  70  70  70  70  70  70  70  70  70  70 

70  70  70  70  70  70  70  70  70  70  70  70  70 


120 
120 
120 
120 
120 


.03  120 

.03  120 

.03  120 

.03  120 

.03  120 


03  120 

03  120 

03  120 

03  120 
03 


03  120 

03  120 

03  120 

03  120 


.03  120 

.03  120 

.03  120 

.03  120 


.03  120 

.03  120 

.03  120 

.03  120 


03 
03 
03 
03 


95 
95 
95 
95 
95 


02  95 

02  95 

02  95 

02  95 

02  95 


.02  95 

.02  95 

.02  95 

.02  95 
.02 


02  95 

02  95 

02  95 

02  95 


02  95 

02  95 

02  95 

02  95 


02  95 

02  95 

02  95 

02  95 


02 

02 
02 
02 


95 
95 
95 
95 
95 


03  95 

03  95 

03  95 

03  95 

03  95 


.03  95 

.03  95 

.03  95 

.03  95 
.03 


03  95 

03  95 

03  95 

03  95 


03  95 

03  95 

03  95 

03  95 


03  95 

03  95 

03  95 

03  95 


.03 
.03 
.03 
.03 


140 
140 
140 
140 
140 


03  140 
03  140 
03  140 
03  140 
03  140 


03  140 

03  140 

03  140 

03  140 
03 


03  140 

03  140 

03  140 

03  140 


03  140 
03  140 
03  140 
03  140 


.03  140 

.03  140 

.03  140 

.03  140 


03 
03 
03 
03 


95 
95 
95 
95 
95 


02  95 

02  95 

02  95 

02  95 

02  95 


02  95 

02  95 

02  95 

02  95 
02 


02  95 

02  95 

02  95 

02  95 


02  95 

02  95 

02  95 

02  95 


.02  95 

.02  95 

.02  95 

.02  95 


02 
02 
02 
02 


95 
95 
95 
95 
95 


03  95 

03  95 

03  95 

03  95 

03  95 


03  95 

03  95 

03  95 

03  95 
03 


03  95 

03  95 

03  95 

03  95 


03  95 

03  95 

03  95 

03  95 


03  95 

03  95 

03  95 

03  95 


03 
03 
03 
03 


120 


Technological .Change 

3  2  3  1989  5 

Off 

0  5 

5  10 

5  10 

Alt. A   0 

Alt.B   0 

.95  .95  .95  .95  .95  .95 

.95  .95  .95  .95  .95  .95 

.95  .95  .95  .95  .95  .95 

.95  .95  .95  .95  .95  .95 
.95  .95 

Satisf icing 
650 

30.0  1.3  7.0  .05  0  0  20    -.2 

100.0  30.3  1.3  7.3  .05  0  0  50.0  -.2 

102.0  30.6  1.3  7.6  .05  0  0  51.0  -.2 
104.4  30.9  1.3  7.9  .05  0  0  52.2  -.2 

106.1  31.2  1.3  8.2  .05  0  0  53.0  -.2 

108.2  31.5  1.3  8.5  .05  0  0  54.1  -.2 

110.4  31.8  1.3  8.9  .05  0  0  55.2  -.2 
112.6  32.2  1.3  9.2  .05  0  0  56.3  -.2 
114.9  32.5  1.3  9.6  .05  0  0  57.4  -.2 

117.2  32.8  1.3  10.0  .05  0  0  58.6  -.2 

119.5  33.1  1.3  10.4  .05  0  0  59.7  -.2 
121.9  33.5  1.3  10.8  .05  0  0  60.9  -.2 

124.3  33.8  1.3  11.2  .05  0  0  62.1  -.2 

126.8  34.1  1.3  11.7  .05  0  0  63.4  -.2 

129.4  34.5  1.3  12.1  .05  0  0  64.7  -.2 

131.9  34.8  1.3  12.6  .05  0  0  65.9  -.2 

134.6  35.2  1.3  13.1  .05  0  0  67.3  -.2 
137.3  35.5  1.3  13.6  .05  0  0  68.6  -.2 

140.0  35.9  1.3  14.2  .05  0  0  70.0  -.2 
142.8  36.2  1.3  14.7  .05  0  0  71.4  -.2 

145.7  36.6  1.3  15.3  .05  0  0  72.8  -.2 
148.6  37.0  1.3  16.0  .05  0  0  74.3  -.2 
151.6  37.3  1.3  16.6  .05  0  0  75.8  -.2 

154.6  37.7  1.3  17.2  ,05  0  0  77.3  -.2 

157.7  38.1  1.3  17.9  .05  0  0  78.8  -.2 

160.8  38.5  1.3  18.7  .05  0  0  80.4  -.2 

164.1  38.9  1.3  19.4  .05  0  0  82.0  -.2 

300.0  15.6  1.3  4.2  .05  0  0  90.0  -.2 

303.0  16.2  1.3  4.3  .05  0  0  90.9  -.2 

306.0  16.9  1.3  4.5  .05  0  0  91.8  -.2 

309.1  17.5  1.3  4.7  .05  0  0  92.7  -.2 

312.2  18.2  1.3  4.9  .05  0  0  93.7  -.2 

315.3  19.0  1.3  5.1  .05  0  0  94.6  -.2 

318.5  19.7  1.3  5.3  .05  0  0  95.5  -.2 

321.6  20.5  1.3  5.5  .05  0  0  96.5  -.2 

324.9  21.3  1.3  5.7  .05  0  0  97.5  -.2 
328.1  22.2  1.3  5.9  .05  0  0  98.4  -.2 

331.4  23.1  1.3  6.2  .05  0  0  99.4  -.2 

334.7  24.0  1.3  6.4  .05  0  0  100.4  -.2 
338.0  25.0  1.3  6.7  .05  0  0  101.4  -.2 
341.4  26.0  1.3  6.9  .05  0  0  102.4  -.2 

344.8  27.0  1.3  7.2  .05  0  0  103.5  -.2 
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348 

.3 

28. 

.1 

.3 

7, 

,5 

.05 

0 

0 

104. 

.5 

-.2 

351, 

.8 

29. 

.2 

,  3 

7, 

.8 

.05 

0 

0 

105, 

.5 

-.2 

355. 

.3 

30. 

.4 

.3 

8, 

.  1 

.05 

0 

0 

106, 

,6 

-.2 

358. 

.8 

31. 

,6 

.  3 

8, 

,4 

.05 

0 

0 

107, 

,7 

-.2 

362. 

.4 

32. 

.9 

.3 

8, 

.8 

.05 

0 

0 

108, 

,7 

-.2 

366. 

.  1 

34, 

.2 

,3 

9, 

,  1 

.05 

0 

0 

109, 

.8 

-.2 

369, 

.7 

35, 

.5 

,  3 

9, 

.5 

.05 

0 

0 

110, 

.9 

-.2 

373, 

,4 

37, 

.0 

.  3 

9, 

,9 

.05 

0 

0 

112, 

.0 

-.2 

377, 

.1 

38. 

.4 

.3 

10, 

,3 

.05 

0 

0 

113, 

.1 

-.2 

380. 

.9 

40, 

,0 

,3 

10, 

,7 

.05 

0 

0 

114. 

.3 

-.2 

384. 

,7 

41, 

.6 

.3 

11, 

,1 

.05 

0 

0 

115. 

.4 

-.2 

Range(nm)        0.4  104.6  -.03 

Accuracyd)      0.3   87.2  -.02 

Avallablllty(%)  0.3   93.5  -.03 

95  95  95  95  95  95  95  95  95  95  95  95  95 

95  95  95  95  95  95  95  95  95  95  95  95  95 

80  80  80  80  80  80  80  80  80  80  80  80  80 

80  80  80  80  80  80  80  80  80  80  80  80  80 

70  70  70  70  70  70  70  70  70  70  70  70  70 

70  70  70  70  70  70  70  70  70  70  70  70  70 


120 
150 
172 
180 
180 


.03  123 

.03  154 

.03  175 

.03  180 

.03  180 


03  128 
03  158 
03  178 
03  180 
03 


.03  134  -.03  140 

.03  162  -.03  166 

.03  179  -.03  180 

.03  180  -.03  180 


.03  145  -.03 

.03  169  -.03 

.03  180  -.03 

.03  180  -.03 


95.0 
95.6 
96.2 
96.7 
96.9 


.02  95.1 

.02  95.7 

.02  96.3 

.02  96.7 

.02  97.0 


02  95.2 

02  95.8 

02  96.4 

02  96.8 
02 


.02  95.3  -.02  95.4  -.02  95.5  -.02 

.02  95.9  -.02  96.0  -.02  96.1  -.02 

.02  96.4  -.02  96.5  -.02  96.6  -.02 

.02  96.8  -.02  96.9  -.02  96.9  -.02 


95.0 
95.6 
96.2 
96.7 
96.9 


.03  95.1 

.03  95.7 

.03  96.3 

.03  96.7 

.03  97.0 


03  95.2 
03  95.8 
03  96.4 
03  96.8 
03 


03  95.3  -.03  95.4  -.03  95.5  -.03 

03  95.9  -.03  96.0  -.03  96.1  -.03 

03  96.4  -.03  96.5  -.03  96.6  -.03 

03  96.8  -.03  96.9  -.03  96.9  -.03 


140 
158 
168 
173 
176 


.03  144 

.03  160 

.03  169 

.03  174 

.03  177 


03  147 
03  162 
03  170 
03  175 
03 


.03  150  -.03  153  -.03  155  -.03 

.03  164  -.03  165  -.03  167  -.03 

.03  171  -.03  172  -.03  173  -.03 

.03  175  -.03  176  -.03  176  -.03 


95.0 
95.6 
96.2 
96.7 
96.9 


02  95.1 
02  95.7 
02  96.3 
02  96.7 
02  97.0 


02  95.2 
02  95.8 
02  96.4 
02  96.8 
02 


02  95.3  -.02  95.4  -.02  95.5  -.02 

02  95.9  -.02  96.0  -.02  96.1  -.02 

02  96.4  -.02  96.5  -.02  96.6  -.02 

02  96.8  -.02  96.9  -.02  96.9  -.02 


95.0 
95.6 
96.2 
96.7 
96.9 


03  95.1 

03  95.7 

03  96.3 

03  96.7 

03  97.0 


-.03  95.2 

.03  95.8 

.03  96.4 

.03  96.8 
.03 


03  95.3  -.03  95.4  -.03  95.5  -.03 

03  95.9  -.03  96.0  -.03  96.1  -.03 

03  96.4  -.03  96.5  -.03  96.6  -.03 

03  96.8  -.03  96.9  -.03  96.9  -.03 
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Requirement .Escalation 
3  2  3  1989  5 


Off 

0  5 

5  10 

5  10 

Alt, 

.A 

0 

Alt, 

.B 

0 

.95 

.95  , 

,95  . 

.95 

i  .95 

.95 

.95 

.95  , 

.95  , 

,95 

i  .95 

.95 

.95 

.95  , 

.95  , 

.95 

i  .95 

.95 

.95 

.95  , 

.95  , 

.95 

i  .95 

.95 

.95 

.95 

Lagrangian 

700 

30, 

.0 

.  3 

7.0 

.05 

0 

0 

20 

-.2 

100. 

,0 

30, 

,3 

,3 

7.  3 

.05 

0 

0 

50, 

.0 

-.2 

102. 

,0 

30, 

,  6 

.3 

7.6 

.05 

0 

0 

51, 

.0 

-.2 

104  . 

,4 

30, 

.9 

.  3 

7  .  9 

.  05 

0 

0 

52  , 

.2 

-.2 

106. 

.  1 

31, 

.  2 

.  3 

8.2 

.05 

0 

0 

53, 

.0 

-.2 

108  . 

2 

31. 

.5 

.  3 

8.5 

.05 

0 

0 

54, 

.  1 

-.2 

110. 

.4 

31, 

.8 

.3 

8  .9 

.  05 

0 

0 

55. 

.2 

-.2 

112. 

.6 

32, 

.2 

.3 

9.2 

.05 

0 

0 

56. 

.3 

-.2 

114. 

.9 

32. 

.5 

.3 

9.6 

.05 

0 

0 

57. 

.4 

-.2 

117. 

2 

32, 

.8 

,3 

10.0 

.05 

0 

0 

58. 

.6 

-.2 

119. 

.5 

33, 

.1 

.3 

10.4 

.05 

0 

0 

59. 

.7 

-.2 

121. 

,9 

33, 

,5 

.3 

10.8 

.05 

0 

0 

60. 

.9 

-.2 

124. 

,3 

33, 

.8 

.3 

11.2 

.05 

0 

0 

62. 

.  1 

-.2 

126. 

,8 

34, 

.1 

,3 

11.7 

.05 

0 

0 

63. 

,4 

-.2 

129. 

.4 

34, 

,5 

,3 

12.1 

.05 

0 

0 

64  . 

,7 

-.2 

131, 

,9 

34. 

.8 

.3 

12.6 

.05 

0 

0 

65. 

.9 

-.2 

134. 

.6 

35. 

.2 

.  3 

13.1 

.05 

0 

0 

67. 

.3 

-.2 

137. 

.3 

35. 

.5 

.3 

13.6 

.05 

0 

0 

68. 

.6 

-.2 

140, 

.0 

35. 

.9 

.3 

14.2 

.05 

0 

0 

70, 

.0 

-.2 

142, 

.8 

36. 

.2 

.3 

14.7 

.05 

0 

0 

71, 

.4 

-.2 

145, 

,7 

36. 

.6 

.3 

15.3 

.05 

0 

0 

72, 

.8 

-.2 

148. 

,6 

37. 

.0 

.3 

16.0 

.05 

0 

0 

74, 

.3 

-.2 

151, 

,6 

37. 

.3 

.3 

16.6 

.05 

0 

0 

75, 

.8 

-.2 

154. 

,6 

37. 

.7 

.3 

17.2 

.05 

0 

0 

77, 

,3 

-.2 

157, 

,7 

38. 

.1 

.3 

17.9 

.05 

0 

0 

78, 

.8 

-.2 

160, 

,8 

38. 

.5 

.3 

18.7 

.05 

0 

0 

80, 

.4 

-.2 

164, 

,1 

38. 

,9 

.3 

19.4 

.05 

0 

0 

82. 

.0 

-.2 

300, 

,0 

15. 

.6 

.  3 

4.2 

.  05 

0 

0 

90. 

.0 

-.2 

303  , 

,0 

16, 

.  2 

,  3 

4  .3 

.05 

0 

0 

90. 

.9 

-.2 

306. 

,0 

16, 

.9 

.  3 

4  .5 

.05 

0 

0 

91, 

.8 

-.2 

309. 

1 

17. 

.5 

.  3 

4.7 

.05 

0 

0 

92, 

.7 

-.2 

312. 

2 

18. 

.2 

.3 

4.9 

.05 

0 

0 

93. 

,7 

-.2 

315. 

3 

19. 

.0 

,  3 

5.1 

.05 

0 

0 

94, 

,6 

-.2 

318. 

,5 

19, 

.7 

.3 

5.3 

.05 

0 

0 

95, 

,  5 

-.2 

321. 

,6 

20, 

.5 

.3 

5.5 

.05 

0 

0 

96, 

.5 

-.2 

324. 

,9 

21, 

.3 

,3 

5.7 

.05 

0 

0 

97. 

,5 

-.2 

328. 

1 

22. 

.2 

,3 

5.9 

.05 

0 

0 

98. 

-.2 

331. 

,4 

23. 

,  1 

.3 

6.2 

.05 

0 

0 

99. 

-.2 

334. 

7 

24. 

,0 

.3 

6.4 

.05 

0 

0 

100. 

-.2 

338. 

0 

25. 

,0 

.3 

6.7 

.05 

0 

0 

101. 

-.2 

341. 

4 

26. 

,0 

.3 

6.9 

.05 

0 

0 

102. 

-.2 

344  . 

8 

27. 

,0 

,  3 

7.2 

.05 

0 

0 

103. 

5 

-.2 

12,3 


348, 

.3 

28, 

.1 

.3 

7. 

,5 

.05 

0 

0 

104. 

.5 

-.2 

351, 

.8 

29, 

.2 

.3 

7. 

,8 

.05 

0 

0 

105. 

.5 

-.2 

355, 

.3 

30, 

.4 

,3 

8. 

.1 

.05 

0 

0 

106. 

.6 

-.2 

358. 

.8 

31, 

.6 

,3 

8. 

,4 

.05 

0 

0 

107, 

.7 

-.2 

362, 

.4 

32, 

.9 

,3 

8. 

.8 

.05 

0 

0 

108, 

.7 

-.2 

366, 

.1 

34, 

.2 

,3 

9. 

.1 

.05 

0 

0 

109, 

.8 

-.2 

369, 

.7 

35. 

.5 

,3 

9. 

,5 

.05 

0 

0 

110, 

.9 

-.2 

373. 

,4 

37, 

,0 

.3 

9. 

,9 

.05 

0 

0 

112, 

.0 

-.2 

377. 

.1 

38, 

.4 

,3 

10. 

.3 

.05 

0 

0 

113, 

.1 

-.2 

380. 

,9 

40. 

.0 

,3 

10. 

,7 

.05 

0 

0 

114, 

.3 

-.2 

384. 

,7 

41. 

,6 

,3 

11. 

.1 

.05 

0 

0 

115, 

.4 

-.2 

Range(nni)  0.4  104.6  -.03 
Accuracy(%)  0.3  87.2  -.02 
Availability(%)  0.3   93.5    -.03 

95   95   95  100  100  100  110  110  110  110  110  110  110 
110  110  110  110  110  110  110  110  110  110  110  110  110 

80  80  80  80  80  80  80  80  80  80  80  80  80 
80  80  80  80  80  80  80  80  80  80  80  80  80 

70  70  70  72  72  72  72  75  75  75  75  75  75 
75  75  75  75  75  75  75  75  75  75  75  75  75 


120 
150 
172 
180 
180 


.03  123 

.03  154 

.03  175 

.03  180 

.03  180 


03  128 

03  158 

03  178 

03  180 
03 


03  134 

03  162 

03  179 

03  180 


.03  140 

.03  166 

.03  180 

.03  180 


03  145 

03  169 

03  180 

03  180 


.03 
.03 
.03 
.03 


95.0 
95.  6 
96.2 
96.7 
96.9 


02  95.1 

02  95.7 

02  96.3 

02  96.7 

02  97.0 


02  95.2 

02  95.8 

02  96.4 

02  96.8 
02 


02  95.3 

02  95.9 

02  96.4 

02  96.8 


02  95.4 

02  96.0 

02  96.5 

02  96.9 


02  95.5 

02  96.1 

02  96.6 

02  96.9 


02 
02 
02 
02 


95.0 
95.6 
96.2 
96.7 
96.9 


03  95.1 

03  95.7 

03  96.3 

03  96,7 

03  97.0 


.03  95.2 

.03  95.8 

.03  96.4 

.03  96.8 
.03 


03  95.3 

03  95.9 

03  96.4 

03  96.8 


03  95.4 

03  96.0 

03  96.5 

03  96.9 


.03  95.5 

.03  96.1 

.03  96.6 

.03  96.9 


03 
03 
03 
03 


140 
158 
168 
173 
176 


.03  144 

.03  160 

.03  169 

.03  174 

.03  177 


03  147 

03  162 

03  170 

03  175 
03 


03  150 

03  164 

03  171 

03  175 


03  153 

03  165 

03  172 

03  176 


03  155 

03  167 

03  173 

03  176 


.03 
.03 
.03 
.03 


95.0 
95.6 
96.2 
96.7 
96.9 


.02  95, 1 

.02  95.7 

.02  96. 3 

.02  96.7 

.02  97.0 


02  95.2 
02  95.8 
02  96.4 
02  96.8 
02 


02  95.3 

02  95.9 

02  96.4 

02  96.8 


02  95.4 

02  96.0 

02  96.5 

02  96.9 


02  95.5 

02  96.1 

02  96.6 

02  96.9 


02 
02 
02 

02 


95.0 
95.6 
96.  2 

96.7 
96.9 


.03  95. 1 

.03  95.7 

.03  96. 3 

.03  96.7 

,03  97,0 


03  95.2 

03  95.8 

03  96.4 

03  96.8 
03 


03  95.3 

03  95.9 

03  96.4 

03  96.8 


.03  95.4 

.03  96.0 

.03  96.5 

.03  96.9 


03  95.5 

03  96.1 

03  96.6 

03  96.9 


03 
03 
03 
03 
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Upgrading . Scheduling 

3  2  3  1989  5 

Off 

0  5 

5  10 

5  10 

Alt. A   0 

Alt.B   0 

.95  .95  .95  .95  .95  .95 

.95  .95  ,95  .95  .95  .95 

.95  .95  .95  .95  .95  .95 

.95  .95  .95  .95  .95  .95 

.95  .95 

Cost-Eff. Ratio 


780. 

,08 

30. 

,0 

.3 

7. 

.0 

.05 

0 

0 

20 

-.2 

100. 

,0 

30. 

.3 

.3 

7. 

.3 

.05 

0 

0 

50. 

,0 

-.2 

82. 

.0 

30, 

.6 

,3 

7. 

,6 

.05 

0 

0 

51. 

.0 

-.2 

104  . 

.4 

30, 

.9 

.  3 

7. 

.9 

.05 

0 

0 

52. 

.2 

-.2 

106. 

.  1 

31, 

.2 

,  3 

8. 

.2 

.05 

0 

0 

53. 

.0 

-.2 

86, 

.2 

31, 

.5 

.3 

8. 

.5 

.05 

0 

0 

54. 

.  1 

-.2 

110, 

,4 

31, 

.  8 

.  3 

8  . 

.9 

.05 

0 

0 

55, 

.2 

-.2 

112, 

,  6 

32, 

.2 

.3 

9. 

.2 

.05 

0 

0 

56, 

.  3 

-.2 

90. 

.9 

32, 

.5 

.3 

9, 

.6 

.05 

0 

0 

57, 

.4 

-.2 

117, 

.2 

32, 

.8 

.  3 

10, 

.0 

.05 

0 

0 

58, 

.6 

-.2 

119, 

,  5 

33. 

.1 

.3 

10, 

,4 

.05 

0 

0 

59, 

.7 

-.2 

95, 

.9 

33, 

.5 

.3 

10, 

.8 

.05 

0 

0 

60, 

.9 

-.2 

124, 

.3 

33, 

.8 

.3 

11, 

.2 

.05 

0 

0 

62 

.1 

-.2 

126, 

.8 

34, 

.1 

.3 

11, 

.7 

.05 

0 

0 

63, 

.4 

-.2 

101, 

.4 

34, 

.5 

.3 

12, 

.  1 

.05 

0 

0 

64, 

.7 

-.2 

131, 

,9 

34, 

.8 

.3 

12, 

.6 

.05 

0 

0 

65, 

.9 

-.2 

134, 

.  6 

35, 

.2 

,3 

13, 

.1 

.05 

0 

0 

67, 

.3 

-.2 

107, 

.3 

35, 

.5 

.3 

13. 

,6 

.05 

0 

0 

68, 

,6 

-.2 

140, 

.0 

35, 

.9 

.3 

14. 

,2 

.05 

0 

0 

70, 

.0 

-.2 

142, 

.8 

36, 

.2 

.3 

14, 

.7 

.05 

0 

0 

71, 

.4 

-.2 

113, 

.7 

36. 

.  6 

.3 

15, 

.3 

.05 

0 

0 

72, 

.8 

-.2 

148, 

.6 

37, 

.0 

.3 

16, 

.0 

.05 

0 

0 

74, 

.  3 

-.2 

151, 

.6 

37, 

.3 

.3 

16, 

.6 

.05 

0 

0 

75, 

.8 

-.2 

120, 

.6 

37, 

.7 

,3 

17, 

,2 

.05 

0 

0 

77, 

.  3 

-.2 

157, 

.7 

38, 

.1 

.3 

17, 

.9 

.05 

0 

0 

78, 

.8 

-.2 

160, 

.8 

38, 

.  5 

.3 

18, 

,7 

.05 

0 

0 

80, 

.4 

-.2 

128, 

.  1 

38, 

,9 

,3 

19, 

.4 

.05 

0 

0 

82, 

.0 

-.2 

300. 

.  0 

15. 

.  6 

,3 

4  , 

.2 

.05 

0 

0 

90, 

,  0 

-.  2 

303. 

.  0 

16. 

.  2 

.  3 

4  . 

.  3 

.05 

0 

0 

90, 

.  9 

-.  2 

306. 

,0 

16. 

.9 

,  3 

4  . 

,5 

.05 

0 

0 

91. 

.8 

-.2 

309. 

1 

17. 

.5 

.3 

4  . 

.7 

.05 

0 

0 

92. 

.7 

-.2 

312. 

2 

18. 

,2 

.  3 

4  . 

,  9 

.05 

0 

0 

93  . 

.7 

-.2 

315. 

3 

19. 

.0 

.3 

5. 

,1 

.  05 

0 

0 

94. 

,6 

-.2 

318. 

.5 

19. 

,7 

,3 

5. 

.  3 

.05 

0 

0 

95. 

.5 

-.2 

321. 

,6 

20. 

.5 

,  3 

5. 

,5 

.05 

0 

0 

96. 

,5 

-.2 

324. 

.9 

21. 

,3 

,3 

5. 

.7 

.05 

0 

0 

97. 

.5 

-.2 

328. 

.1 

22. 

.2 

,3 

5. 

.9 

.05 

0 

0 

98. 

.4 

-.2 

331. 

.4 

23. 

,1 

,3 

6. 

,2 

.05 

0 

0 

99. 

,4 

-.2 

334  . 

,  7 

24. 

,0 

.  3 

6, 

,4 

.05 

0 

0 

100. 

,4 

-.2 

338. 

,0 

25. 

.0 

.  3 

6. 

,7 

.05 

0 

0 

101. 

,4 

-.2 

341. 

,4 

26. 

.0 

.3 

6. 

,9 

.05 

0 

0 

102. 

,  4 

-.  2 

344  . 

,8 

27. 

,0 

3 

7. 

,2 

.05 

0 

0 

103. 

.5 

-.2 

125 


348, 

.3 

28. 

.1 

1.3 

7.5 

.05 

0 

0 

104.5 

-.2 

351, 

.8 

29, 

.2 

1.3 

7.8 

.05 

0 

0 

105.5 

-.2 

355, 

.3 

30, 

.4 

1.3 

8.1 

.05 

0 

0 

106.6 

-.2 

358, 

.8 

31, 

.6 

1.3 

8.4 

.05 

0 

0 

107.7 

-.2 

362. 

.4 

32. 

.9 

1.3 

8.8 

.05 

0 

0 

108.7 

-.2 

366, 

,1 

34. 

.2 

1.3 

9.1 

.05 

0 

0 

109.8 

-.2 

369, 

.7 

35. 

,5 

1.3 

9.5 

.05 

0 

0 

110.9 

-.2 

373. 

.4 

37. 

.0 

1.3 

9.9 

.05 

0 

0 

112.0 

-.2 

377. 

,1 

38. 

.4 

1.3 

10.3 

.05 

0 

0 

113.1 

-.2 

380. 

.9 

40. 

,0 

1.3 

10.7 

.05 

0 

0 

114.3 

-.2 

384. 

,7 

41. 

,6 

1.3 

11.1 

.05 

0 

0 

115.4 

-.2 

Range(nin)  0.4  104.6  -.03 
Accuracy(%)  0.3  87.2  -.02 
Availability(%)  0.3   93.5    -.03 

95   95   95  100  100  100  110  110  110  110  110  110  110 
110  110  110  110  110  110  110  110  110  110  110  110  110 

80  80  80  80  80  80  80  80  80  80  80  80  80 
80  80  80  80  80  80  80  80  80  80  80  80  80 

70  70  70  72  72  72  72  75  75  75  75  75  75 
75  75  75  75  75  75  75  75  75  75  75  75  75 


120 
150 
172 
180 
180 


.03  123 

.03  154 

.03  175 

.03  180 

.03  180 


.03  128 

.03  158 

.03  178 

.03  180 
.03 


03  134 

03  162 

03  179 

03  180 


03  140 
03  166 
03  180 
03  180 


03  145 
03  169 
03  180 
03  180 


03 
03 
03 
03 


95.0 
95.6 
96.2 
96.7 
96.9 


02  95.1 

02  95.7 

02  96.3 

02  96.7 

02  97.0 


02  95.2 

02  95.8 

02  96.4 

02  96.8 
02 


02  95.3 

02  95.9 

02  96.4 

02  96.8 


02  95.4 

02  96.0 

02  96.5 

02  96.9 


.02  95.5 

.02  96. 1 

.02  96.6 

.02  96.9 


02 
02 
02 
02 


95.0 
95.6 
96.2 
96.7 
96.9 


03  95.1 

03  95.7 

03  96.3 

03  96.7 

03  97.0 


.03  95.2 

.03  95.8 

.03  96.4 

.03  96.8 
.03 


03  95.3 

03  95.9 

03  96.4 

03  96.8 


03  95.4 

03  96.0 

03  96. 5 

03  96.9 


.03  95.5 

.03  96. 1 

.03  96.6 

.03  96,9 


03 
03 
03 
03 


140 
158 
168 
173 
176 


03  144 

03  160 

03  169 

03  174 

03  177 


03  147 

03  162 

03  170 

03  175 
03 


03  150 

03  164 

03  171 

03  175 


03  153 
03  165 
03  172 
03  176 


03  155 

03  167 

03  173 

03  176 


.03 
.03 
.  03 
.03 


95.0 
95.6 
96,2 
96.7 
96.9 


02  95.  1 

02  95.7 

02  96.3 

02  96.7 

02  97.0 


02  95.2 
02  95.8 
02  96.4 
02  96.8 
02 


02  95. 3 

02  95.9 

02  96.4 

02  96.8 


.02  95,4 

.02  96.0 

.02  96. 5 

.02  96.9 


02  95.5 

02  96.1 

02  96.6 

02  96.9 


02 

02 
02 
02 


95.0 
95.6 
96.  2 
96.7 
96.9 


03  95.1 

03  95.7 

03  96.3 

03  96.7 

03  97.0 


03  95.2 

03  95.8 

03  96.4 

03  96.8 

03 


03  95.3 

03  95.9 

03  96.4 

03  96.8 


.03  95,4 

.03  96,0 

.03  96.5 

.03  96.9 


03  95.5 
03  96.1 
03  96.6 
03  96.9 


03 
03 
03 
03 
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Development . Lead . Time 
3    2    3    1989    5 


On 

0  5 

5  10 

5  10 

Alt, 

,  A 

1 

Alt. 

.B 

0 

.95 

.95  . 

.95  . 

.95 

i  .95 

.95 

.95 

.95  . 

.95  . 

,95 

i  .95 

.95 

.95 

.95  . 

.95  . 

,95 

i  .95 

.95 

.95 

.95  , 

.95  . 

.95  .95 

.95 

.95 

.95 

Lagrangian 

634. 

,52 

30. 

,0 

,3 

7.0 

.05 

0 

0 

20 

-.2 

100. 

,0 

30. 

.3 

.3 

7.  3 

.05 

0 

0 

50. 

.0 

-.2 

82. 

,0 

30. 

.6 

.3 

7.6 

.05 

0 

0 

51. 

.0 

-.2 

104. 

,4 

30. 

.9 

.  3 

7.9 

.05 

0 

0 

52. 

.2 

-.2 

106. 

.1 

31. 

.2 

.3 

8.2 

.05 

0 

0 

53. 

.0 

-.2 

86. 

,2 

31, 

.5 

.  3 

8.5 

.05 

0 

0 

54  , 

.1 

-.2 

110. 

.4 

31, 

.8 

.  3 

8.9 

.05 

0 

0 

55, 

.2 

-.2 

112. 

.6 

32 

.2 

.3 

9.2 

.05 

0 

0 

56, 

.3 

-.2 

90. 

,9 

32, 

.5 

.3 

9.6 

.05 

0 

0 

57, 

.4 

-.2 

117. 

,2 

32 

.8 

.3 

10.0 

.05 

0 

0 

58, 

.6 

-.2 

119. 

,5 

33, 

.  1 

.3 

10.4 

.05 

0 

0 

59, 

.7 

-.2 

95. 

.9 

33 

.5 

.  3 

10.8 

.05 

0 

0 

60, 

.9 

-.2 

124. 

.3 

33, 

.8 

.3 

11.2 

.05 

0 

0 

62, 

.1 

-.2 

126. 

.8 

34, 

.1 

.3 

11.7 

.05 

0 

0 

63, 

.4 

-.2 

101. 

.4 

34, 

.5 

.3 

12.1 

.05 

0 

0 

64, 

.7 

-.2 

131. 

.9 

34, 

.8 

.3 

12.6 

.05 

0 

0 

65, 

.9 

-.2 

134. 

.6 

35, 

.2 

,3 

13.1 

.05 

0 

0 

67, 

.3 

-.2 

107. 

.3 

35. 

.5 

,3 

13.6 

.05 

0 

0 

68, 

.6 

-.2 

140. 

,0 

35. 

.9 

.3 

14.2 

.  05 

0 

0 

70, 

.0 

-.2 

142. 

,8 

36. 

,2 

.3 

14.7 

.05 

0 

0 

71, 

.4 

-.2 

113. 

,7 

36. 

.6 

.3 

15.3 

.05 

0 

0 

72, 

.8 

-.2 

148. 

,6 

37, 

.0 

.3 

16.0 

.05 

0 

0 

74, 

.3 

-.2 

151. 

,6 

37. 

,3 

,3 

16.6 

.05 

0 

0 

75. 

.8 

-.2 

120. 

,6 

37. 

.7 

,3 

17.2 

.05 

0 

0 

77. 

.3 

-.2 

157. 

,7 

38. 

,1 

,3 

17.9 

.05 

0 

0 

78. 

.8 

-.2 

160. 

.8 

38. 

.5 

,3 

18.7 

.05 

0 

0 

80. 

.4 

-.2 

128. 

,1 

38. 

.9 

,3 

19.4 

.05 

0 

0 

82. 

.0 

-.2 

300. 

,0 

15. 

,6 

.  3 

4  .2 

.05 

0 

0 

90. 

.0 

-.2 

303. 

,0 

16. 

,2 

.3 

4  .3 

.05 

0 

0 

90. 

.9 

-.2 

306. 

,0 

16. 

,9 

.3 

4  .5 

.  05 

0 

0 

91. 

.8 

-.2 

309. 

,1 

17. 

,5 

.3 

4  .7 

.05 

0 

0 

92. 

.7 

-.2 

312. 

,2 

18. 

,2 

.3 

4  .9 

.05 

0 

0 

93. 

.7 

-.2 

315. 

3 

19. 

,0 

.3 

5.1 

.05 

0 

0 

94. 

.6 

-.2 

318. 

,5 

19. 

.7 

,  3 

5.3 

.05 

0 

0 

95. 

,5 

-.2 

321. 

6 

20. 

.5 

.3 

5.5 

.05 

0 

0 

96. 

,5 

-.2 

324. 

9 

21. 

,3 

,3 

5.7 

.05 

0 

0 

97. 

,5 

-.2 

328. 

1 

22. 

,2 

.3 

5.9 

.05 

0 

0 

98. 

,4 

-.2 

331. 

,4 

23. 

,1 

.3 

6,2 

.05 

0 

0 

99. 

,4 

-.2 

334. 

,7 

24. 

,0 

.3 

6.4 

.05 

0 

0 

100. 

.4 

-.2 

338. 

,0 

25. 

,0 

.3 

6.7 

.05 

0 

0 

101. 

,4 

-.2 

341. 

,4 

26. 

.0 

.3 

6.9 

.05 

0 

0 

102. 

.4 

-.2 

344. 

,8 

27. 

.0 

,  3 

7.2 

.05 

0 

0 

103. 

,5 

-.2 
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348. 

.3 

28. 

,1 

.3 

7.5 

.05 

0 

0 

104.5 

-.2 

351. 

.8 

29. 

.2 

.3 

7.8 

.05 

0 

0 

105.5 

-.2 

355, 

.  3 

30. 

.4 

,3 

8.1 

.05 

0 

0 

106.6 

-.2 

358. 

.8 

31. 

,6 

.3 

8.4 

.05 

0 

0 

107.7 

-.2 

362, 

.4 

32. 

.9 

.3 

8.8 

.05 

0 

0 

108.7 

-.2 

366. 

.1 

34. 

.2 

.3 

9.1 

.05 

0 

0 

109.8 

-.2 

369. 

.7 

35. 

.5 

.3 

9.5 

.05 

0 

0 

110.9 

-.2 

373. 

.4 

37. 

.0 

,3 

9.9 

.05 

0 

0 

112.0 

-.2 

377, 

.1 

38. 

.4 

,3 

10.3 

.05 

0 

0 

113.  1 

-.2 

380. 

.9 

40. 

,0 

.3 

10.7 

.05 

0 

0 

114.3 

-.2 

384. 

,7 

41. 

,6 

.3 

11.1 

.05 

0 

0 

115.4 

-.2 

Range(n)n)  0.4 
Accuracy(%)  0.3 
Availabillty(%)  0.3 


104.6  -.03 
87.2  -.02 
93.5    -.03 


95   95   95  100  100  100  110  110  110  110  110  110  110 

110  110  110  110  110  110  110  110  110  110  110  110  110 

80  80  80  80  80  80  80  80  80  80  80  80  80 

80  80  80  80  80  80  80  80  80  80  80  80  80 

70  70  70  72  72  72  72  75  75  75  75  75  75 

75  75  75  75  75  75  75  75  75  75  75  75  75 


120 
150 
172 
180 
180 


.03  123 

.03  154 

.03  175 

.03  180 

.03  180 


03  128 

03  158 

03  178 

03  180 


03  134  -.03  140 

03  162  -.03  166 

03  179  -.03  180 

03  180  -.03  180 


,03  145  -.03 

.03  169  -.03 

.03  180  -.03 

.03  180  -.03 


-.03 


95.0 
95.6 
96.2 
96.7 
96.9 


.02  95. 1 

.02  95.7 

.02  96. 3 

.02  96.7 

.02  97.0 


.02  95.2 

.02  95.8 

.02  96.4 

.02  96.8 
.02 


.02  95.3  -.02  95.4  -.02  95.5  -.02 

.02  95.9  -.02  96.0  -.02  96.1  -.02 

.02  96.4  -.02  96.5  -.02  96.6  -.02 

.02  96.8  -.02  96.9  -.02  96.9  -.02 


95.0 
95.6 
96.2 
96.7 
96.9 


.03  95.1 

.03  95.7 

.03  96.3 

.03  96.7 

.03  97.0 


.03  95.2 

.03  95.8 

.03  96.4 

.03  96.8 
.03 


.03  95.3  -.03  95.4  -.03  95.5  -.03 

.03  95.9  -.03  96.0  -.03  96.1  -.03 

.03  96.4  -.03  96.5  -.03  96.6  -.03 

.03  96.8  -.03  96.9  -.03  96.9  -.03 


140 
158 
168 
173 
176 


.03  144 

.03  160 

.03  169 

.03  174 

.03  177 


03  147 
03  162 
03  170 
03  175 
03 


.03  150  -.03  153  -.03  155  -.03 

.03  164  -.03  165  -.03  167  -.03 

.03  171  -.03  172  -.03  173  -.03 

.03  175  -.03  176  -.03  176  -.03 


95.0 
95.6 
96.2 
96.7 
96.9 


.02  95. 1 

.02  95.7 

.02  96.3 

.02  96.7 

.02  97.0 


02  95.2 

02  95.8 

02  96.4 

02  96.8 
02 


.02  95.3  -.02  95.4  -.02  95.5  -.02 

.02  95.9  -.02  96.0  -.02  96.1  -.02 

.02  96.4  -.02  96.5  -.02  96.6  -.02 

.02  96.8  -.02  96.9  -.02  96.9  -.02 


95.0 
95.6 
96.  2 
96.7 
96.9 


.03  95.1 

.03  95.7 

.03  96. 3 

.03  96.7 

.03  97.0 


.03  95.2 

.03  95.8 

.03  96. 4 

.03  96.8 

.03  . 


03  95.3  -.03  95.4  -.03  95.5  -.03 

03  95.9  -.03  96.0  -.03  96.1  -.03 

03  96.4  -.03  96.5  -.03  96.6  -.03 

03  96.8  -.03  96.9  -.03  96.9  -.03 
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0. 

.025 

.025 

.029 

.029 

.032 

.029 

.032 

27.04 

18.8 

11.4 

66.7 

27.04 

11.4 

11.4 

66.7 

27.04 

11.4 

11.4 

66.7 

27.04 

11.4 

11.4 

11.4 

APPENDIX  C:  SAMPLE  OUTPUT  FILES 

Dynamic  End-Of-Lif e-Cycle  Planning  Model 

Case  of  Technological  Change 

Approach  -  Lagrangian  Budget  =  $       601.74 

Note  :  Budget  is  used  only  for  the  Lagrangian  relaxation  approach  and  is 
ignored  in  the  other  approaches. 

Iterations  of  the  Lagrangian  Relaxation  Method 

Lamdal         Laiiida3         Lamda2           Slopel  Slope3      Slope2           LI  L3              L2 

2.080       -306.09  -32.79 

2.080          -32.79  -.01 

2.080               -.01  27.04 

.032               -.01  27.04 

Network  of  Possible  Solutions  in  Nodal  Form 

Note  -  The  criterion  function  of  a  node  will  be  very  low  if  all  paths 

leading  to  it  fail  to  satisfy  lead  time  reg.  and/or  minimum  MOE  levels. 

No.   Sys.   Year   EP   Crit.Fn.   Best . Pred. Node   Best . Pred . Arc   Best. Alt. 

1  0  1989     0       0. 

2  1  1989     1       1.411  1  10 

3  1  1990     2         .649  1  2             0 

4  1  1991     3        -.311  1  3             0 

5  1  1992     4      -1.374  1  4             0 

6  1  1993     5-1.000E+007  1  5             0 

7  1  1994     6-1.000E+007  1  6             0 

8  2  1994     7      -4.107  2  7 

9  2  1995     9      -4.476  3  11 

10  2  1996    13      -5.013  4  17 

11  2  1997    19      -5.659  5  25 

12  2  1998    27      -6.302  5  33 

13  2  1999    37      -6.876  5  43 

14  2  2000    49      -7.378  5  53 

15  2  2001    59-1.000E+007  4  60 

16  2  2002    67-1.000E+007  5  67 

17  2  2003    73-1.000E+008  0  0 

18  2  2004    77-1.000E+008  0  0 

19  3  1999    79      -7.772  8  79 

20  3  2000    81      -7.772  9  83 

21  3  2001    85      -7.962  10  89 

22  3  2002    91      -8.285  11  97 

23  3  2003    99      -8.598  12  107 

24  3  2004   109      -8.860  13  119 

25  3  2005   121      -9.075  14  131 

26  3  2006   133      -9.316  14  141 

27  3  2007   145      -9.510  14  151 

28  3  2008   157      -9.657  14  161 

29  3  2009   169      -9.755  14  171 

30  3  2010   181-1. OOOE+007  14  182 
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31 

3 

2011 

191-1. OOOE+008 

15 

32 

3 

2012 

199-1. OOOE+008 

16 

33 

3 

2013 

205-1. OOOE+008 

0 

34 

3 

2014 

209-1. OOOE+008 

0 

191  1 

199  1 

0  0 

0  0 

Solution  Path 

Transition        Alternative   Life-Cycle   Life-Cycle  Criterion 

Node     Year     System    Employed       Cost     Effectiveness  Func. Value 

1          273.366           5.029  -7.772 

1          321.339           4.702  -4.107 

0          -20.000            .775  1.411 

The  total  life-cycle  cost  and  effectiveness  for  the  optimal  solution 
are      574.70  and      10.506  respectively. 

Optimal  value  of  the  Lagrangian  multiplier  is       .03180 


19 

1999 

3 

8 

1994 

2 

2 

1989 

1 
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Dynamic  End-Of-Lif e-Cycle  Planning  Model 

Case  of  Upgrading. Scheduling 

Approach  -  Satisficing  Budget  =  $ 


780,08 


Note  :  Budget  is  used  only  for  the  Lagrangian  relaxation  approach  and  is 
ignored  in  the  other  approaches. 


Sys.Code   Min.Life  Max. Life  Base.Yr, 

10          5  1989 

2  5         10  1994 

3  5         10  1999 


Alternative  Name 

Alt. A 

Alt.B 


Alt. Code 
1 
2 


Measure  of  Effectiveness 
Code       Name        Weightage 

1  Range(nm)  .400 

2  Accuracy(%)  .300 

3  Availability(%)     .300 

Candidate  Arcs  For  Solution  Path 


No. Nodes  No. Arcs 
6  6 

11         72 
16        132 

Lead  Time 
0 
0 

Present  System 
Level    Exp. RateOf Change 
104.600         -.030 
87.200         -.020 
93.500  -.030 


Note  -  LCC  is  made  to  be  unfavorably  high  if  the  minimum  requirement 
level  for  any  component  MOE  of  LCE  has  not  been  met.  The  limiting  MOE 
is  given  as  the  highest  index  of  the  MOE(s)  that  failed  to  meet  the 
minimum  requirement  level (s)  during  the  life-span  covered  by  the  arc. 


Arc 

Alt. 

Tail 

Head 

Life-Cycle 

Limiting 

No. 

Code 

Node 

Node 

Effectiveness    Cost 

MOE 

1 

0 

1 

2 

.7745     -20.00 

0 

2 

0 

1 

3 

1.5287      27.67 

0 

3 

0 

1 

4 

2.2632      80.95 

0 

4 

0 

1 

5 

2. 97839999999. 00 

1 

5 

0 

1 

6 

3.67489999999.00 

1 

6 

0 

1 

7 

4 . 35309999999.00 

2 

7 

1 

2 

8 

4.7024     302.34 

0 

8 

2 

2 

8 

4.9501     385.58 

0 

9 

1 

2 

9 

5.41609999999.00 

1 

10 

2 

2 

9 

5.7008     415.63 

0 

11 

1 

3 

9 

4 .74309999999.00 

1 

12 

2 

3 

9 

5.0031     374.09 

0 

13 

1 

2 

10 

6. 11099999999.00 

1 

14 

2 

2 

10 

6.4318     443.29 

0 

15 

1 

3 

10 

5.46269999999.00 

1 

16 

2 

3 

10 

5.7617     403.70 

0 

17 

1 

4 

10 

4.8083     299.01 

0 

18 

2 

4 

10 

5.0436     362.25 

0 

19 

1 

2 

11 

6.78779999999.00 

3 

20 

2 

2 

11 

7 .  14369999999.00 

3 

21 

1 

3 

11 

6.16359999999.00 

1 

131 


22 

2 

3 

11 

6.5004     430.97 

0 

23 

1 

4 

11 

5.53779999999.00 

1 

24 

2 

4 

11 

5.8083     391.28 

0 

25 

1 

5 

11 

4.8860     270.43 

0 

26 

2 

5 

11 

5.0842     351.22 

0 

27 

1 

2 

12 

7.44699999999.00 

3 

28 

2 

2 

12 

7.83679999999.00 

3 

29 

1 

3 

12 

6.84609999999.00 

3 

30 

2 

3 

12 

7.21969999999.00 

3 

31 

1 

4 

12 

6.24819999999.00 

1 

32 

2 

4 

12 

6.5529     418.01 

0 

33 

1 

5 

12 

5.6271     313.96 

0 

34 

2 

5 

12 

5.8550     379.75 

0 

35 

1 

6 

12 

4.9637     278.45 

0 

36 

2 

6 

12 

5.1247     340.96 

0 

37 

1 

2 

13 

8.08899999999.00 

3 

38 

2 

2 

13 

8.51169999999.00 

3 

39 

1 

3 

13 

7.51089999999.00 

3 

40 

2 

3 

13 

7.92009999999.00 

3 

41 

1 

4 

13 

6.93999999999.00 

3 

42 

2 

4 

13 

7.27799999999.00 

3 

43 

1 

5 

13 

6.34879999999.00 

1 

44 

2 

5 

13 

6.6054     406.06 

0 

45 

1 

6 

13 

5.7164     320.43 

0 

46 

2 

6 

13 

5.9016     369.09 

0 

47 

1 

7 

13 

5.0290     269.00 

0 

48 

2 

7 

13 

5.1529     330.40 

0 

49 

1 

3 

14 

8.15829999999.00 

3 

50 

2 

3 

14 

8.60209999999.00 

3 

51 

1 

4 

14 

7.61369999999.00 

3 

52 

2 

4 

14 

7.98389999999.00 

3 

53 

1 

5 

14 

7.05159999999.00 

3 

54 

2 

5 

14 

7.33629999999.00 

3 

55 

1 

6 

14 

6.4494     359.70 

0 

56 

2 

6 

14 

6.6580     395.03 

0 

57 

1 

7 

14 

5.7915     309.51 

0 

58 

2 

7 

14 

5.9340     357.99 

0 

59 

1 

4 

15 

8.26979999999.00 

3 

60 

2 

4 

15 

8.67139999999.00 

3 

61 

1 

5 

15 

7.73599999999.00 

3 

62 

2 

5 

15 

8.04779999999.00 

3 

63 

1 

6 

15 

7.1632     396.23 

0 

64 

2 

6 

15 

7.3944     418.86 

0 

65 

1 

7 

15 

6.5340     347.41 

0 

66 

2 

7 

15 

6.6945     383.44 

0 

67 

1 

5 

16 

8.40249999999.00 

3 

68 

2 

5 

16 

8.74059999999.00 

3 

69 

1 

6 

16 

7.85829999999.00 

3 

70 

2 

6 

16 

8.11159999999.00 

3 

71 

1 

7 

16 

7.2570     382.67 

0 

72 

2 

7 

16 

7.4349     406.84 

0 

73 

1 

6 

17 

8.53509999999.00 

3 

74 

2 

6 

17 

8.80989999999.00 

3 

75 

1 

7 

17 

7.96109999999.00 

3 

76 

2 

7 

17 

8.15599999999.00 

3 

132 


77 

1 

7 

18 

8.64659999999.00 

3 

78 

2 

7 

18 

8.85799999999.00 

3 

79 

1 

8 

19 

5.0290 

255.72 

0 

80 

2 

8 

19 

5.1529 

337.32 

0 

81 

1 

8 

20 

5.7915 

296.82 

0 

82 

2 

8 

20 

5.9340 

365.87 

0 

83 

1 

9 

20 

5.0943 

263.82 

0 

84 

2 

9 

20 

5.1934 

327.16 

0 

85 

1 

8 

21 

6.5340 

335.28 

0 

86 

2 

8 

21 

6.6945 

392.24 

0 

87 

1 

9 

21 

5.8667 

303.43 

0 

88 

2 

9 

21 
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Network  of  Possible  Solutions  in  Nodal  Form 

Note  -  The  criterion  function  of  a  node  will  be  very  low  if  all  paths 
leading  to  it  fail  to  satisfy  lead  time  req.  and/or  minimum  MOE  levels, 
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Minimum  Average  Annual  Cost  Solution  Path 

(  Time  value  of  money  is  ignored,  ie.  interest  rate  is  assumed  to  be  0  ) 

Transition         Alternative   Life-Cycle   Life-Cycle     Criterion 
Node     Year     System    Employed       Cost     Effectiveness   Func. Value 
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Minimum  Annualised  Cost  Solution  Path 
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